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Abstract. This review reports on the research done during the past years on 
violations of the fluctuation-dissipation theorem (FDT) in glassy systems. It is 
focused on the existence of a quasi-fiuctuation-dissipation theorem (QFDT) in 
glassy systems and the currently supporting knowledge gained from numerical 
simulation studies. It covers a broad range of non-stationary aging and 
stationary driven systems such as structural-glasses, spin-glasses, coarsening 
systems, ferromagnetic models at criticality, trap models, models with entropy 
barriers, kinetically constrained models, sheared systems and granular media. 
The review is divided into four main parts: 1) An introductory section explaining 
basic notions related to the existence of the FDT in equilibrium and its possible 
extension to the glassy regime (QFDT), 2) A description of the basic analytical 
tools and results derived in the framework of some exactly solvable models, 3) 
A detailed report of the current evidence in favour of the QFDT and 4) A brief 
digression on the experimental evidence in its favour. This review is intended for 
inexpert readers who want to learn about the basic notions and concepts related 
to the existence of the QFDT as well as for the more expert readers who may be 
interested in more specific results. 



PACS numbers: 64.40.-i, 64.60.Cn, 75.10.Nr 
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1. Introduction 

The search for a general theory of non-equihbrium processes has been a primary goal in 
modern statistical physics. Despite of many efforts in this direction we have a limited 
understanding of the basic principles behind non-equilibrium theories. Compared with 
ensemble equilibrium theory, a general principle like the equal probability Boltzmann 
principle (that forms the basis of equilibrium statistical mechanics and provides a 
statistical foundation of thermodynamics) is still lacking. During the last century, 
the field of non-equilibrium phenomena has grown in two directions: 1) by developing 
new statistical models as an inspiring source of fruitful new concepts and ideas 2) by 
establishing partial links among different, apparently disconnected, non-equilibrium 
phenomena. 

Although much progress has been made in the first direction, the second one 
remains less unexplored. While a general principle governing non-equilibrium systems 
probably does not exist, substantial progress could be done following the second route 
in the search for basic principles governing a restricted category or class of systems. 
The applications of such basic principles may be very important because a priori many 
different systems can fall into the same category. Hence the interest in the research 
on the existence of such restricted formulations. 

During the past years it has become increasingly clear that glassy systems 
may constitute one of these large categories where their physical behavior can be 
rationalized within a restricted formulation. Glassy systems are rather common in 
nature and many systems such as structural glasses, spin glasses, disordered and 
granular materials or proteins present what is called glassy behaviors. This means 
a dramatic slowing down of relaxational processes when some control parameters 
are varied. A typical signature of glassy behaviour is a power-law or stretched 
exponential behavior of correlation functions, as opposed to exponential decay. As 
the characteristic relaxation time may change by several orders of magnitude it can 
easily exceed the observation time. As a consequence the system ages: the observed 
static and dynamic properties depend on the age of the system defined as the time 
since the system was prepared (also called waiting time) . For this reason this residual 
very slow non-equilibrium phenomena is commonly known as aging. 

Aging systems include a large variety of materials. In fact, nearly all physical 
systems, within an appropriate set of conditions and observed during a specific time- 
window, display glassy properties. The origin of glassy behaviour, however, can 
vary from system to system. The most important class of glassy systems (which 
include window glasses) are glass forming liquids where glassy behaviour is due to 
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the appearance, as some external parameter is changed, of a long-lived complex 
pattern of interacting bonds between their microscopic constituents which strongly 
inhibits relaxation toward equilibrium. Aging follows from the very-slow motion of 
such a complex pattern of interacting bonds which induces a slow change of the 
atomic structure of the liquid. For this reason glass-forming liquids are usually called 
structural glasses. Our current understanding of the slow glassy relaxation dynamics 
is greatly limited by the lack of a general non-equilibrium theory that accounts for 
these phenomena. 

Glasses can be generated by the fast cooling of a liquid. Upon cooling from 
high temperatures down to the melting transition temperature Tm, sometimes 
crystallization does not occur and the liquid continues its way down in temperature 
beyond Tm by following a line (called supercooled liquid line) which is the continuation 
of the liquid line. As the liquid line is thermodynamically stable only above Tm, 
the supercooled liquid line is metastable with locally equilibrated properties, so its 
lifetime can be extremely large. As cooling proceeds it is observed that the supercooled 
liquid falls out of equilibrium (i.e. departs from the supercooled liquid line), below a 
temperature T* (r) which depends on the cooling rate r. The state reached below T* (r) 
is called a glass and the corresponding relaxational regime is indistinctly termed as 
aging or glassy. For small values of r a sharp transition is observed at T*(r), usually 
referred as structural arrest, where the heat capacity jumps down, indicating the 
freezing of degrees of freedom. Contrarily to the supercooled state, the glass state is 
of non-equilibrium nature and T*{r) is observed to decrease with r. As T*{r) depends 
on the cooling rate, no equilibrium phase transition occurs at that temperature. This 
means that the liquid will eventually equilibrate back to the supercooled state. The 
equilibration process may take an extremely long time (even for temperatures only 
a few degrees below T*) being inaccessible from any practical point of view. Under 
some conditions the equilibration time can be larger than the age of the universe!. In 
these conditions the glass state is the only observable state. 

Long equilibration times imply that the glass state is characterized by very low 
energy dissipation rates, also called entropy production. This may give the false 
impression that the glass is in a stationary state. For instance, a piece of silica glass 
at room temperature looks pretty stable, indeed its optical, electrical and mechanical 
properties appear constant in time. However a more careful examination reveals that 
the physical properties are constant only if observed on time-scales much smaller than 
the time elapsed since the glass was prepared or formed. Beyond that timcscale, the 
physical properties change revealing that the glass is aging. 

Although aging was identified long time ago in the study polymers X. it has 
received a renewed interest in connection with the study of spin glasses. Measurements 
of the magnetization in spin glasses have shown that aging is a general property 
of the low-temperature spin-glass phase. There are several types of spin glass 
materials, the most common ones are metallic spin glasses. These are random 
diluted magnetic systems where glassy behaviour arises from the disordered pattern of 
exchange interactions, rather than being self-generated as in structural glasses. Indeed 
random dilution generates exchange interactions with random competing signs, the 
system is then frustrated since a finite fraction of bonds cannot be satisfied. Aging 
is consequence of the slow evolution of the pattern of satisfied bonds which becomes 
strongly inhibited as the temperature is lowered. 

Another class of systems with glassy properties are driven systems which, under 
certain conditions, reach a stationary state characterized by non-Gibbsian probability 
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distributions. After applying a time-dependent perturbation of frequency w > l/t^q, 
upon an initially equilibrated system of relaxation time icq, a new stationary state is 
reached which for many aspects is similar to the aging state of the relaxational system 
of age ~ 

Another important aspect of glassy systems that has received considerable 
attention for a long time |2I is the idea about the existence of an effective temperature 
(sometimes also called fictive temperature) describing the non-equilibrium properties 
of the glassy state. During the last decades it has emerged that a possible way to 
rationalize the existence of an effective temperature is by measuring violations of the 
fluctuation-dissipation theorem (FDT). In glassy systems, a new modified relation 
between correlations and responses that goes under the name of quasi FDT (QFDT) 
provides a description of the dynamics in the glassy state by quantifying the violations 
of the FDT. In this new theorem the effective temperature plays the role of the 
temperature of the bath. Related to the concept of the effective temperature is the 
idea of the existence of a heat flow from the glass to the thermal bath put in contact 
with the system. As the glass has an effective temperature larger than that of the 
bath, the heat flows from the glass to the bath. However, the energy dissipation rate 
from the glass toward the bath is extremely low (hardly measurable) and, in general, 
this flow can be also understood in terms of an effective very low thermal conductivity. 
The reader should be aware that using QFDT is only one among other possible ways of 
introducing an effective temperature for the description of the glassy state. In general 
other definitions which use a generalization of different equilibrium relations to the 
non-equilibrium regime are possible. This arises the problem of the equivalence of all 
possible definitions. We shall not discuss this point in this review and we will stick to 
the QFDT definition of an effective temperature. 

This review will concentrate on the existence of a QFDT, its physical meaning, 
in what conditions it can emerge and the numerical evidence reported in favour of its 
validity. This is a rapidly growing area of research which is attracting new condensed 
matter and statistical physicists. We will report here the most important results 
obtained until the summer of the year 2002. Although we have tried to cover most 
of the published work some contributions may have been overlooked. We apologize in 
advance to those colleagues. Although some of the results here reviewed are currently 
well understood many others still lack a full comprehension so it is not exaggerated to 
say that some of the ideas and suggestions described in this review could be modified in 
the future to adapt to the forthcoming theoretical, numerical or experimental evidence. 
Most of the results here reported deal with relaxational aging systems (as compared 
to driven systems) since these are the ones that have mostly attracted the attention 
of the researchers in the field. However, future developments in this exciting area 
of research might compensate this original unbalance as driven systems appear more 
amenable to experimental research than aging systems. Moreover, in this review 
we shall only consider the FDT in its classical version. Although most of the the 
ideas can be extended to the quantum regime, to our knowledge there are neither 
numerical or experimental works challenging FDT violations in the quantum aging 
regime. Therefore we shall not address them, the interested reader is referred to a 
recent review Pj. 

Many textbooks and article reviews can be useful to complement the contents 
of this review. Basic reviews on the glass transition since mid 80 's until now can be 
found in Other accounts deahng with aspects of the glass transition include: 

thermodynamic theories of the glass transition El , mode-coupling theory (MCT) 
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| 1UI [TT] and numerical simulations 12 . For spin glasses a good selection of review 
articles can be found in ^^l- A clear discussion of mode-coupling approximations 
in the context of disordered systems can be found in jJJ. A recent discussion of 
several aspects of concerning FDT violations can be found also in a recent review on 
kinetically constrained models ^B]. Finally, a throughout compendium of analytical 
methods for glassy dynamics has been recently collected in fjl and a review of granular 
systems in Proceeding articles covering several aspects of glasses and spin-glasses 
can be found in |17l IT^ and for kinetically constrained models in jTHI . 

The contents of this review have been written with two kinds of readers in mind: 
inexpert and expert. Those inexpert readers who want to understand the most basic 
ideas as well as the interest of investigating FDT violations must read Sections I2I3I4I 
These sections have been written at an introductory level, so expert readers who know 
about the subject may directly start reading from Section|3 However, a careful reading 
of section 0] is recommended to those readers who want to have a more physically 
appealing description of the possible origin of FDT violations. Section |S1 deals with 
some of the thermodynamic consequences of FDT violations. Sections I6I7I constitute 
the core of the review. Section El describes our knowledge of FDT violations gathered 
from several exactly solvable models where many aspects of their non-equilibrium 
behavior can be understood by analytical means. Sectiondcovers all evidence collected 
in the past years in favour of the existence of a QFDT in glassy systems. Many of 
the model systems described in this section correspond to realistic as well as model 
systems for which analytical solutions are hardly known. The expert reader who wants 
to grasp the state of the art concerning these questions will be mainly interested on 
these two sections. Finally, a brief account of some experimental results on FDT 
violations is described in Section |S1 Section |51 presents some conclusions. 

2. Basic definitions and concepts 

In this section we recall some concepts of equilibrium theory which will be needed 
later for the description the glassy state. 

2.1. The microcanonical and canonical ensembles 

The foundations of equilibrium statistical mechanics rely on the maximum entropy 
postulate and the Boltzmann's equal probability hypothesis. An introduction to the 
basic postulates can be found in the classical books on statistical mechanics, rather 
excellent are that by Ma |^ and Callen [21]. Good discussions also comes from the 
Information Theory, see for example the book of Beck and Schlocgel "22". 

In what follows we shall denote with C a generic system configuration in the 
phase space. The phase space can be either continuous or discrete depending on the 
particular system. For example, for a system of N particles C are the positions and 
momenta in a continuous GA^-dimensional space, while for a system of N 1/2-spins C 
is a point in a a discrete N dimensional space with 2^ points. The system evolves 
in time following a dynamical rule which generally speaking is a rule that for each 
configuration C associates a new configuration C . The set of configurations which can 
be visited given a dynamical rule, defines the region of motion in the phase space. 
Let r be the volume of the region of motion allowed by the invariant quantities. The 
basic assumption of statistical mechanics asserts that the entropy is the logarithm of 
r. This makes the entropy computable without having to solve the dynamics. If we 
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assume that the system is described by an energy function E(C) then the motion is 
confined to a region in phase space of constant energy. The calculation of entropy is 
then reduced to 

S{E)=lnY,S{E-E{C)) (1) 
c 

where in the case of continuous variables the sum must be read as an integral. This 
equation defines the microcanonical ensemble. Since all allowed states are included 
S{E) is clearly a maximum over all possible regions of constant energy E into which 
the phase space can be divided. 

Let us consider now an observable A{C) which we will assume to be neither a 
constant of motion or a univocal function of the energy E. We also assume that 
A{C) is extensive, i.e., it is proportional to the system size (volume or the number of 
constituents). We can then divide the phase space according to the value of A{C) and, 
defining the degeneration il.{E, A) of the partition as the total number of configurations 
C of energy E and observable value A: 

f](£;,A)-^<5(£;-£;(c))<5(A-A(c)), (2) 

c 

introduce the entropy in analogy with : 

S{E,A) ^liin{E,A). (3) 
Using the integral representation of delta functions fl{E,A) can be rewritten as 

exp {laiE + ia2A) ^ exp [-iaiE{C) - ia2A{C)] 
-00 (.^"'j 1^ 

^^l^e^p[S{E,A,a,,a2)] (4) 

where the function S{E, A, ai, 02) is given by, 

S{E, A, ai, a2) — iaiE + ia2A + lnZ(ai, 02) (5) 
and Z{ai, 0,2) is the partition function given by, 

Z{ai, a2) — exp [— iQ;i£'(C) — ia2A{C)\ (6) 
c 

Since both energy and the observable are extensive quantities the sum in (0 is 
dominated, in the limit of large system size, by the largest contribution, and Z 
is exponentially large in the system size. The function iS(£', A, ai, 012) is then an 
extensive quantity and, in that limit, the integrations can be done selecting the 
dominant contribution using the saddle point method: 

dS{E,A,ai,a2) dS{E, A,ai,a2) 



dai da2 
which leads to the saddle point equations 



(7) 



E = V E{C) exp [-iaiE{C) - ia2A{C)] = {E) (8) 



Z{ai,a2) ^ 
1 



A - V A{C) exp [-iaiE{C) - ia2A{C)\ = (A) (9) 

Z(ai,a2) ^ 



Violation of FDT in glassy systems 



8 



Reality of i^{E,A) implies that the solution al{E,A), al{E,A) of the saddle point 
equations must be pure imaginary: = —ifi, a2 = — i/i with /3 and fi real numbers. 
The entropy Q is given by the value of S{E, A, ai, evaluated at the saddle point: 

S{E, A) = S{E, A, -i(3, -i^i) ^ l3E + fiA + lnZ{f3, /i). (10) 

We can now ask the following question. What is best choice of the value of A{C) for 
which S attains its maximum? Stationarity of S with respect to ai and a2 at the 
saddle point implies 



dS{E, A) 

dE 
dSiE, A) 



a (11) 

dA ^ . . ^''^ 

so that the maximum entropy assumption requires = 0, i.e., the entropy S(E^A) 
must be stationary with respect to variations of A. For any energy E the best choice 
of A is then: 

1 



2(/3)y 

where 

Z(/3) = = exp [-/3F(/3)] (14) 

c 

The value of the entropy is 

S(E) ^ iiE + \^Z(,0) (15) 

and is independent on A as required from stationarity. Finally from l|ll|) it follows that 
can be identified with the temperature T of equilibrium thermodynamics, while 
insertion of ((1^ into ifTCjl yields the thermodynamic relation F(/3) — E{T) — TS{E) 
identifying F{(3) with the Helmholtz free energy. Equations Hill) . (|13|l . (|14|l and (|15|l 
define the canonical ensemble. 

At difference with the microcanonical ensemble the measure of the canonical 
ensemble is not restricted on states of constant energy. All possible states C enter 
but with a weight proportional to exp[— /3i?(C)]. For a given temperature, however, 
only states with energy E = (E) given by ((T^ for A{C) = E{C) [See significantly 
contribute to the measure. The temperature can be seen as a Lagrange multiplier 
used to fix the value of the energy. Vice versa each value of E in equilibrium selects 
a temperature T through Hll|l . 



2. 2. Einstein fluctuation theory 

A key contribution in the development of equilibrium statistical mechanics is the 
statistical theory of fluctuations developed by Einstein [21]. In the previous section 
we have seen that in equilibrium the value of any observable A{C) is given by H13(l . 
which for the purpose of this section will be denoted by A^q. The equilibrium value 
corresponds to the most probable value of A{C), i.e., to the value of A{C) which has the 
overwhelming probability to be seen in equilibrium. The same considerations apply 
to the energy E{C) in the canonical ensemble. 

We may then ask what is the probability of observing a value of A{C) different 
from the equilibrium value. This probability is simply proportional to the number of 
configurations with A{C) = A which from is 

Pi6A) cx n{E,q, A = Aeq + SA) cx exp [S{E,q, A.^ + 6A)] (16) 
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where foUowing our notation E^q is equilibrium energy. For sniaU value of the 
fluctuations 6 A the exponent can be expanded and using stationarity of the entropy 
with respect to variation of A we get, 

SiE,q,A)^SiE,q,A,q) + ^f^\ +0{SAy (17) 



If Aeq is a maximum, then a necessary condition is 



\^A2} =-7^<0- (18) 

\dA^JA=A^^ TXA 



From \i<6\ we finally obtain. 



P(M) = -=^cxp 



{5Af 



(19) 



As xa is an extensive quantity only subextensive fluctuations SA ^ VV, where V is 
the system size, have finite probability in equilibrium. This justifies the most probable 
character of the equilibrium value Acq. 

The quantity XA is called susceptibility, and from l|19() is related to fluctuations 
of A through 

TxA = {A') - {Af. (20) 

This relation is the simplest form of the static Fluctuation-Dissipation Theorem (FDT) 
which relates the magnitude of thermal fluctuations with the response of the system 
to a (small) perturbation. 

Suppose we add a constant perturbation ~eA{C) to the energy E{C). Then in 
the new equilibrium value (A)^ of A is [See (|13f) ] 

^ Ec m exp[^/3£;(C)+/3eA(C)] 

Ecexp[-/3i?(C)+/3eA(C)] ' ^ ' 

The susceptibility XA is defined as the variation of (A) induced by a small perturbation 
d{A), 



XA 



de 



(22) 

=0 

and hence measures the response of the system to the perturbation. Inserting (|21|l 
into (|22|l a straightforward calculation leads to (|20|l . 

The FDT formula (|20|) is non-trivial result, since it relates different physical 
processes: the susceptibility describes an extensive 0{V) variation of the observable 
A while the r.h.s. of (|20|l describes subextensive 0{VV) thermal fluctuations. This 
fact is at the basis of Onsager regression principle discussed in the next section. 



2.3. The Onsager regression principle: a simple derivation of the 
fluctuation-dissipation theorem (FDT) 

Onsager proposed [23 1^ a simple derivation of FDT for time-dependent 
perturbations. The derivation bypasses the more cumbersome analytical developments 
using linear response theory formalism, the Fokker-Planck equation or the generalized 
master-equation approach. 

Onsager derivation is based on the following regression principle: If a system 
initially in an equilibrium state 1 is driven by an external perturbation to a different 
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<A(t)> 



e(t) 



t=0 t 

Figure 1. Perturbation e{t) and typical evolution curve for {A{t))c 



equilibrium state 2, then the evolution of the system from state 1 toward state 2 in the 
presence of the perturbation can be treated as a spontaneous equilibrium fluctuation 
(in the presence of the perturbation) from the (now) non-equilibrium state 1 to the 
(now) equilibrium state 2. 

Suppose that the system is initially in equilibrium with a thermal bath at 
temperature T, then the probability distribution of system configuration C in state 1 
is given by the canonical ensemble H13|l: 



PoiC) 



exp 



(23) 



The subscript 
At time t 



exp [-mC)] 
0' indicates that the system in unperturbed. 

= a constant perturbation coupled to the observable B{C) is applied 
to the system changing its energy into 

E,{C)^E{C)-e{t)B{C), (24) 

where e(t) = e if i > 0, and zero otherwise. The effect of the perturbation can 
be monitored by looking at the evolution of the expectation value {A{t))^ of an 
observable A{C), not necessarily equal to B(C), from the equilibrium value in state 1 
{A{t = 0))e = {A)o toward the new equilibrium value in state 2. The shape of e(t) 
and a typical evolution of {A{t))e are shown in Figure H 

The expectation value of (A(i))e is given by the average over all possible dynamical 
paths originating from initial configurations weighted with the probability distribution 



{A{t)), = J2 A{C)P,{C,t\Co,0)Po{Co) 

C,Co 



(25) 



where P^{C,t\Co,0) is the conditional probability for the evolution from the 
configuration Co at time t — to the configuration C at time t. If e = the 
expectation value becomes time independent since the initial state is in equilibrium 
and Po(C,i|Co,0) describes spontaneous equilibrium fluctuations. 

The Onsager's regression principle asserts that the conditional probabilities 
after having applied the perturbation are equal to those of spontaneous equilibrium 
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fluctuations in state 2. Hence since the state 2 is still described by the canonical 
ensemble (|23f) . but with the energy H24|) now including the perturbation term, then 

P,{C,t\Co,0) = Po{C,t\Co,0) exp{/?e[B(C)-B(Co]}. (26) 

where the r.h.s of this equations is just the product of the spontaneous equilibrium 
fluctuations conditional probabilities Pa{C, t\Co, 0) in state 1 corrected by the presence 
of the perturbation term eB{C) %. Inserting into and expanding the 

exponential up to the linear order we get, 

{A{t)), ~ (A)o = /36^ A{C) [B{C) ~ B(Co)] Po(Co) 

C,Co 

= /3e[(A(i)SW)o-(A(t)i3(0))o]. (27) 
If we define the correlation function and time-dependent susceptibility as 

CAs(i,s) = (A(t)B(s))o (28) 

e-»0 e 

then from l|27|) we get the integrated form of the FDT relation 

XA,i3(0 =/9[C^AB(i,<)-CA,B(i,0)] (30) 

The static form of the FDT 12U|I is easily obtained from (|30|l by taking A ^ B and 
the limit t oo. In this case XA,A{t) — > Xa [dr. (|22|) ] and CAA{t,t) — {A'^)o while 
CAA{t,0) {A)q as correlations factorize for infinitely separated times. 

Sometimes the FDT relation is written in a differential form by considering the 
two-times response or retarded Green function: 

i?As(i,5) = ^^^, t>s (31) 
Se{s) 

which gives the response to an impulsive perturbation e(s) acting at time s. Causality 
imposes that the response function RA,B{t, s) is zero for t < s: perturbations cannot 
propagate backward in time. The susceptibility XA.nit) is the integral of the response 
function Rji B{t,s) then, using (|30|l : 

RA,B{t,s)ds = XA,B{t) 

^ f3[CA,B{0) - CA^Bit)] 

= 13 ^^CA.B{t,s)ds. (32) 
The last equality, and the arbitrarily of time t, implies that 

RA.B{t,s)= P-^CA.B{t,sW-s) (33) 
OS 

This is the differential form of the FDT relation. 

I The relation 1261 is only valid for Co ^ C. Indeed, the normalization condition of conditional 
probabilities J]^ Pj(C, t|Co, 0) = 1 implies that P^{C,t\C,0) + Po{C,t|C,0) so relation does not 
hold for Co = C. The difference between both probabilities Pe(C, t|C, 0), Po(C, 4|C, 0) does not matter 
as the transition C — > C does not contribute to the response function in 1271 
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3. The Master Equation Approach 

In this section we introduce the Master equation for the dynamical evolution of a 
generic system and show how the FDT arises within this approach. Besides the 
previous derivation in Section 12.31 many other derivations of the FDT exist. We 
give a few collection of references where these are presented. Derivations can be 
classified in two families: deterministic or stochastic. Deterministic approaches are 
linear response theory |25| , operator formalism for master equations '26' and quantum 
statistical mechanics ^7\. Stochastic approaches are the Langevin and Fokker-Planck 
equations j28| . Here we present a stochastic derivation which is convenient for the 
purpose of the present review. 



3.1. The Master Equation (ME) 

Any dynamical law describing the evolution of a system is a rule which for each 
system configuration C associates a new configuration C . The time is just a label for 
bookkeeping the sequence so generated. Therefore to simplify the presentation and 
the notation we shall consider the time as an integer variable giving the equivalent 
expressions for the continuous limit when needed. This picture has also the advantage 
of being more closely related to numerical simulations since all numerical methods use 
discrete time schemes. 

The dynamics can be encoded into the conditional or transition probability 
W{C,t\C' ,t — 1) of going from configuration C at time i — 1 to configuration C at 
time t. Indeed if P{C,t) denotes the probability that the system at time t is in the 
configuration C, then from the Bayes theorem it follows that 

P{C, t) = Y, * - 1) ^ - !)■ (34) 

c 

and the W{C,t\C' ,t — 1) together with the initial condition P{C,0) fully define the 
dynamical evolution of the system in the phase space. Equation H34I) is an identity 
valid for all processes and is the first of a hierarchy of equations for joint probabilities. 
Only if the process is Markovian, i.e., only if the conditional probability is determined 
entirely by the knowledge of the most recent past, then the hierarchy can be closed. 
Equation (|34|l and probability conservation at all times {J^c ^i^'^ ^) ^ 1) implies that 
W{C, t\C' , t — 1) must satisfy the normalization condition 

W{C', t\C, i - 1) = 1 for all C and t. (35) 

c 

In the continuous time limit (|34|l is not well defined. To have an expression valid 
in this limit one then considers the variation of P(C,t) between two successive times 
which, using the normalization condition ()35|1 . reads 

P{C,t+l)-P{C,t) = 

W{C,t+l\C',t)P{C',t) - J" W{C',t + l\C,t)P{C,t). (36) 



c c 



Dividing both sides of this equality for the time increment At and and taking it to 
zero we get the Master Equation (ME) 



dP{C,t) 

at 



J2 W{C\C'; t) P{C', t)^J2 ^(^'1^; ^) ^(^' ^) (37) 
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where W{C'\C;t) — \im^t->o ^ , t + At\C' ,t)/ At is called transition rate and gives 
the transition probability per unit of time. Solving the ME is often an extremely 
difficult and unfordable task, even in the Markovian case. 

Loosely speaking the transition probability W{C, t\C' , t — I) can be seen as a 
transition rate for a unit time interval {At — 1), thus in what follows we shall not 
make distinction between transition probabilities and transition rates and shall call 
them generically transition rates using for both the notation W{C'\C; t). Which one is 
the appropriate will be clear from the context. 

Transition rates depend on the specific dynamical rules and hence by the 
Hamiltonian and eventual constraints (holonomic or non-holonomic) . Let the system 
under consideration be described by an Hamiltonian which can depend on time 
through a set of time-dependent external parameters Aj. For instance, Xt may denote 
an time-dependent external pressure applied to a liquid or a time varying electric or 
magnetic field applied in a dielectric or a magnetic medium. We shall denote the set of 
these parameters by the vector At and the Hamiltonian by Ti.x^ (C) to indicate the time 
dependence through At. Accordingly we shall also denote W{C\C',t) by W\f{C\C'). If 
the Hamiltonian is time-independent either At or just the subindex t will be dropped, 
depending on the context. 

Regardless of their form the transition rates must satisfy the following 
requirements: 

• Non-negativeness and normalization. The Wx^{C\C';t) are probabilities so they 
must be non-negative and satisfy the normalization condition 1)35(1 . 

• Ergodicity. Transition rates must be such that starting from any configuration 
any other configuration of a finite system can be visited in a finite time. For 
continuous variables the condition is stated by considering an arbitrary finite 
phase space region around a given point ('neighborhood'). 

• Detailed balance. If the A are time-independent the equilibrium distribution 
P^'^{C) is the stationary solution of the ME (|37|l . A sufficient condition for this 
is that the transition rates be time-independent and 

wx{c'\c) ^ pnc) 

Wx{C\C') P^^(C) ^ ^ 
This condition receives the name of detailed balance. Different equilibrium 
ensembles are thus encoded into the different forms of transition rates. The 
Perron-Frobenius theorem assures that this condition, together with non- 
negativeness, normalization and ergodicity guarantees that the equilibrium 
distribution in a finite system is reached in a finite time. 

• Causality. This is an important assumption for the time-dependent transition 
rates and means that future is only determined by the past, i.e., a perturbation 
applied at a given time can only propagate forward in time and not backward. 
Consequence of this is that the transition rates must depend only on the values 
of At taken at the lowest time t. 

For arbitrary time-dependent At the stationary solution of ME will in general 
also depend on time. In this case, however, the previous conditions, and in particular 
the detailed balance condition, are not enough for determining the stationary state 
which in general will not be Gibbsian, i.e., not described by the Boltzmann-Gibbs 
distribution. Only when At can be treated as small perturbations some predictions 
can be obtained from the linear response theory. 
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3.2. Correlations, responses and the FDT 



Consider two arbitrary observables A{C), B{C) which for simplicity are assumed time- 
independent. These can be either local or global quantities defined over a microscopic 
or macroscopic region respectively. For instance, in a liquid an observable could be 
the local density at a point or the total mass of a given macroscopic region. In systems 
with discrete variables such as magnetic systems, it can be a spin of a given magnetic 
atom or the magnetization of a macroscopic part of the system. 

The two-times correlation function CA,B{t,s) between A{t) and B{s) is defined 
as the average {A{t)B{s)) over all possible dynamical paths from time to time t and 
all possible initial conditions weighted by the probability distribution P(C, 0), 



CAAt^s) - {A{t)B{s)) 



E 

Cfi' 



A{C')P{C',t\C,s)B{C)PiC,s) (39) 



where P{C',t\C, s) is the conditional probability to evolve from C at time s to C at 
later time t. Unless otherwise stated, in what follows we shall adopt the convention 
that t > s. 

To simplify the notation we switch to discrete (integer) time variable so that 
a dynamical path from time to time t is given by a sequence of t configurations 
{Co,Ci, . . . ,Ct}. along which A takes the sequence of values {Aq, Ai, . . . , A(}. In this 
case, using H34II , the correlation can be easily rewritten in term of transition rates as 



CA,B{t,s) 



Cs,...,Ci 



A{Ct) 



Yl WxACk+i\Ck) 



_k—s 



B{Cs)P{Cs,s) (40) 



where we have used the short-hand notation W\^{Ck+i\Ck) = Wx^{Ck+i\Ck; k). 

In equilibrium correlations satisfy the time-translational invariance (TTI) 
property: CA,B{t, s) = CA,B{t — s). Indeed in this case P{Cs, s) is replaced by P'^'i{Cs) 
and the transition rates satisfy the detailed balance condition ()38|): 



Wx{Ck+i\Ck) ^ Wx{Ck\C, 



k+l) 



k+l, 



(41) 



Inserting this relation into (|40|) the factors P'^'^{Ck) in the numerator and denominator 
of the product cancel one by one and, exchanging the indexes i <-!■ s, we obtain 

"t-i 

l[Wx{Ck+i\Ck) 



CA^Bit,s) = 



E 



..,Ct 

{t ^ s) 



.k—s 



B{Cs) 



E B{C, 



Ct,-fis 



\{Wx{Ck\Ck^ 



.k=t 



A{Cs)Pl\Cs) 



(42) 



= {B{t)A{s)) 

= CB,Ait,s). 

which implies that CA,B(t, s) = CA,B{t — s). 

The correlation function CA,B{t, s) is a measure of how the system loses memory of 
its previous past history and hence decays for large time separations. To measure how 
a system responds to external perturbations one introduces the response functions. 
Similar to correlations, responses also tend to decay with time because the effect of 
the perturbation is progressively forgotten in a thermal environment. However there 
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is an important difference between correlations and responses: causality. While two 
observables can be correlated forward or backward in time, a perturbation cannot 
propagate backward in time and the response of the system for times before the 
perturbation is applied must be zero. Nevertheless despite this difference response 
and correlations can be treated on equal footing by employing a supersymmetric 
formalism. The interested reader can find more details, e.g. in the classical book 
by Zinn- Justin j28| . 

To study the response of the system to an external perturbation we assume that at 
time s an impulsive perturbation of small intensity e is applied to the observable B{C) 
and measure the variation of the average value of an observable A{C) at later times. 
The response function RA.Bit,s) is defined in the limit of vanishing perturbation 
strength as 

(A(O)e. - {A{t))o 



RA.sit, s) 



lim 



t > s 



with 



{A{t)),^=Y,A{C)P,^{C,t), 
c 

(A(t))o-5] A(C)Fo(C,t). 



(43) 

(44) 
(45) 



where Pe,(C, t) and Po{C, t) are the probability that the system is in the configuration 
C at time t > s in the perturbed and unperturbed case, respectively. If the system is 
described by the unperturbed Hamiltonian 7io(C), then in presence of the perturbation 
the Hamiltonian becomes: 

n,XC)=nn{C)-5t,stB{C). (46) 
By using 1)34(1 and the fact that the Hamiltonians only differ at time s when the impulse 
is applied, the probabilities P^^{C,t) and Po{C,t) can be written as: 
t-i 



P.ACut) 



Po{Ct,t) 



E 



W Wo(Cfc+i|Cfc) 



E 



.k=s+l 
t-l 



W,{C,+i\Cs)Po{Cs,s) 



(47) 



n Wo{Ck+i\Ck) 



^ Wo(C,+i|a)Fo(C.,s) (48) 

where Wo,£ denotes the transition rates in the unperturbed/perturbed case, and we 
have used the short-hand notation W{Ck+i\Ck) = W{Ck+i\Ck] k). 

Because we are interested in the e — > limit the transition rates in the perturbed 
and unperturbed case can be related by expanding the detailed balance condition ((38|l 
for the perturbed state around e = up to the first order in e, 
W,{C'\C) P^%C') 



W,{C\C') 



P!\C) 
^ W^{C'\C) 
Wo{C\C') 



1 



^ 1 

e — In 



P^\C') / P°,'^{C) 



(49) 



where Pf'^{C) is the equilibrium probability distribution in the perturbed state f . In 
general we can write, 

\og[P-,\C)] = log[P°^(C)] + 00 - 0. + peB{C) (50) 

t An observation concerning 1491 and its relation with the Onsager postulate 1261 is important. 
Ideally, in order to demonstrate FDT, one would like to have a relation similar to 1261 relating 
the unperturbed and the perturbed rates rather than a relation between the forward and backward 
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where (p denotes the corresponding thermodynamic potential. For instance, in the 
canonical ensemble it corresponds to minus the Helhmoltz free energy F while in the 
grandcanonical ensemble it corresponds to the grandcanonical potential given by the 
pressure times the volume. 

Using H49() we finally obtain to the leading order in e, 

'W,{C\C') 



W,{C'\C) - Wo{C'\C) = 



- 1 



Wo{C'\C) 



+ e 



Wo{C\C') 

P°^(C')/P°'^(C) 



^ 1 



Wo{C'\C) (51) 



c=0 



The physical meaning of the two terms appearing in this expression is different. 
The first term, absent in the Onsager's postulate H26(l . accounts for the variation 
due to the change in the transition rates and does not directly depend upon the 
particular form of the equilibrium distribution. The second term depends directly on 
the equilibrium distribution. This distinction is important since they give different 
contributions to the response. For the response function H43|) to be well defined the 
difference {A{t))f^ — (yl(t))o must be at least linear in e. This requirement is at the 
roots of the applicability of linear response theory and implies that both terms in (|51|l 
must be at least linear in e. Concerning the first term, it is required that the difference 
We(C|C')/tyo(C|C') — 1 to be linear in e, so the transition rates must change linearly 
with e when the system is perturbed. This means that, during the dynamics, and after 
applying the perturbation, one cannot switch arbitrarily from one class of dynamics to 
another class of dynamics in a random fashion. To better illustrate what this means let 
us consider a Monte Carlo stochastic dynamics. The are different possible algorithms 
or transition rules which fulfill detailed balance (e.g. Metropolis, heat-bath, Glauber). 
One could imagine of switching randomly from one algorithm to another one while 
doing the dynamics. Nothing forbids this quite artificial choice. But, when measuring 
the response function, it is required that the same time sequence of algorithms must 
be used for the unperturbed and perturbed dynamical evolutions. Usually, the same 
algorithm is chosen in a given simulation so one does not care about this subtlety. For 
the second term, we require that (d/de) \n[P°'^{C') / P^^{C)]\f=o be finite. Inspection of 
H50|l reveals that this holds if the observable B{C) does not jump discontinuously when 
the perturbation is switched on. This condition requires that the system is not at a 
first order transition point. This situation is encountered, for instance, by perturbing 
the Ising model at zero temperature with a uniform magnetic field. 

Inserting (|51|) into (|43|l the response function decomposes into two parts. 



i?^,s(i, s) - R^aU*, s) + R%{t, s) 



(52) 



where. 



i?^^^ (t 



1 



s I = lim — > 



A{Ct) 



t-i 

n 



WoiCk+i\Ck) 

.k=s+l 

rates as given in 1491 . However, such a relation does not exist as it depends upon the type of 
dynamics through the particular form of the transition rules. For instance, by considering rates of 
the type WxiC'\C) oc P^(C'), i.e. depending only upon the final configuration, one finds that i26l 
automatically holds. However such rates cannot be used to derive the FDT 1331 . 1551 as they lead to 
the trivial identity = because any impulse does not affect dynamics at later times R{t, s) = and 
there are no time correlations. In other words, the Onsager postulate extended to the non-equilibrium 
regime generates violation terms different to those obtained in the present approach (for instance the 
term in the r.h.s of IS^I. 
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(53) 



II Wo{Ck+i\Ck) 

.k=s+l 



P^''{C')/P°,\C) 



Wq{Cs+i\Cs)Po{Cs,s) 



(54) 



Inserting (|50|l for the second term we get 



Rfsii^s)^ {A{t) 



d\nP^^%Cs+i 



de 
A{t) 



£=0/ 

d\nP^'i{Cs 



de 



e=0/ 



which using (|5()(l becomes 

R%{t, s) = l3{A{t)B{s + l))o - (3{A{t)B{s))o 
In the limit of continuous time R^'^^ in 156(1 must be replaced by 



Rfj,{t,s)=p^^{A{t)B{s))o 



d 



(55) 
(56) 
(57) 



The first term i?*^^-* cannot be expressed in a simple form. Only in equilibrium it is 
possible to show that it vanishes. To prove it requires the following steps: first use 
the identity 



t-i 



W T^o(Cfe+i|Cfe) = 



fc=s+l 



t-1 



n wo{Ck\Ck+i) 



_k=s+l 



(58) 



and then the normalization 1(351) . Collecting all terms we finally obtain for the response 
function 



i?^,s(i, s) = i^^^^i, s) + {A{t)B{s))^ 



(59) 



In equilibrium, besides i?*^^^ = 0, the correlation function satisfy TTI, so from 
1(59(1 we get the equilibrium FDT 



t > s (60) 



The term R^^^ may also vanish in the non-equilibrium state if the first term in ((51(1 
vanishes faster than linearly with e for e ^ 0. In this case ((59(1 reduces to the usual 
FDT relation (jSHJ- 



RA.B{t,s)^ PS{t-s)^CA,Bit,s) 
OS 



(61) 



Therefore, the lowest time s has a special role in the relation 1(59(1 between 
correlations and responses. This is not a surprise since the relation has been obtained 
by assuming causality which privileges the lowest time. In equilibrium, the role of the 
lowest time disappears because the system is TTI. 
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3.3. The Component Master Equation 

Let us now divide the phase space into different non-overlapping subsets TZ that can 
be called regions, phases, components or domains. In what follows, if not stated 
otherwise, we shall use the term component. Later in Section 14.21 we will see that the 
reduction of the phase space by a suitable partitioning can be relevant for the study of 
the non-equilibrium regime in glassy systems and in particular for the understanding 
of the FDT relations. For the moment, however, we do not attach any physical 
meaning to such a partitioning, and assume it to be completely arbitrary postponing 
the identification of a suitable partitioning scheme for glassy systems to Section E31 
For each partition of the phase space the probability V{Tl,t) that system be in 
the component TZ at time t is given by, 

cell 

Y,V{TZ,t)^l (63) 
n 

where the normalization condition follows from normalization of P{C,t) and the 
non-overlapping assumption on TZ. The probability distribution T^{TZ,t) obeys the 
master equation obtained by projecting the (microscopic) master equation H37() over 
the component space. The Markovian character of the dynamics is preserved under 
projection. Summing H37(l over all configurations C belonging to a given component 
TZ we get the Component Master Equation: 

^^gl^ = W(7^|7^'; t) P{TZ',t) - W(7^'|7^; t) r{TZ, t) (64) 
w w 
where yV(TZ'\TZ;t) are the component transition rates which in terms of the original 
transition rates W{C'\C;t) read, 

>V(7^'|7^;t)^ (65) 

The component transition rates satisfies the same normalization conditions as W . For 
example in the case of discrete time we have [Cfr. 

^ W(7^',^|7^,^- 1) = 1 for aU7^ audi. (66) 
w 

However there is an important difference between the (microscopic) master equation 
H37|l and component master equation H64|l : the transition rates H65fl are time-dependent 
even in for time-independent Hamiltonians since are computed with the run-time 
configuration probability distribution function P{C,t). Direct consequence of this is 
that, while the properties of non-negativeness, normalization, ergodicity and causality 
discussed in section lXTl for W do apply to W, the detailed balance is not valid anymore 
in the component space. Indeed from (|65|l it follows that 

WiTZ'\TZ;t) ^ j:ceu.c'^Ti'W{C'\C)P{C,t) ^ v(n' ,t) 

W{TZ\TZ'-t) Y.cenfi'en'W{C\C')P{C',t) V{TZ,t) ^ > 

The detailed balance condition is recovered, however, at equilibrium where W, 
accordintg to l|65|) . becomes time- independent: 

.,p^,M^(C'|C)P°q(C) 

W•=^(7^'|7^) = ^ceR.c en \\ J 

V I J Vi{TZ) 
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where T"^'^{TZ) denotes the equihbrium probabiUty distribution function in the 
component space associated through to the phase space equihbrium probabihty 
distribution function V^'i{C). Using the detailed balance condition (|38|l we obtain, 

>v"q(7^^|7^) _ v'lin') 

which is the detailed balance condition in the component space. 

In conclusion, in the component space the rates W satisfy the same set of 
conditions as the rates W except for the detailed balance condition which in the 
component space only holds at equilibrium. In general, rates W do not satisfy 
detailed balance so the FDT derived in the previous sections for the microscopic 
master equation will not be valid in the component space. They are, however, 
valid in equilibrium so, for example, after an appropriate redefinition of correlations 
and responses the equilibrium fluctuation-dissipation theorem (|60ll still holds in the 
component space. 

4. FDT extensions to the non-equilibrium regime 

In this section we discuss how to extend the previous ideas to the non-equilibrium 
glassy regime. After a brief introduction of aging (intended for the non-specialist) we 
discuss how to derive a free-energy master equation by introducing the configurational 
entropy and the notion of the effective temperature. We then discuss how to extend 
the FDT beyond equilibrium by introducing the fluctuation-dissipation ratio and 
the quasi-FDT. This requires the notion of neutral observables. Finally, a possible 
partitioning scheme is presented. 

An intermezzo on aging 

In this section we discuss one of the main signatures of non-equilibrium regime of 
glassy systems, i.e. the existence of aging and the quantifying of FDT violations 
through a modified FDT in terms of a set of effective temperatures. This discussion 
is intended for the inexpert reader who wants to have a glance on the key aspects of 
glassy systems before entering into the more detailed exposition. Therefore, the level 
of this discussion is highly introductory and the expert reader can move directly into 
the next section. 

In relaxational glassy systems the two basic properties of correlations and 
response, time-translational invariance (TTI) H42|I and the FDT (|6Uf) . do not hold 
anymore. In particular, TTI is observed to be violated in the following way: both 
correlations and responses decay slower as the system gets older, i.e. the system 
is aging. This fact stems from several experimental observations in polymers and 
deformable materials PP, structural glasses PHI and spin glasses ^3). In driven systems 
the situation appears less complicated as TTI holds and only the FDT is violated. 
Some of the aspects described below carry over also to driven systems, however for 
sake of clarity we will stick in what follows to the case of aging systems. 

Experimentally aging is manifest through the measurement of the so called 
integrated response function (IRF) or time-dependent susceptibility described in 
Section 14.31 In aging systems the weak long term memory makes response functions 
R{t, s) hardly measurable as they asymptotically decay to zero, instead it is easier 
to measure the cumulative response or time-dependent susceptibility x(^i = 
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R{t, t')dt' (defined after perturbing the system at the waiting time s -often denoted 
as t^-) which in general are finite: in dielectric measurements of glasses the IRF 
corresponds to the polarizability of the sample after cutting or applying an electric 
field; in magnetic systems it corresponds to the thermoremanent or zero-field cooled 
susceptibility; in mechanical systems aging is observed by measuring the deformation 
of the sample after applying a tensile load. All these measurements reveal that s) 
is well approximated by the sum of two contributions, 

X(t,5) =Xst(t-s) + Xag(i,s) (70) 

where the first contribution Xst(^^s) stands for an stationary part that asymptotically 
decays to a finite value or plateau and Xag(^;'S) is the aging part that is well 
approximated by the scaling relation, 

Xag(i,s)=x(^) (71) 

This scaling behavior, that is obtained within many solvable models, is known as full 
aging or simple t/ s scaling %. However, deviations from full aging have been reported in 
many experiments suggesting that this simple scaling behavior does not fully account 
for the experimental data. In particular, spin glass measurements clearly favour a 
subaging scenario where x has as scaling argument the variable p- with 5 < 1 (we 
note, however, that experiments report values around 0.95, i.e. very close to 1). The 
physical origin of these deviations is still unknown. For both correlation and response 
functions similar decompositions as in (|70|l , H71|l are expected to be valid but replacing 
X by C or i? (however, for the aging part of the response there is an additional factor 
1/t multiplying i?ag(t, s)) . Correlations and responses are difficult to experimentally 
access. In theoretical or numerical simulation studies, the calculation of correlation 
functions is always preferred. Typical curves for the susceptibilities or correlations are 
depicted in Figure |2| Within the full aging scenario (|71|l it is usually shown that, in 
the asymptotic limit where both s and t are large, the FDT (|60|1 is violated depending 
on the ratio it - s)/s. If {t - .s)/s < 1, C(i, s) ~ Cst(t - s) and FDT holds, 

^^^^ = mt(t~.). (72) 

OS 

However, in the other case {t — s)/s> 1, then C(i, s) ~ Cs,g{t, s) and the FDT is 
violated according to the new relation, 

^^|M=T,ff(.)i?,g(t,.) (73) 

where Teg{s) is a new parameter that enters into the new relation playing the role of an 
effective temperature. The particular way in which the FDT is violated suggests calling 
the new relation 1)73(1 a quasi fiuctuation-dissipation theorem (QFDT). Although other 
possible terms have been used to refer to the modified FDT, here we will adhere to 
the term QFDT, as this is the one originally introduced by Horner 32 that expresses 
the idea of partial equilibration among a subset of degrees of freedom. 

Aging carries associated the decomposition of time into time sectors. In the 
previous example of full aging 1(71(1 there are two time sectors depending on the ratio 
{t — s)/s as described in 1(72(1 . 1(73(1 . However, as we already said, deviations from the 
full aging behavior are expected to be present in general. In those cases, the stationary 

f Coarsening system also deviate from the simple form 1711 in favour of XagityS) = s^^xi^/^) with 
a > 0. 
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Figure 2. Typical sliape of a two-time correlation function C(t, tm) (the same 
would be valid for the integrated susceptibility x(*i*™)): plotted as a function 
of t — tm (in log-scale) with system age increasing from left to right. The 
decomposition 1701 shows the two separate time sectors; in the first fast part of 
the relaxation C{t, i^j) — Cst(t — t-w) is independent of s and obeys TTI while 
the slow second part Cag {t — tio) decays from the plateau in a timescale growing 
with tw ■ This Figure has been taken from a lattice gas coarsening model [Ml with 
properly normalized correlations at equal times C{t,t) = 1. 



result (|7^ for the short-time sector still holds but (f75)l is replaced by the more general 
relation, 



ds 



Tcff(ag(t,s))i?ag(i,s) 



(74) 



where the new effective temperature Tcff (Cag(i, s)) = Tcs{C{t, s)), since C = Cag 
when {t — s)/s ~ 0(1), and depends on both times only through the value of C. In 
this case, the aging part develops time sectors defined as those values of t, s where 
t/s ~ 0(1) U3, 34 . Each sector is then labeled by the value of the correlation function 
C {t, s) and many effective temperatures arise in the description of the non-equilibrium 
regime, the QFDT (|74|l quantifying FDT violations within each sector. Glassy systems 
are often classified into three different groups according to the dependence of Tcs{C). 
For coarsening systems Tes{C) only takes two values: T and infinite for the stationary 
and aging regimes respectively. §. For structural glasses Teff(C) also takes two values: 
T and Tcsis) > T for the stationary and aging regimes respectively. These are often 
referred as two-timescales systems. Finally for spin-glass systems rcfr(C') takes a 



§ We note however that in some special cases coarsening systems can display more complex behaviour 
of Tg{i{C), see Sections 16. 617. 3l and the discussion in Section 13.1 in |3j. 
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continuous spectrum of values extending from a lower bound T* > T up to infinity. 
These are known as many-timescales systems. All these three limit cases correspond to 
a well known static low-temperature description in the framework of spin-glass theory 
| ci5| in terms of replica-symmetry breaking (RSB) : coarsening systems are those where 
replica symmetry (RS) is unbroken, structural glass systems correspond to one-step 
of RSB while spin-glass systems corresponds to full RSB. As a particular example of 
one-step systems there are some models (such as entropy barrier models, see section 
lf).5|l that display glassy behavior only at zero temperature. The stationary regime is 
then absent in these models and their non-equilibrium dynamics is characterized by 
a single effective temperature, the bath temperature being zero. Along this review, 
we will often refer to them as one-timescale models. The three possible scenarios are 
depicted in figure 01 (see section IT^ for a more detailed exposition). The experimental 
challenge of these ideas remains one of the most awaited results. 

4-2. The unbiased component ensemble and the master free energy equation 

One of the key ideas behind the existence of a quasi fluctuation-dissipation theorem 
(QFDT) in aging systems is the emergence of a non-equilibrium ensemble in the 
asymptotic long-time regime of the relaxation process. A related non-equilibrium 
ensemble could also emerge in driven stationary systems. Although the nature of this 
ensemble is yet to be understood we can anticipate some of its main properties. Some 
of these ideas have been already presented Here we present more elaborated 

work. 

In a description of glassy phenomena where the system is kept in contact with a 
thermal bath at temperature T the energy is not constant and cquipartition does not 
necessarily hold. Therefore this ensemble is neither microcanonical or canonical but 
of a more complicated nature. 

The possible existence of a non-equilibrium ensemble traces back to Palmer |37| 
who introduced the unbiased component ensemble to characterize the equilibrium 
sampling of phase space components. Let us consider a given partition of phase space 
into components (see Section IX^ and let us define the free energy F{Tl) of a given 
component TZ by. 



It is possible to extend this idea to the non-equilibrium regime by assuming an equal 
probability hypothesis: different components with identical free energy F(TZ) do have 
the same probability. 



In what follows, we will use the letter T (as opposed to F{TZ)) to denote component 
free energies defined in (|75|l after dropping off the explicit argument TZ. The existence 
of the unbiased ensemble tantamount to the appearance of a new measure based 
on free energy rather than on energy. Contrarily to Palmer 37, (who assumes 
an equilibrium probability distribution for P(J^, i), see H87|l below) the probability 
distribution P{J^,t) is unknown, time dependent and must be found as a solution 
of a master equation (ME) as follows. To derive the free-energy ME we define the 
probability density, 




(75) 



ceTZ 



v{'Ji,t) = p{F{n),t). 



(76) 




(77) 
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where we have used (|76|) and the definition, 

n{J=-,T) = Y,S{^ - F{T^))- (78) 

and we have introduced exphcitly the temperature dependence in fl to stress the 
temperature dependence of the free energy H75|l. Although consistency requires to 
add the T dependence also to J- here we drop this dependence in order to lighten the 
notation. Equation (|77|) describes the probability for the system to be in a component 
of free energy at time t. We have indicated it in bold to distinguish it from the 
probability V{TZ,t). Summing (|64|l over components having identical free energy T 
we get, 

= P(r , t)Z,{T\T') - ^ V{T, t)Z,{T'\T) (79) 
with the conditioned probabilities Zt defined by. 



Z,{T\T') = — — - Y W(7^|7^'; t)S{J^{'JZ} - T)5{T'{n') - T')P{n' , t) (80) 

where the W have been already defined in 1)6711 . Note that both W and Zt are 
time-dependent rates. Again, as for the transition probabilities (I65II . the new rates 
Zt{J-\J-') do not satisfy detailed balance but satisfy the other requirements (non- 
negativeness, ergodicity and causality). Expression H8U|) is exact but intractable. As 
we are postulating the existence of the unbiased component ensemble (|76ll , consistency 
in the component master equation H64() implies that the rate W{TZ,Tl' ;t) is a time- 
dependent function of the initial and final components, only through the value of their 
free-energies TiJZ'),T{TZ), 

win,n';t) = Wt{T{n),T{n')) . (si) 

The transition rates (|8U|) can be further simplified, 

Zt{T\T') =WtiT\T')n{T,T) (82) 

The quantity f2(J^, T) is exponentially large with the volume of the system and defines 
what we will denote as the configurational entropy or complexity 5c(.F, T) f, 

n{J^,T)^exp{Sc{J',T)) (83) 

Therefore, all the information on the master equation (|79|l goes into the density of 
components Sc{F, T) and the rates 'Wt{J-\J-')- These contain all the information about 
the properties of the unbiased ensemble. 

The description of glassy dynamics in terms of a master free-energy equation 
such as H79|l has been wandering around for many years in the literature of the field. 
Several equations have appeared scattered in the literature during the last decades, 
but generally written in terms of the energy instead of the free energy, see for instance 
These equations describe what are usually known as trap models (see Section 
I6.4|l . Examples master equations proposed by Dyre [3^] and Bouchaud 00]. Other 
attempts include granular media j41|. 

t The term configurational entropy has often been used with difli'erent meanings, leading to confusion. 
Originally, as used by Adam and Gibbs for their thermodynamic theory, it denotes the part of the 
total entropy including only the configurational degrees of freedom. More recently, in the context 
of spin-glass theory, this concept has been coined to denote that part of the configurational entropy 
that counts the number of metastable states rather than configurations. It is with this last meaning 
that we understand it here. For a throughout discussion of this concept see section f5.il 
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4-2.1. Complexity and the effective temperature Before closing the present discussion 
let us note that, in equilibrium, both the transition rates Zt{J-\J-') and Wt{^\^') are 
time-independent and satisfy detailed balance, 



= exp(-/3(<f>(^',r)-$(^,T))) (84) 
= exp(-/3(J^' - T)) (85) 



with 

^{T,T) ^ T ~TSc{T,T) (86) 

where $(JF, T) is a new potential where the free energy appears balanced by the 
complexity. The probability (|77|l assumes the simple form, 

P.(^, = ,87) 

where Z{P) is given by p4|) . We will describe in Section [5.1.21 how this relation 
provides us with a tool to obtain the configurational entropy as function of both the 
free energy and the temperature. 

The hints at the existence of this new potential were found by Kirkpatrick, 
Thirumalai and Wolyness |42[ I43L IH| who identified the marginal transition Ta 
("A" standing for activated) as the temperature below which different metastable 
states concur in such a number to compensate the lower equilibrium free energy 
of the paramagnetic or liquid state. The subindex in Ta stands for the fact that 
below that temperature activation is dominant and relaxation occurs in the form 
of activated jumps from one metastable state to another. Ta is identified with the 
mode-coupling transition temperature Tc of mode-coupling theories (see section [6. Ill 
and corresponds to a spinodal instability 02 I43L E| . The mathematical argument 
behind the compensation of the free-energy of metastable states by the complexity is 
as follows. Let us decompose the canonical partition function of the system as a sum 
over a set of non-overlapping components TZ (as explained in section 3(1 . 

Z{T) ^Y.'^M-mC)) = ^exp(-/3.F(7^)) 
c n 

= J2 ^i^^ T) exp(-/3.F) = T)) (88) 

where we have used the definition H75|) for the free energy of components J'^JV) and 
H83|l . (18611 . Due to the extensive character of the variables T^Sc{T,T^ and $(J^, T), 
the dominant contribution to the sum in H88() is evaluated through the saddle point 
method. At each temperature T there is a free energy F*{T) such that its contribution 
to the exponent 1(88(1 is dominant, i.e. Z{T) ~ exp(— /9(f>* (T)) where we have defined 
$*(T) = <i>(F*(r), T). The behavior of this solution depends on the shape of the 
function S'c(.F, T). In general this function is a monotonically increasing function 
of the free energy J-. Because the exponential is a positive definite function we 
have $(F,T) > Fpara(r) where Z{T) = exp(-/3Fpara(T)) and Fp,,,{T) denotes the 
paramagnetic free energy. 

Above Ta there is no solution F*{T) which can compensate the equilibrium 
paramagnetic free-energy and ^{F, T) > Fpara(r). Below Ta a solution appears F*{T) 
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such that $*(T) = Fpara(r) and gives the dominant contribution to so F*{T) 
satisfies the saddle-point relation, 

1 _ dSc{T,T) 



T dT 



(89) 

The identity $*(r) = Fpara(T) implies Fpara(T) = F*{T) - TSc{F* {T),T). This 
means that for T < Tc there is a band of components with free energy F*{T) > 
FpaiaiT) (therefore with free energy above the equilibrium one) whose difference with 
Fpe^raiT) IS Compensated by the complexity S*(T) = S'c(F*(T), T). This solution ex- 
ists as long as S*{T) > 0. Because Sc{F,T) is a monotonically increasing function 
of F and both F*{T), S*{T) decrease with T there is a temperature Tk at which 
S*{Tk) = 0. Below this temperature, the complexity vanishes and the solution F*{T) 
ceases to change with temperature (so equation does not hold anymore) but sticks 
to its minimum value F*(Tk). This is the entropy crisis scenario where Tk corresponds 
to the Kauzmann temperature j45| . 

In mean-field models it has been shown ^] that the complexity Sc{F, T) defines 
a free-energy dependent effective temperature through the relation, 

reff(^) ~ dJ^ ^^"^ 

From this relation, it emerges that the configurational entropy plays the role of 
the thermodynamic potential associated to the effective temperature. Similarly, 
the entropy is the potential conjugated to bath temperature in the microcanonical 
ensemble. However, (|9U|I has been only derived in mean-field models close to 
the asymptotic free-energy threshold where lower free-energy states are inaccessible 
| 4t)l 14 7| . In models that do not have a marginal free-energy threshold above the 
equilibrium value, it is possible to show that free-energies are uncorrelated random 
variables exponentially distributed 0S1 ISH]' The extension of these results beyond 
mean-field where all free energy states are accessible through activated processes is 
at the roots of the existence of the unbiased component ensemble. Some attempts 
have been proposed in 0^1 HHI In particular it would be very interesting to 

understand the general form of the transition probabilities Zt as these lead to very 
specific predictions amenable of numerical checks. The exponential character of the 
free-energy distribution of the lower free-energy states below the threshold, and the 
fact that this distribution is time dependent (as shown by the fact that Zt is itself 
time dependent) appear as the two crucial ingredients to understand the emergence of 
effective temperatures in glassy systems. These features are present in many models of 
glasses such as, mean-field spin-glass models, trap models (see Section or entropy 
barrier models (see Section FS. 5. 



^.3. The integrated response function (IRF) and fluctuation- dissipation (FD) plots 

When a system is in a non-equilibrium state its response to a external perturbation 
cannot be described, in general, by FDT relations such as (|60|) since this has been 
derived assuming that the system is in a stationary state. However, the glassy state is a 
particular non-equilibrium state characterized by extremely slow relaxation processes. 
Hence, while the system is not in a real equilibrium state, it may be thought as 
being in a sort of quasi-equilibrium regime over time-scales much longer than the 
microscopic time-scales but still smaller than the typical relaxation time-scales of 
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the slow processes. In this quasi-equihbrium regime the evolution of the system is 
quasi-stationary since non-stationary effects are seen only for times of the order of 
the time-scales of the slow processes. In this situation we may think that relations 
similar to (|60|l can still be valid, even if TTI cannot be assumed anymore. Thus a 
possible generalization of FDT to the glassy regime requires to introduce the following 
nondimensional quantity: 

TRA^B{t,s) 

XA^sit, s) = -g— — — , t > s (91) 

■§^CA.B{t,s) 

where XA.B(t,s) is called the Fluctuation-Dissipation Ratio (FDR). In equilibrium 
XA,B{t, s) = 1 whatever times t, s and observables A, B are used [Cfr. (|60|l ]. Thus X 
is a measure of the violation of the true equilibrium in the quasi-equilibrium glassy 
state. The validity of H91|l . i.e. the proportionality between the response and the 
time derivative of correlation function, can only be checked a posteriori since it is 
based on the quasi-equilibrium hypothesis that up to now it has been proved only 
in mean-field models. Eliminating one time in favour of the correlation function, the 
time dependence of XA,B{t, s) can be recast in the form 

XA,B{t. s) EE XaMCaM*^ s): s]- (92) 
The FDR was first studied in spin glasses where analytical results have shown that 
glassy systems in general satisfy the weak ergodicity breaking scenario _40; , discussed 
in Section 16.21 For the present purpose it is enough to note that calculations in 
mean-field spin-glass models |52 H2] have shown that in the limit s — oo the FDR 
is a non-trivial function which depends on the relation between the times t and s 
only through the correlation function CA,B{t, s). The following specific form of FDT 
violations has been proposed to be generically valid in the non-equilibrium regime of 
glassy systems, 

lim XA,B{i^s)= XaACaA^^ s)] (93) 

S — fCfO 

Using (inSJ and we obtain the differential form of the Quasi-FDT (QFDT) relation 
PI, 



XaAC) 



TRA.B{t,s) 



CA,B{t, s) 



t > s. (94) 

CA.Bit,s) = C 

which describes the response of the system in the (quasi)-equilibrium state to a 
impulsive perturbation at time s < t. Experiments and numerical simulations 
usually measure integrated response functionsf (IRF), i.e., the response at time t 
to a perturbation switched on or off at time s < t. According to the definition (|43|l 
the variation of the observable {A{t))f^ to linear order in the perturbation intensity e 
is given by, 

(A(t)),^ = (^(t))o + ei?A,B(i,s)+0(e'), t > s. (95) 
Assuming that the perturbation acts for all times t > s and that its intensity is small 
enough for the accumulated response to be linear in e we get 

t 

{A{t)),^ = {A{t))o + Y,^t'RAMt,t') 

t'=s 

^{A{t))o + J^ dt'e{t')RA,B{t,t') (96) 

j The integrated response functions are also called time dependent or non-equilibrium susceptibilities. 
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where Ct' (or e{t')) is the perturbation at time t'. In the particular case of e(i') = 
e9(t' — s), i.e., of a constant perturbation, one gets 

XAMt^s) = ^^^imU^Smi^ f dt' RAMt,t'). (97) 

which is also called zero-field cooled (ZFC) susceptibility. The name zero-field 
cooled follows from the experimental protocol used in spin glass measurements to 
distinguish it from the thermoremanent magnetization (TRM). To measure the TRM 
susceptibility a constant external perturbation is applied to the system at time t — Q 
and removed at time s > and the subsequent decay of {A{t)) is recorded. The TRM 
susceptibility is given by 

xlTit^s) = (^lWk^M!))£ = r dt'RA.BM (98) 

In the large s limit the ZFC H97|l and TRM H98|l susceptibilities are equivalent since 
they are related by, 

XA,B (i , s) -f xYb (i , s) = (i , 0) (99) 

and XA,B{t ^ oo, 0) = x^'Zb- 

Inserting the QFDT relation in we obtain the formula 



1 



CA,B{t,t) 



XA.B{t,s) = - dC'XAAC) (100) 

^ JCA,B(t.s) 

which relates XA,B{t,s) to the FDR Xa,b and provides a simple way to calculate 
Xa.b from measurements of CA,Bit, s) and XA,B{t, s) in the time sector t > s. Suppose 
indeed we fix the lowest time s and plot XA,Bit, s) as function of CA,B{t, s) for different 
values of t, then the value of the FDR can be obtained from the slope of the resulting 
curve. In many of the examples considered in this review the equal times correlation 
function is time- independent, for instance CA,Bit,t) = 1. In this case the slope can 
be simply obtained by derivation of XA.Bit^s) with respect to Ca,_b(^, s) for fixed s 
which from l|l()()(l yields, 

dxA,B{t,s) 



Xa.b{C) = -(3 



(101) 

Cyi,s(f,i)— const, s fixed 



dCA,B{t, s) 

Typical FD plots are shown in FigureOlfor the three possible scenarios (see section^^. 
In the general case in which C^i ^(i) t) is time-dependent one needs to be more careful 
in computing the FDR. SoUich and coworkers have proposed |^ to construct FD 
plots X versus C with t kept constant and varying the lowest time s. From (|91|l we 
have, 

dxA,B{t,s) 



XaMC) = -/3 



(102) 

t fixed 



dCA,B{t, s) 

If C(t, t) changes with time it is convenient to normalize correlations and the IRF by 
the equal times correlation CA.B{t,t): 

^ ^ Ca,b{^,s) ^ XA,B{t,s) 

CA,B[t,s)^— — — , XA,B(t,s)^- — — . (103) 

CA,B[t,t) CA,B[t,t) 



With these definitions, 

dxA,B{t,s) 



XA^BiC) = -(3 



dCA,B{t, s) 



(104) 

t fixed 
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Figure 3. FD plots for the three possible scenarios. From bottom to top: 
1) RS models or coarsening systems with two-timescales and two-temperatures 
(one identical to the bath, the other infinite), 2) one-step RSB models with 
two-timescales and two-temperatures (one identical to the bath, the other finite 
but higher than that of the bath), 3) full RSB models with many timescales 
and temperatures. For C > Qea (the Edwards- Anderson parameter), i.e. the 
stationary regime, all models satisfy FDT (the dashed line). 



The importance of normalizing the raw FD plots (|l()2|l is well appreciated in trap 
models discussed in section l^^ or in kinetically constrained models discussed in section 
17.51 In this last case, for example, raw FD plots can lead to awkward representations 
as the one shown in Figure OBI 

4- 4- The concept of neutral observables 

If the FDR H91|l . (19411 has the physical interpretation of a temperature (as has been 
suggested, see the discussion in Section IF?^ then one would expect the FDR to be 
independent on the observables A, B used to construct correlations and responses. In 
fact, this is true in equilibrium where Xa,b = 1 whatever A,B. However, although 
Xa,b 7^ 1 observable independence is not at all required in the glassy regime. In this 
section we present a brief digression on which conditions the observables A, B must 
satisfy for the FDR to be observable independent. This issue is yet unresolved, so the 
present discussion is quite speculative. 

Albeit restricted, for simplicity we will consider here the case of a glassy system 
with only two-timescales where A — B (so we will denote Xa,a simply by Xa) 
in the time sector where (t — s)/s ^ or Xa ^ 1- In equilibrium one could 

argue that the equality Xa = 1 is related to the fact that the entropy S{E,A), as 
defined in the microcanonical ensemble Section f2.1l is maximum (|12|l for A equal to its 
equilibrium value. This property, is observable independent as well as it is the identity 
Xa = 1. a similar argument, but extended to the glassy regime, would require to 
define the configurational entropy Sc{^,A) (i.e. the equivalent generalization ^ of 
(IHSl) where J^,A denote the component averaged values of the free-energy l(7^ and 
the corresponding restricted Gibbs average for the observable A. We could then say 
that A is a neutral observable if its dynamically averaged value at all times < A{t) > 
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t 

Figure 4. Illustration of the neutrality property of observables. The value of a 
neutral observable remains fixed to its initial (equilibrium) value despite of the 
fact that the initial state is non-equilibrium. 

coincides with the stationary maximum of the function Sc{J',A) as A is varied. Of 
course, if this were not true the value of Xa would then depend on the value of A 
in the same way that the value of the temperature l/T — P in the microcanonical 
ensemble, and for a given value of the energy E, would depend on the value of the 
observable A if /i in (|83|l were not zero. 

For instance, the magnetization in mean-field spin-glass models is known to 
be a neutral observable and nearly all computations of the FDR have used this 
observable (see Sections 16. 2() . In fact, in the framework of the TAP approach it can 
be shown that the configurational entropy evaluated as a function of the free-energy 
and magnetization of the TAP states is maximum at zero magnetization. Indeed, 
the fact that the magnetization is not a good order parameter in these models (it 
vanishes in both the paramagnetic and the spin-glass phase) is related to its neutral 
character. Not by chance the majority of numerical studies in glassy systems use the 
magnetization as the central observable to investigate and measure FDT violations. 

A non-neutral observable A would correspond to a situation where, starting at 
time zero from a non-equilibrium initial state where the value of A is taken to coincide 
with its equilibrium value (this is not a contradiction, since A could coincide with its 
equilibrium value but not the value of any other observable) the subsequent evolution 
of A deviates from its initial equilibrium value. In this case A is not neutral because 
its time evolution is correlated to that of the other observables. On the contrary, if 
A were a neutral observable, then it would stick to the value A forever (as happens 
for the magnetization in mean-field spin-glass models). Fig 0] illustrates this behavior. 
The definition of neutrality can be easily adapted in trap models by assuming that 
observable values of traps are uncorrelated with their energies (see also Section 

4-5. Numerical approach to component dynamics: Stillinger-Weber decomposition 

The increase of computational power and the recent developments in the theory of 
disordered systems has pushed forward an approach to the glass transition based on 
the analysis of a reduced dynamics in the component space. The underlying ideas date 
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back to more than 30 years ago to a seminal, but talkative, paper of Goldstein |56| . 
The glass transition is of purely dynamical origin and hence must reflect the properties 
of the dynamical evolution of the system. Goldstein suggested that the dynamics of a 
supercooled liquid can be understood in terms of a diffusive process between different 
basins of the potential energy surface. At low temperature the dynamics slows down 
since the system gets trapped for a long time in a basin. This approach, which focus on 
the topological properties of the energy surface is rather appealing since it naturally 
leads to a convenient framework for understanding the complex phenomenology of 
glassy systems. 

The implementation of these (qualitative) ideas came some years later by 
StilHnger and Weber [571 1551 ISU) (SW) who formahzed the concept of basin in 
configuration space identifying it with a component in the component space, and 
proposed a procedure for identifying them: the set of all configurations connected 
to the same local energy minimum by a steepest-descent path on the energy surface 
uniquely defines the basins of the minimum. Stillinger and Weber (SW) called the local 
minimum inherent structure (IS) to stress its intrinsic nature. Since the identification 
of ISs is unique, the mapping from configurations to local minima gives a unique 
well defined decomposition of the phase space into a disjoint set of basins. The SW 
decomposition defines a mapping of the phase space to the component space in which 
each basin, usually labeled by the energy Eis of the local minimum, is a component. 

This decomposition does not cover completely the configuration space since it 
leaves out the boundaries between different basins. However, under the assumption 
that those configurations do non contribute to the thermodynamic of the system, e.g., 
the boundaries between basins are sub-extensive, it does cover almost all the phase 
space and the partition function can be written as a sum of contributions from different 
components: 

^(T)~^ Zis(Sis,r) (105) 

Bis 

Let il{E) denote the number of IS with energy Eis — E, then collecting all components 
with the same value of Eis 

ziT)^yn{E) V Z,siT)/n{E) 



E Eis=E 



e^p[SciE)-pFi,iT,E)]. (106) 



The term 

Sc{E) = \ogn{E), (107) 

which accounts from the entropic contribution arising from the number of different 
basins with the same IS energy, is called the SW configurational entropy of complexity. 
This quantity is strongly related to the partitioning, so we add the adjective SW to 
distinguish it to other definitions of configurational entropy taken from mean-field 
concepts. 

The second term Fb{T,E) is defined as 



FbiT,E)^-T log 



^ ' Eis=E 



(108) 



In general this quantity differs from the average free energy of components with Eis = 
E, however if all these components have similar statistical properties, then Fh{T,E) 
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is the free energy of the system when constrained to any one of the components with 
Eis = E. In the thermodynamic hmit the system populates components with energy 
Eis — Eis{T) fixed by the condition 

- PF{E) = Sc{E) ~ f3 Fb{T, E) = maximum over E (109) 
and the free energy of the system can be calculated using 

F{T) = Ft,[T,Eis{T)] ~TS,[Eis{T)] (110) 
The condition (|109|l is equivalent to that of F{T) being minimal, i.e., 

dF _ dF,{T,E) dS^jE) _ 

dE dE dE ' ^ ' 

Note that the minimum condition follows from the balance between the contribution 
from the change with the energy of the shape of the basins {dFb (T, E) / dE) and its 
corresponding number {dSc [E] / dE) . Often the free energy is written as 

i^b(T, Ei^) = Eis + F,{T, Eis) (112) 

The first term in (|112|l takes into account the average energy of IS visited in 
equihbrium at temperature T, as can be seen from (I110|l : U{T) = d{f3Fi,)/df] = 
EisiT) + d{PFy)/df3. It can be shown (HOI ED that if the density of states n{E) is 
Gaussian and the basins have approximately the same shape then Eis oc 1/T. The 
second term in H112|l describes the volume of the corresponding components and is 
called the "vibrational" contribution. 

To understand the success and limitation of the IS approach we have to analyze 
the idea behind the SW approach. It is clear that even if the phase space can be always 
partitioned, not all possible partitions will lead to a physically relevant dynamics in the 
component space. This is a well known problem in the theory of dynamical systems, 
where the component dynamics is called symbolic dynamics, see e.g. Ref. |22|. To 
prove that the SW is a physically good partition for a given system is a problem of 
the same hardness as proving ergodicity. One then adopts a constructive point of 
view, along the same lines of equilibrium statistical theory: based on some reasonable 
hypothesis one first assumes that the SW partition is a good partition and then check 
if this reproduces the desired features of the dynamics. 

The physical motivation behind the SW proposal follows from the observation 
that the potential energy surface of a super-cooled liquid contains a large number of 
local minima and that the time evolution can be separated into two different processes: 
thermal relaxation into basins [intra-hasin motion) and thermally activated potential 
energy barrier crossing between different basins (inter-basin motion). This scenario 
has been recently confirmed from numerical analysis [601 1611 1621 163L ISl] . The time- 
scales separation of the two processes strongly depends on temperature. When the 
temperature is lowered down to the order of the critical mode-coupling theory (MOT) 
temperature Tc the typical barrier height is of the order of the thermal energy ksTc, 
and the slow inter-basin motion dominates the relaxation dynamics. If the temperature 
is further reduced the relaxation time eventually becomes of the same order of the 
physical observation time and the system falls to a non-equilibrium state since there 
is not enough time to cross barriers and equilibrate. With this picture in mind it is 
natural to view the IS partitioning as the natural elements to describe the slow glassy 
dynamics. This approach is rather appealing since naturally leads to universality: all 
glassy systems with similar IS dynamics must have similar glassy behaviour. Recent 
IS analysis performed on disordered spin systems displaying a transition of fragile 
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glass type do support this conclusion [SSI IMl IHH EHl HHI ■ It should be noted [021 that 
the definition of IS for spin systems is more subtle than for systems with continuous 
variables. Indeed usually for spin systems IS are defined as one-spin fiip stable states, 
however these may be not stable for two-spin or higher number of spin flips. One 
possibility of making IS well defined also for spin systems is to define them directly 
from the T = limit of the dynamics, i.e., as states which are stable under the T = 
dynamics |7()j . This is the definition used in this review when discussing IS for spin 
systems. 

5. Thermodynamic description of the aging state 

We saw in section [4.21 how the self-generated dynamical measure allows a description 
of the aging dynamics in terms of a probabilistic master equation with transition 
rates characterized by an extensive quantity that was defined as a configurational 
entropy or complexity (|83|l . This quantity has received considerable attention in 
studies of spin-glasses since the seminal paper of Thouless, Anderson and Palmer 
(TAP) [^ on the SK model where a way to compute the configurational entropy 
was proposed [221 ■ Later studies in the context of structural glasses 021 QHl 121 have 
shown its importance as the mechanism for an entropy crisis of the supercooled liquid 
as proposed by Kauzmann many years ago j45| . 

5.1. Methods to compute the complexity 

In this section we present a schematic overview of some of the analytical and numerical 
methods that have appeared in the literature to compute the configurational entropy. 
In the absence of a full solution of the dynamics in many systems, and under the 
assumption that there is a connection between the effective temperature and the 
configurational entropy (see the discussion in Section I4.2.1() , the calculation of the 
later, by using equilibrium methods taken from statistical physics, appears as an 
alternative way of quantifying FDT violations. In mean-field theories metastable 
states give a natural partition of the phase space since their life-time diverges in the 
thermodynamic limit. For systems with short-range interactions, however, metastable 
states can be defined unambiguously only referring to some reference time-scale. 
Therefore the identification of metastable states for real systems can be a very hard 
task. In section [4. 51 we have presented a partition scheme, proposed by Stillinger and 
Weber, which in principle can be applied to any system. The scheme essentially uses 
a zero-temperature dynamics and thus it is free from the ambiguities due to the finite 
metastable life-time. The results described in next sections must be seen as instructive 
attempts to evaluate a quantity (the complexity) that governs the slow dynamics of 
relaxational glassy systems. The extension of these equilibrium concepts to other 
non-equilibrium systems beyond aging systems (e.g. driven systems) remains an open 
problem. 

5.1.1. Analytical methods. Bray and Moore [721 calculated Sc{^,T) for the 
Sherrington-Kirkpatrick model within the TAP approach. The TAP equations give 
the local magnetization nii in a system confined to a metastable state, which for 
mean-field models have infinite life-time. As a consequence the number of metastable 
states (i.e., components) can be readily obtained just counting the number J\fs{J-', T) 
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of solutions of the TAP equations at temperature T with a free energy 

N 

Afs{T,T)^ I W dmMeiH{{m,})\ 



,(^,T)== / Y[dm,\detH{{r 



x5{j^-J^{{m,},T))S{g,{{m,})) (113) 

where J-{{mi},T) if the TAP free energy at temperature T as function of the local 
magnetizations TOj, gi{{mi}) = d!F{{mi\,T)/dmi = arc the TAP equations and 
Hij — dgi{{mi]) / druj the Hessian. This type of calculation has been done for other 
mean-field models, such as p-spin [72| and Random Orthogonal Model (ROM) |74[I75| '). 
finding in all cases that J\fs {J^, T) increases exponentially fast with the system size 
system N. This remains true if the number of free-energy minima Afmi^,T), instead 
of the number of stationary points Ms T) , is considered |7B] . Although these type 
of calculations can be done only in exactly solvable mean-field models the exponential 
growth with the system size of the number of free energy local minima or stationary 
points is generally applicable to any system (mean-field or not) displaying glassy 
behaviour. The knowledge of the number of minima allows to define the complexity 
lf7S)l as, 

Sc{T,T)^l0gMr,^{T,T) (114) 

and hence the thermodynamic potential $(JF, T) = T — TSc{J-,T) as described in 
section ^21 

A general framework to evaluate the complexity has been devised by Monasson 
|77). The starting point in his procedure is to consider m interacting copies or 
replicas of the original system, with an attractive interaction term of the form 
e^™^^-^ (3(Co,Cb) where Q{Ca,Cb) is a suitable overlap function which takes its 
maximum value only if Ca — Cb- The free energy of the replicated system is then, 

m m 
g-/3F(T,™)^ ^ CXp[-/3^W(Ca)+/3e ^ Q(Ca,Cb)] (115) 
Ci,...fim a=l a,b=l 

If the thermodynamic limit is taken before the limit e — )■ 0^ then the configurations 
Ca, Cb tend to lie as close as possible since maximization of the coupling term minimize 
the global free energy F^™'\T) and hence, given a phase space partition, the replicas 
tend to "condensate" into the same component. Thus collecting all components with 
the same free energy the partition function Z{T, m) can be decomposed as 

Z{T, m) = ^ T) eyip{~m(3T) 

= ^exp[-/3$(.F,r,m)] (116) 

where ^{J-,T,m) = mJF — TSc{^,T) is basically the potential $(JF, T) discussed in 
section l^^ with the term multiplied by m (the order of limits, first volume — > oo and 
then e — !■ 0+, enforces the m replicas to occupy the same component TZ). In the limit 
m ^ 1 we recover the potential $(J^, T): ^{J^,T,m = 1) = $(JF, T). The knowledge 
of ^{J^ ,T,m) allows to compute the configurational entropy. In the thermodynamic 
limit the sum in H116() is dominated by the free energy F* (T, m) that satisfies the 
relation. 



T dT 



(117) 

T=F-'(T,m) 
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Inserting the solution F*{T,m.) into ^{!F,T,m) we obtain the free-energy potential 
$*(T, to) = (^{F* (T,m),T^m). If m (originally an integer value) is continued to real 
values then it can be show that the following relations are satsified, 



Varying m allows to compute the configurational entropy Sc{P , T) as function of the 
two variables T and T. The potential $(^, T, to) can be explicitly evaluated with the 
sole knowledge of the microscopic Hamiltonian of the system and using the replica 
method. Although this procedure was initially applied only to mean-field disordered 
systems [77| |7S], more recently it has been extended to more realistic interacting 
potentials such as Lennard- Jones liquids '75', "SD" and binary mixtures |81l I82| . 

This method of computing the configurational entropy can be easily implemented 
in the framework of the standard replica method for mean-field disordered systems. 
This has been worked in some detail in |83[ . The starting point is free energy 
at one-step level of replica symmetry breaking F{qo, qi^m), where qo, qi and to 
are the parameters that describe the Parisi matrix |35) in the one-step replica 
symmetry breaking scheme. By expanding the free energy around to = 1 one 
gets, F{qo,qi,m) = FRs{qo) + F'-^'^ {qo, qi){m - 1) -l- 0{{m - 1)^) where qo stands 
for the overlap among replicas belonging to different subboxes and Frs(qo) is the 
free energy in the replica symmetric approximation, i.e., in the limit to — > 1. 
Extremization of Fusiqo) yields qo{P) which inserted into F{qQ,qi,m) allows to 
find Ffis{0),F'^^\(3,qi). The knowledge of these functions fully determines the 
configurational entropy of the system for temperature Trsb < T < (where 
corresponds to the MCT, see Section [6. If) . Indeed the dynamical transition is 
found solving the equations {d/dqi)F^^'> {(3, qi) = {d"^ /dql)F'^^\p, qi) = while the 
static transition Trsb, where the configurational entropy vanishes, follows from the 
solution of gi) — {d/dqi)F^^''{P,qi) — 0. Finally the complexity in the region 

Trsb <T < Tc is given by the value of F(^^(/3, qi) evaluated for qi{P) solution of the 
equation {d/ dqi)F'-^'> {f}, qi) = 0. This approach gives a detailed description of the 
metastable properties in the range Trsb < T < Tc- Below Trsb more sophisticated 
methods are needed to describe the metastable behavior. 

The potential method has been proposed by Franz and Parisi |85| in the framework 
of the replica approach. The starting point in this procedure is to write down the 
partition function of a generic system at temperature T whose configurations C are 
constrained to have an overlap Q{C,Cq) with a reference configuration Cq. The free 
energy of the constrained system is then averaged, using the replica method, over 
the reference configuration Co thcrmalized at a temperature T' in general different 
from T. This yields the potential V{Q,T,T'). For T = T' , and in a given range 
of temperatures, the potential V , as a function of Q, has two local minima. The 
difference of the potential at these two values yields the configurational entropy at that 
temperature. The method has been applied to evaluate the configurational entropy in 
the hypernetted chain approximation usually employed for liquids |86[ I57| . 




F*{T,m), 



— S,{F*,T) 



T 
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5.1.2. Numerical methods Among numerical approaches, Speedy '88' has proposed a 
method that consists in estimating what he defines as the statistical entropy (basically 
identical to the intrastate entropy in the inherent-structures approach discussed 
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Figure 5. Stillinger and Weber decomposition 



in section I4.5|l and comparing it to the thermal or total entropy obtained from 
integration of the specific heat. The difference between the thermal entropy and the 
statistical entropy is the complexity. To compute the statistical entropy the method 
considers different reference configurations representative of an amorphous glass state 
and introduces a coupling term between a reference configuration and the system 
that forces it to stay within a given distance of that reference configuration. By 
progressively slowly changing the intensity of the coupling the energy of the system 
can be evaluated for each value of the coupling. The entropy associated to a particular 
reference state is then estimated by integrating the energy as function of the intensity 
of the couphng. Speedy has applied this approach ^SS, 8211201 to hard sphere systems 
where the center of the hard spheres are tethered to a spherical region with a variable 
diameter that regulates the intensity of the coupling. 

Probably up to now the most powerful method to compute numerically the 
configuration entropy numerically is the one based on the IS formalism. Moreover, due 
to its relatively simple implementation, the IS formalism has become an important 
tool in the numerical analysis of models. For this reason we shall give a more detailed 
presentation. The calculation of IS, summarized in FigurejSl follows directly from the 
definition. First the system is equilibrated at a given temperature T, then starting 
from an equilibrium configuration the system is instantaneously quenched down to 
T = by decreasing the energy along the steepest descent path. The procedure 
is repeated several times starting from uncorrelated equilibrium configurations. In 
this way the IS are identified and their probability distribution can be computed. In 
equilibrium at temperature T the system explores the IS of energy Eis — E with 
probability [see l(TUB|) . ifTHl 



where F{T) is the equilibrium free energy. Then the SW configurational entropy can 
be computed just inverting this relation: 



and using the computed IS probability distribution. If the energy dependence of 
Fy(T, E) can be neglected, then curves for different temperatures can be superimposed 
and the resulting curve is, except for an unknown constant, the SW complexity Sc{E). 
The unknown constant can be fixed either comparing the numerical results with 



V{E, T) = exp [Sc{E) - f3E - (3F^{T, E) + (3 F{T)] 



(119) 



S,{E) = \nV{E,T) + (3E + f3 F,{T,E) - (3 F{T) 



(120) 
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known theoretical predictions, or using the method described below. This method 
works rather well for some disordered spin systems as the Random Orthogonal Model 
(ROM), where the data collapse is rather good for a quite large energy interval |65j . 
The vibrational contribution Fy follows from the motion inside the component. Its 
independence from the energy of IS means that all components are equivalent, i.e., 
have similar shapes. In general this is not the case and its contribution must be 
included. For systems with continuous variables Fy can be calculated at low T in the 
harmonic approximation by expanding the energy about the IS configuration [SI ESI • 
An alternative numerical method consists in computing directly the configura- 
tional entropy as function of temperature. This method is free from unknown constants 
but does not resolve the configurational entropy as function of IS energy. Considering 
Hll()(l and H112(l we have 

F{T) = E{T) -TS{T) 

= EisiT) + E4T) - r5c(T) ^TS4T). (121) 

where Eis{T) is the average energy of IS seen at equilibrium at temperature T. The 
total entropy is then the sum of two contributions 

S{T) ^ S,{T) + Sy{T) (122) 

The first term accounts for the multiplicity of components of energy Eis (T) while the 
second for their "volume" . The SW configurational entropy can then be computed as 
difference between the total and the vibrational entropy. 

The total entropy S can be evaluated via thermodynamic integration of the total 
energy at temperature T from a known reference point: 

AS = S{T) - SiT*) ^ —. (123) 
Jt* ^ 

To compute the vibrational contribution is more difficult, however at low temperature 
the system mainly explores the bottom of the components, near the IS. If the system is 
described by continuous variables then the vibrational contribution can be computed 
in the harmonic approximation by expanding about the IS. This leads to 

AT 

5,(T)~5wm(T)=AA-^ log 

where uji is the (average) frequency of the z-th normal mode and Af the number of 
normal modes. It is possible to refine this approximation by adding terms which 
take into account the basins anharmonicities, however usually these are negligible 
when compared with H124|l j92| . For systems with discrete variables, as for example 
disordered Ising-spin systems, the vibrational contribution can be estimated from the 
T expansion of the Thouless Anderson Palmer entropy, which leads to 

N 

Sy{T) ~ 2(3\h,\ exp(-2/3|/i,|) (125) 

where N is the number of spins, and hi is the local field acting on z-th spin evaluated 
at the IS configuration |93| . 

These methods have been successfully applied to several model systems with both 
continuous variables such as Lennard- Jones glasses EH or discrete variables 
like the ROM (HS] or the SK model [5511^5] . 



hwi{T) 
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Recently in |96| it has been introduced a numerical method to compute directly 
Sc{^,T) within the IS decomposition scheme based on the probabilistic definition of 
the component free energy. The dynamical evolution of the system in equilibrium 
defines a probability measure p-jz over the components. In the case of an ergodic 
dynamics, and assuming that the observation time Tobs is larger than the equilibration 
time, the statistical weight of the component TZ is: 

PniT) = ^ = cxp [-/3F(7^) + /3F(T)] (126) 

'^obs 

where t-jz denote the time spent by system in the component TZ during the total 
observation time Tobs, FiT) the equilibrium free energy and F(TZ) the component free 
energy [See ()75|l ]. The probability to find at temperature T a component with free 
energy equal to is 

n 

= exp[Sc{J',T)-f3J^ + l3F{T)] (127) 

so that 

SciT, T) ^InViT, T) + P[T ~ F{T)] . (128) 

If the number of different components is not too large F{TZ) can be estimated directly 
using (|126|l and the frequency with which a given component TZ appears in a (long) 
simulation at temperature T: 

F{TZ) = -T\n( — ] + F{T). (129) 

\Tohs J 

The equilibrium free energy F{T) can be computed by performing simulations at 
different temperatures and integrating the energy E(T) of the system from T = oo 
down to T: 

(3F{T)^ / d(3' E,^{[3')-S{f3 = Q). (130) 
Jo 

where S{j3 = 0) is the infinite temperature entropy of the system. From the value of 
F{TZ) it is now easy to construct the histogram T'{!F,T) and using H128(l compute 
Sc{J',T). Because the system is equilibrated, in this approach components with 
identical free energy are sampled with the same probability. This differs from the 
previous method where components with the same energy i?ig are assumed to be 
equally probable which is clearly an approximation. The two methods coincide only 
if components have similar volumes so that the component entropy is the same. 

This method has been successfully applied in (SS] to the study of the ROM and 
the SK model, two cases with completely different critical behaviour. In both cases 
the computed potential $ allows for a very precise calculation of critical temperatures 
using relatively small systems giving confidence on Sc{T , T). Moreover the form of $ 
clearly discriminates between the two different types of transitions. 

5. 2. The concept of the effective temperature 

There is the long standing idea that the non-equilibrium regime in aging or driven 
systems can be characterized by the FDR Xa,b{C) ()94|I that has the meaning of a 
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temperature in the thermodynamic sense. This suggestion stems from the observation 
that H94|) can be recast in the foUowing form, 



XA,BiC) ^ 



defining an efi'ective temperature through the relation 



t>s(131) 

CA.B{t,s) = C 



T 



T^^'''\C) = ^—^^ (132) 

As defined in p31|) . p32|l the effective temperature is nothing else than a suitable 
parameter which tells that the QFDT becomes the usual FDT by replacing the bath 
temperature with the effective temperature. In many cases Xa.b{C) < 1 so the 
effective temperature T^g'^'^ is larger than the bath temperature. 

The idea that some concepts of thermodynamic systems can be applied also to 
non-equilibrium systems has been wandering around in the literature for a long time (in 
the context of turbulence see [HZ| or in the context of structural glasses OEHIIMI)- The 
statement that the effective temperature (I132|l has indeed a thermodynamic meaning 
faces some conceptual problems and difficulties not found in equilibrium theory. T^^'^^ 
should satisfy the following properties, 

• Observable independence. T^g'^\c) must be independent of the observables 
A, B used to construct correlations and responses. If this is not always possible, 
as the present numerical evidence suggests, at least one would still like to know 
beforehand which set of "good" observables endow H132|) with a physical meaning. 
These observables have received the name of neutral observables and have been 
discussed in section 14.41 From many perspectives, this condition appears quite 
strong. It could be relaxed by only requiring independency of Tch upon the 
measured observable A for a given perturbation B (rather than on both A, B). 

• Zero**^ law. If the slow degrees of freedom of a system described by effective 
temperature Tcff(C) (we have dropped the A,B dependence) are put in contact 
with a thermal bath at temperature T, the net heat flow between the system and 
the bath should vanish only if Toff(C) = T, where C determines the relevant 
timescale (as described in section 14. 1|) at which the thermal bath, acting as 
a thermometer, responds. The same conclusion must hold between two glassy 
systems described by two effective temperatures T^g{C), T^g\c). After putting 
them in contact the net heat current between them, at the relevant timescale 
defined by the correlation C, vanishes only if T^g\c) = T^g\c). This definition, 
apparently reasonable, encounters some difficulties that we will describe below. 
In particular, systems with identical effective temperatures T^g\c) — Tl^g\c') 
but at different timescales (C ^ C") cannot be in mutual equilibrium. 

• Existence of a non-equilibrium measure. The zero*^ law carries associated a 
maximum principle. In standard thermodynamics the zero*** law establishes that 
after putting in contact two systems at different temperatures the global system 
reaches a stationary state with a unique temperature. This stationary state is the 
one that maximizes the global entropy of the compound system compatible with 
a given total energy. Moreover fluctuations around this maximum entropy state 
are ruled by the temperature. By the same token, the aging state of relaxational 
systems and the stationary state of driven systems must exhibit some fluctuations 
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or deviations around the aging (or driven) state that are described by the effective 
temperature T'^^{C). The full characterization of these fluctuations is presently 
unknown. 

In Section 14.21 we tried to fortify the idea that a thermodynamic description 
is indeed possible and that the effective temperature shares some properties of 
thermodynamic temperatures. Cugliandolo, Kurchan and Peliti |l()()j have emphasized 
these aspects showing that the eS'ective temperature can be defined only regarding 
the timescale under consideration. They considered a small thermometer that can be 
mimicked by a single harmonic oscillator of frequency to that is put in contact with 
the original system. For definitiveness let us consider that x denotes the oscillator 
coordinate and 0{y) an observable of the system (described by the variable y) to which 
the oscillator is coupled by an interaction term, —exO{y) where e is the intensity of 
the coupling. If e is small enough, then the interaction between the oscillator and 
the system can be treated within the linear-response theory and the energy of the 
oscillator evaluated in the stationary state. The equipartition theorem relates this 
energy to the temperature measured by the oscillator acting as a thermometer. In 
aging systems, the effective temperature is given by the FDT in the frequency domain 
(also called Nyquist formula), 



where S'o(w, t^) is the power spectrum or correlation {0{t)0{tw)) expressed in Fourier 
space (see H275|) ') and Xo('^'^tu) ^he corresponding out-of-phasc susceptibility. A 
similar expression is valid for the stationary non-equilibrium state of driven systems, 
however because TTI holds the t^ dependence in H133|) drops off. The connection 
between H132f) and (|133|l appears when translating the meaning of uj and t^ into the 
many timescales scenario. According to that lo corresponds to l/{t — tw) and therefore 
we can define C* = C{t, t^) = C{tji, -|- l/w, tw)- This means that a thermometer put in 
contact with the system at time t^ and responding at a given frequency w measures 
the temperature Teff(C*) = T/X{C*). Equivalently, to measure the temperature 
Tcff(C*) = T/X{C*) in an aging system a thermometer responding to a time scale 
t* — 1/lu with C* = C{t^ + l/uj,t^) should be used. In aging systems with two- 
timescales (such as structural glasses) characterized by the full aging t/t^^, in order to 
measure the effective temperature associated to the slow process, the frequency of the 
thermometer must he lo ^ 1 jt^ . The thermometer must respond in a timescale of the 
order of the waiting time!!. In this scenario, effective temperatures can be extremely 
difficult to measure and this raises the question about their true meaning as the system 
drifts away from that state in the time required for a single measurement. To cope 
with this problem it has been proposed IIOII that an ensemble of small thermometers 
is needed for the measurement. However, this does not solve the problem of how to 
measure, using this procedure, the effective temperature of a vitrified piece of glass 
quenched one thousand years ago. These difficulties are inherent to aging systems 
but not necessarily in driven systems that reach a stationary TTI state. For these 
reasons, experimental measurements of FDT violations and effective temperatures 
could be more suitable in driven rather than aging systems (concerning experiments 
see discussion in section IH)). 

Zero"' law aspects of the effective temperature have been considered in |1(J(JI ll(J2) 
within the spherical p-spin model. It has been shown that in the large t^ ^ oo 
limit and low frequency limit w — > there are only two possible Ansatz solutions 



■K LoSo(y>,tw) 



(133) 



2 X^(w,i») 
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that close the dynamical equations in the aging state: either the two systems remain 
uncoupled with different effective temperatures or they thermalize and reach a common 
temperature. These results have been endorsed by a systematic study of these 
coupled solutions in the framework of the oscillator (OSC) model |l(J3j . The OSC 
model is characterized by a single timescale corresponding to the slow process at 
zero temperature (see section ffi.5.1|l . It has been shown that, in the presence of an 
interaction term in the Hamiltonian that describes a compound system formed by two 
ensembles of the OSC model, the dynamics behaves in two different ways: either the 
effective temperatures equalize or they are different as if the systems were uncoupled. 
Dynamics remains always uncoupled (independently of the coupling intensity) if 
dynamics is sequential on both systems, i.e., the updating is done sequentially over 

the two ensembles 1 — > 2 — > 1 ^ 2 In this case the effective temperatures of 

the two systems differ, each one corresponding to that of the non-interacting OSC 
ensemble. If dynamics is parallel, i.e., updating is done simultaneously over the two 
OSC ensembles 1 + 2^ 1 + 2— > 1 + 2 .... and so on, then the effective temperatures 
of both models coincide even for a zero value of the coupling e in the Hamiltonian. 
This result shows that, in general, two glassy systems interacting through a coupling 
term in the Hamiltonian do not necessarily reach the same effective temperature on 
timescales of the order of the waiting time. 

In all these studies the same question remains always unanswered: why fast 
and slow degrees of freedom decouple into different effective temperatures (in two- 
timescales systems, one is the bath temperature, the other the (higher) effective 
temperature). A necessary condition is that the relaxation rate of the energy 
(or entropy production) decays to zero slow enough jl()4| . To better understand 
this question, in |105j the thermal current between the oscillator model and a 
thermometer was analyzed. There it was shown that the measured temperature 
Tm, which makes the net current flow between system and the thermometer vanish, 
coincides with the effective temperature ii cot ~ 1. However, in the limit ujt <ti I 
the measured temperature is much smaller than the effective temperature, while in 
the other extreme U!t ^ 1 the thermometer measures tolerably well the effective 
temperature. Discrepancies between the measured and the effective temperature have 
been also reported in the SK model in the presence of asymmetric interactions as 
an example of a driven system [lOfij . Another important aspect is that the zero*'' 
law is hardly effective as transport coefficients (such as the thermal conductivity) 
are exceedingly small, in agreement with the uncoupling of degrees of freedom 
occurring in glassy systems with many timescales. From another point of view 
Nieuwenhuizen |107[ 1491 1108| has formulated a theory to describe the aging regime 
of glassy systems with two-timescales assuming, right from the beginning, that the 
effective temperature is indeed associated to a thermodynamic potential. In this 
formulation, the configurational entropy is an extensive thermodynamic potential 
conjugated to the effective temperature. The first law of thermodynamics, that 
expresses energy conservation reads dE = dW + TdScq + T^ffdSc where S = S^q + Sc 
is the total entropy that receives contributions from the equilibrium (or intrastate) 
entropy and the configurational entropy. The configurational entropy and the effective 
temperature manifest in the reported experimental failure of the second Ehrenfest 
relation (while the first is automatically satisfied by construction) leading to values of 
the Prigogine-Defay ratio larger than 1 jl09| . 



Violatior 



41 




10 10° lo' itf lo"' lo" , lo' lo' 

t/x t/x 



Figure 6. Left panel: Relative difference between tlie measured temperature 
by a thermometer responding in a timescale r = l/ui and tlie glassy 
OSC model aged at t, plotted against t/r for different values of r = 
10'',105,106,10'^,108,109,10i0, IQil (from bottom to top). The measured 
temperature coincides with T^ff for t/r ~ 0{l). In the inset, the same relative 
difference but plotted against t for t = 10'^, 10*. Right panel: Measured 
temperature plotted as a function of t/r. From I105| . 

5. 3. The Edwards ' measure for granular materials 

Granular materials are systems made of a large number of individual grains such 
as sands or powders. At first sight granular systems look quite different from 
thermodynamical systems since for example they interact mainly through frictional 
forces and hence the energy is not conserved. Moreover power is supplied to them 
by tapping, shearing or shaking, all mechanisms quite different from thermal contact 
with a thermal bath. However, despite these differences, granular materials share with 
thermal systems the property that their properties are reproducible given the same 
set of extensive operations, i.e., operations acting upon the system as a whole and not 
on individual grains. For example if some sand is poured uniformly and at low density 
into a container one expects to have a sand of a certain reproducible density. Based on 
these facts it is then reasonable to hypothesize that granular systems can be described 
at a macroscopic level by a small number of parameter analogous, e.g., to temperature 
or pressure, using some ideas of statistical mechanics |ll(J[llllllll2[[TT5llll4| . 

The most important variable describing the state of a granular system is its 
density, or equivalently its volume V . The thermal energy of a granular system at 
room temperature is indeed negligible. The volume V is the actual volume occupied 
by the system, for example measured by the position of a piston in the container, and 
hence it depends on the configurations (position and orientation) of the grains. In 
principle other variables should be necessary to describe the state of the system, but if 
we assume that the grains are rigid then all other microscopic details can be neglected. 
Thus the only valid configurations of grains are stable arrangements where grains can 



Violation of FDT in glassy systems 



42 



remain at rest under the influence of confining forces, and with no overlap among 
them. We are then led to a statistical description of all possible stable configurations 
of grains in the real space, i.e., to a configurational statistical mechanical theory for 
the random packing of grains. Since in a stable configuration grains cannot move these 
configurations are also called blocked states. 

The next step to develop a statistical mechanical description for granular systems 
consists in the introduction of a volume function W which plays the role of 

energy in statistical mechanics. The function W specifies the volume of the system in 
term of the configuration of grains: V = W{q), where q denotes all the grain positions 
and orientations in a blocked state. Under the assumption that for a given volume V 
all configurations for which V = W{q) are equally probable a statistical description 
can be developed through a process completely analogous to that of conventional 
statistical mechanics. It is then possible to introduce a microcanonical ensemble with 
distribution function 

e''^ 6{V -W{q)) (134) 

where 5* = S{V) is the entropy defined as usual in terms of the total number of blocked 
states 

n = J dqS{V ~W{q)), S^lnn (135) 

Similar to its thermodynamic counterpart, the entropy S is an extensive quantity as 
can be seen, for example, in simple toy models |lllllll,')] . The measure (jl34ll - H135(l is 
usually called Edwards' measure for granular systems. 

To define a canonical ensemble it is necessary to define a parameter analogous 
to temperature which characterizes the state of the system. This parameter is the 
compactivity X defined as 

1 dS , , 

X - BV ^'''^ 
and thus the partition function is 

Z = e-^(^)/^ ^ I dqe-^('')/^ (137) 

where potential Y{X) is called the effective volume and plays the role of a free-energy: 

Y{X)^V{X)-XS (138) 

The analogy can be pushed forward and many other relations similar to that of 
conventional statistical mechanics can be derived. We shall not pursuit this here, 
however, before concluding we shall spend some more words on the compactivity X . 
The compactivity X measures the packing of grains, indeed from its definition it is 
clear that it may be interpreted as being characteristic of the number of possible ways 
of arranging the grains of a system by changing the volume in an amount A.V, the 
change in entropy being equal to AS. Consequently the two limits of X are and 
oo corresponding to the most compact the least compact arrangements, respectively. 
The compactivity X also describes the balance between the tendency of the system 
to increase or decrease its volume and the tendency to increase or reduce its entropy. 

At first sight blocked states in granular systems resemble IS discussed for glasses. 
Indeed a connection between the two can be drawn [1161 11171 1118L ET] introducing a 
"tapping" dynamics for glasses, i.e., a dynamics in which each tap consists in raising 
the temperature and, after a short time, quenching it to zero. Similar to what has 
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been done for other glassy systems one can try to describe the dynamics of the slow 
degrees of freedom through an effective temperature defined from the FDR |1(J(JI I119| . 
At the mean field level this temperature turns out to coincide with the Edwards 
compactivity, which is related to the derivative of the entropy of blocked configurations 
of a given density. The Edwards ensemble immediately leads to the definition of an 
entropy SEdwif) as the logarithm of the number of blocked configurations of given p. 
The soundness of the Edwards approach encounters difficulties reminiscent of those 
present in the IS approach. At the present the correspondence between the Edwards' 
construction and the long-time slow dynamics for non mean-field models can only be 
checked "a posteriori" and it is presently unknown how to derive it from first principles. 

6. QFDT from exactly solvable models 

In the structural glass problem the spatial randomness is self-generated rather than 
put in by hand as in random spin glass models. This suggest that there should be 
a connection also with frustrated but regular models. In the last years several spin 
models, both with and without randomness, displaying structural glass transition like 
properties have been found. Interestingly, some of them can be solved in closed form, 
offering an important tool for understanding the glass transition. 

In this section we shall summarize the main results on violation of FDT obtained 
from the exact solution of some spin models. 

6.1. The Mode- Coupling Theory 

A model Hamiltonian or an effective Lagrangian capable of describing relaxation 
processes in supercooled liquids and structural glasses is difficult to obtain. Early 
studies based on both dynamical mode-coupling theories or equilibrium density- 
functional theories suggested that there may be a close connection with mean-field 
spin-glass models I12()j . They thus provide a set of microscopical models where glassy 
dynamics can be studied analytically. The basic simplification occurring in mean-field 
models is that after averaging over the disorder and making the number of spins very 
large {N oo) one is left with a closed set of equations for the two-times correlation 
and response functions. Above the a critical temperature Tc those equations admit 
a TTI solution satisfying the FDT. In this regime they are basically equivalent to 
the schematic mode-coupling equations introduced by Leutheusser, Gotze and others 
[ 1211 1122L I123| as a model for the ideal glass transition. Below Tc the ergodicity 
is broken and the FDT is violated. This is signaled by the appearance of a finite 
Edwards- Anderson order parameter qea at Tc. 

The fundamental quantities in the dynamical mode-coupling theory (MCT) are 
the local particle density correlation functions {Sp{x, t) Sp{x)) where Sp{x) = p{x)—po, 
with 



JV 




(139) 



the local particle density, and po the uniform fluid density 
{p{x)) = pq (homogeneous state) 



(140) 



where the angular brackets denote an ensemble average. In the glassy phase the 
system is trapped into metastable states with nonuniform (average) local density field 
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{p{x)) 7^ pq and the density-fluctuations correlation functions do not decay to zero for 
t ^ oo: 

\im{5pix,t)Sp{x)) j^O. (141) 

t — ^oo 

The complete mode-couphng theories lead to the time-evolution equations for the 
normalized correlation functions 

{Sp{q,t)* Sp{q,0)) 

'^'^'^ = nJ, 

where Sq = {\S p{q)\'^) / N is the static structure factor, and p{q) are the Fourier 
components of the density field p{x) 

N 

p{q) — / exp(— i<7 • a;) p(a;) da; = ^ exp(— iq • aji). (143) 

i—l 

The basic idea of MCT is to derive the equations of motion for the slow relaxing modes 
integrating out all fast modes. This leads to a set of self-consistent equations involving 
only slow modes variables in which all informations from fast modes are buried into 
density-fluctuations memory kernels of the form 

Mq{t) = iVq + ^lmq{t) (144) 

where Vq is a (white-noise) frictional term arising from fast modes, VLq > Q gives 
the frequency or time-scale of microscopic motion and vriqit) accounts for slow modes 
couplings arising from the integration of the fast modes. The general form of the MCT 
equations is (123| 

/■* 

^2 J, i,\ , ,, a J, , n2 . 



dtcl)qit)+iyqdtMt)+^tMt)+^q ds TO, (< - s) 9, 0g (s) = (145) 

Jo 

which must be solved with initial conditions: 

0,(i = O) = l, dt^t = 0)=0. (146) 

The fundamental mechanism for the glass transition in the MCT is the feedback 
between slow density fluctuations expressed through mq{t). The solution of these 
equations is a formidable task since the kernel mq{t) involves higher-order correlations 
between density-fluctuation modes. Therefore when these theories are implemented 
approximations are generally made. The simplest approximation consists of replacing 
the average of products with products of averages to obtain a set of closed equations. 
This is some sort of mean-field approximation. Indeed within this scheme the memory- 
kernel mq(t) can be expressed as a functional of the 

mq{t)=J^q[V,{c^q}]. (147) 

with some vertex functions V. Despite this rather strong approximation, similar to a 
mean-field approach, the theory contains the basic mechanism of the glass transition. 
We note that due to this approximation, the MCT is not capable of describing 
activated process, in the same way they cannot be discussed within mean-field theories. 
Therefore the appearance of activated-process dominated regimes is signaled in this 
theory by the divergence of some quantities. Activated process could in principle be 
included as perturbative terms, however consistent theories which account for them 
are not yet available. 
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The main properties of the MCT can best be seen using a simphfied version of the 
theory called schematic mode- coupling theory in which only one relaxation function is 
considered [11111112111231 : 

dt(t){t) + V dt^{t) + (f>{t) +0^ I ds m{t - s) ds(f>is) = 0. (148) 

The simplest model describing an idealized structural glass transition is the one 
specified by the two coupling constants (wi,f2): 

m{t) = VI (l){t) + V2 (l){t)^ . (149) 

This theory predicts a transition from an ergodic liquid phase, where 4'{t — s- oo) — > 0, 
to a glass phase, where the ergodicity is broken and (j){t oo) ^ / > 0, as the 
parameter (wi,W2) are varied. Depending on the values of {vi,V2) the nature of the 
transition can be either continuous (type A) with / growing continuously from zero 
or discontinuous (type B) with / jumping from zero to a finite value as the transition 
line is crossed. 

6.2. Disordered spin- glass models 

Mean-field spin-glass models can be classified into two broad classes depending on the 
value of the Edwards- Anderson parameter {qea) at the transition (for the structural 
glass transition this can be identified with the long-time limit of the density correlation 
functions). The first class, called discontinuous models, includes models for which a 
finite gEA appears discontinuously at Tc. The prototype model in this class, which 
we discuss here, is the spherical p-spin model. Other models included in this family 
are. Potts-glasses with more than four components |124lll^ . quadrupolar glass models 
[ElinSI and p-spin interaction spin-glass models [ESI EH EHl ESI Ej ■ The second 
class includes models, such as the Sherrington-Kirkpatric (SK) model |129l IliiO] . for 
which gEA starts continuously from at Tc. Those models are termed continuous 
models. 

6.2.1. p-spin spherical model. Among the mean-field spin-glass model with a 
discontinuous spin glass transition an important role has been played by the spherical 
p-spin spin-glass model. Spin-glass models with multispin interactions were first 
considered in the eighties for both Ising [1261 11271 1128| as well as soft spins 021 ■ 
However, while the static properties could be computed for arbitrary values of p, 
dynamical properties were limited to values of p close to 2 021 An important step 
forward came with the introduction of the spherical p-spin spin-glass model |i:illli:^2j . 
since its statics and dynamics can be solved in closed form for any value of p. 

The spherical p-spin interaction spin-glass model is defined by the Hamiltonian 

N 

H^- ^ Jii...ipO-ti ■ ■ -a-i^ - h^ai (150) 

l<ji<-<jp<A' j=l 

where h is an external field, which in the following we shall take equal to zero for 
simplicity. It describes a system of N continuous spins ai interacting via randomly 
quenched infinite range p-spin interactions Ji^...ip which are taken to be Gaussian with 
zero mean and variance 

IJ^^^T^^ . l<ii<---<ip<N . (151) 
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The overbar stands for the average over the couphngs. The scahng with N is chosen 
such that there is a well defined thermodynamic limit N ^ oo. The spins can vary 
continuously from — oo to +00, but are subject to the global spherical constraint 

N 

J2a^ = N (152) 

which must be satisfied at any time. The relaxational dynamics for ai{t) is given by 
the set of Langevin equations 

T^^dMt) = ~r{t) a,{t) - ^ + r^,{t) (153) 

The kinetic coefficient Fq sets the time scale of the microscopic dynamics, and 
will be henceforth set to 1 without loss of generality, while (3 — 1/T. The last term in 
H153|l rji{t) is a Gaussian random field with zero mean and variance 

{V^{t)vAt'))='2^o'S^,S{t-t') (154) 

representing the effects of thermal noise. The average over thermal noise is denoted 
as usual by angular brackets (■■•). The first term in p53(l enforces the spherical 
constraint at any time, and must be fixed self-consistently. In the mean-field limit 
N —> 00 the sample-averaged dynamics is entirely described by the evolution of the 
two-times correlation and linear response functions 



C{t,s) = {a,{t)a,{s)) (155) 

(156) 



6hi{s) 



hi=0 



which vanishes for t < s due to causality. The dynamical equations for C and R are 
obtained from l|153|l )132| through standard functional methods |134llinilH] . At high 
temperatures the system is ergodic, thus for initial time ti —00 both TTI and FDT 
hold. Using the FDT relation the equation for C{t — s) reads |i;-{2j 

drC{T)+C{T)+ f dT'm{T~T')dr'C{T')^0 (157) 



JO 

where t = t — s and 

m{T) = nC{Tf-^ (158) 

which has the same structure of the schematic MCT equation (|148|l . The correlation 
always decays to zero for large t. However slightly above Tc is develops a plateau at 
C ~ qea before eventually decaying to zero, see Figured The length Tp{T) of the 
plateau increases as a power of T — Tc and diverges at Tc. Near the plateau one finds 
that: 

CiT)^qEA+CaT-\ C>qEA (159) 

C{t) ^ QEA - Cb t\ C < (ZEA (160) 

where the exponent a and b are related by: 

r^(l-a) ^ r^(l + 6) ^ (p-2)(l-gEA) 

r(l-2a) r(l + 26) 2gEA ^ ^ 

The plateau length scale Tp(T) sets the equilibration time-scale. Therefore, as the 
temperature is lowered down to Tc the system undergoes a transition since the length of 




Figure 7. C'(t) as function of time for the spherical p-spin model with p = 3. 
Left to right: T = 0.7, 0.65, 0.62, 0.615, 0.613, 0.612 and 0.61. The horizontal 
line is gEA(T = Tc) ~ 0.509034 and Tc = 0.612372 



the plateau diverges and the correlation fails to decay to zero. In the low-temperature 
spin-glass phase (T < Tc) the system cannot equilibrate and the ergodicity is broken 
so the state of the system may depend on its initial state. In this scenario it is clear 
that even a mean-field theory can be highly non-trivial. To discuss the dynamics 
in the low-temperature phase it is convenient to have an understanding of the low- 
temperature structure of the phase-space. A standard method to deal with such a 
problem is the so called replica trick, see e.g. jciSj . The breaking of ergodicity results 
in a breaking of the permutation symmetry between replicas |135) . The general form 
of the breaking is, however, not known. For p-spin models it has been found that 
the solution is given by the so called one-step replica symmetry breaking (IRSB) 
form [1271 11281 144L [TM] , Physically this means that the phase space is broken into 
equivalent ergodic components separated by infinite (for N ^ oo) barriers and the 
equilibrium solution is described by only three parameters: the overlap between two 
different ergodic components, the overlap inside an ergodic component (the Edwards- 
Anderson order parameter) and the probability that two different replicas will be found 
in the same ergodic component. For the p-spin spherical model this solution is valid 
everywhere in the low-temperature phase |131| . Ising-spin p-spin models present an 
even lower temperature phase with a more complex structure described by an infinite 
step replica symmetry breaking (oo-RSB) jl28| . A similar scenario is also found in 
Potts glasses with more than four components |124| . 

Besides the replica approach, which gives essentially information on the lowest- 
lying states which mostly contribute to the equilibrium measure, a good understanding 
of the landscape topology is given by the TAP approach. Using these methods we have 
now a rather good knowledge of the landscape [ElIZlESlinZIIISHlinillini ■ Roughly 
speaking the picture that emerges for the p-spin spherical model is that equilibrium 
states, whose energy and free energy difference is 0(1), are separated by 0{N) barriers. 
As the temperature is changed the solutions neither merge nor bifurcate, and their 
free energy smoothly changes. The TAP equations for the p-spin spherical model have 
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non-trivial solutions in the free energy range (firsbi ^thr)j where the threshold free 
energy Fthr is larger that the equilibrium free energy -Firsb by an 0{N) quantity, the 
difference being the complexity, see section 15.11 which is maximal for F = F^hr f^nd 
vanishes when F = Fi^sh ■ Below the threshold the equilibrium states are local minima 
of the TAP free energy separated by 0{N) barriers, while above the threshold there 
are no minima. 

To have a meaningful investigation of the non-ergodic phase some regularization 
scheme of the dynamics on a very long time-scale is required. One possibility is to 
refer to a large finite system. The finiteness of N guarantees ergodicity by allowing 
the penetration of barriers whose height would diverge for N —^ oo limit, and the 
system can equilibrate. This approach was first proposed by Sompolinsky for the 
Sherrington-Kirkpatrick model |141j . In the Sompolinsky's scheme TTI holds and 
FDT is satisfied up to some time-scale to, which diverges as iV ^ oo, related to the 
(free)energy barrier crossing but it is violated for time-scales larger than to where it is 
replaced by a modified form called quasi fluctuation-dissipation theorem (QFDT). See 
section 15. 2. 21 for a more detailed discussion of the Sompolinsky's scheme. In order to 
consider the motion on the two different time regime (t <ti to and t 3> to) one writes 

C{T)=Ci{T)+a{0: R(t) = Ri{t) + ^ RoiO (162) 

to 

where ^ = r/to. The functions C\{t) and Ri{t) describe the motion in a single ergodic 
component and vary on time-scales <C to, while the functions Co(^) and i?o(C) describe 
the motion among different minima and hence vary on time-scales ^ to. Continuity 
imposes 

Ci(t = 0) = l-QEA, Ci(t^oo)-0 (163) 

Coi( = 0)=qEA: Co(e-^oo) = (164) 

thus Co describes the slow decay of correlation function from ^ea to zero. The initial 
conditions for Rq,Ri are obtained from the FDT, the QFDT and the continuity 
condition for C. 

On time-scales r ^ to FDT holds, so Ci(t) obeys a dynamical equations similar 
to (|157|) . The equations for Co(^) and i?o(C) ^-re more complex, however it can be 
proved |132| that they are related by the QFDT relation 

RoiO^-PxOiOd^CoiO (165) 
with < a: < 1 given by: 

Within the replica formalism the parameter x corresponds to the location of the 
discontinuity in the order parameter q{x). However at difference with the static 
calculation, where x is fixed by the requirement of stationarity of the replica free 
energy F with respect to x, the dynamical calculation requires dxF be maximal 
(marginal condition) jl32| . This condition is equivalent to the condition of a maximal 
configurational entropy |42[ 1731 [77| , so that the dynamics is dominated by the states 
with the largest degeneracy (threshold states). 

The Sompolinsky approach has several similarities with the static approach, and 
indeed in the static limit it correctly reproduces the static results obtained within the 
Parisi scheme. However it suffers of some problems which are difficult to amend since 
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it would require the detailed knowledge of the dynamics for a large but finite system. 
One among the most serious inconsistencies of the Sompolinsky dynamics is that 
FDT violations, as given in (|l(j5|l . satisfy TTI. This is untenable in the aging regime 
of purely relaxational systems (although not necessarily in driven systems) where TTI 
is clearly violated. It must be noted though, that the Sompolinsky approach was never 
proposed to explain aging, since it is eliminated from the theory at the beginning. 

For this reason there have been various attempts to amend the solution. One 
possibility, proposed by Horner ^21 11421 114HI \1'A'2\ , consists in cooling the system, at 
finite cooling rate the system from T > Tc to T < T^. This introduces a regularization 
time scale, the inverse of cooling rate, which is eventually sent to infinity at the end. 
Another approach |32| , consist in making the disorder time-dependent hence restoring 
ergodicity on time-scales larger than the disorder typical time-scale. 

The above methods assume in one way or another equilibrium, thus cannot 
describe non-equilibrium properties typical of glasses such as aging, see section ETI To 
tackle them a different scheme has been proposed by Cugliandolo and Kurchan |52j . 
The main difference lies in that the thermodynamic limit is taken before any large 
time limit, including the initial limit ti —oo. At difference with the Sompolinsky's 
approach this is a non-equilibrium scheme since the system evolves from a non- 
equilibrated initial configuration. As the system evolves in time the dynamical free- 
energy density decreases, and the systems explores an ever decreasing portion of phase 
space. The weak ergodicity breaking ^Uj describes this non-equilibrium regime before 
equilibrium is reached. In this scenario the important parameter is the waiting time 
tyj, i.e., the time elapsed since the quench into the low-temperature phase. For longer 
waiting times the system can explore deeper minima becoming less susceptible to 
external perturbations an hence ages. The weak ergodicity breaking scenario can be 
summarized in the following assumptions. 

(i) After any time tyj the system continues to drift away and asymptotically reaches 
the maximum allowable distance. Thus the correlation functions satisfy 

drC{T + t^,t^)<Q, d,C{t,s)>Q, it>s) (167) 

and in the absence of external magnetic fields 

lim C(t + tw^t^) = for any fixed t^ (168) 

r — *oo 

(ii) The response to a constant small magnetic field applied from s = Q to s = tw - 
i.e., the TRM M'^™^{t,tw)- decays to zero after long enough times 

lim / dt' R{t,t') ^ Vnn M'^^^{t,t^) = Q (169) 

t — 'OO Jq t — *oo 

for any fixed i^,. 

(iii) The evolution of the two-times correlation function presents two distinct regimes. 
After a long time s, but r — i — s small, the correlation decay from 1 at equal 
time to a plateau value (/ea defined as: 

9ea = lini lim C{t + s,s). (170) 

r — >oo s — >oo 

This fast decay corresponds to a fast relaxation toward a local minimum. In 
this time-sector the system behaves as if it were in a local equilibrium, and both 
TTI and FDT hold. The value of gEA measures the size of the local minima or, 
equivalcntly, the width of the channel through which the system evolves. This 
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fast relaxation is followed by a slow decay of C below ^ea and the exploration of 
different minima. Since the depth of minima increases with time, C decays from 
the plateau in a manner that depends on both s and r. To show the two processes 
(i.e. the contribution from fast and slow motion) the response and correlation 
functions are split in two different terms, in a fashion similar to that used in the 
Sompolinsky scheme: 

C{t,s) = Cst{t-s) + C,g{t,s), R{t,s)^R,t{t-s) + R,g{t,s) (171) 

with 

Cst(i-s = 0) = 1-gEA, Cst{t - s ^ (x^) = (172) 

Cag(i, t) = QEA, lim Cag(i, s) = 0. (173) 

t~*oo 

The assumption of local equilibrium implies that FDT is satisfied by the fast motion: 
Rstit -s)^l3e{t- s) d, C,t{t ~ s). (174) 
On long time-scales, however, FDT is violated and replaced by [cfr. eq (|165|l ] 

i?ag(i,s) =/3X[Cag(t,s)]0(t-s)a.Cag(t,s) (175) 

with the Ansatz that X depends on times only through Cag. The two forms of FDT 
can conveniently condensed into one extending the definition of X as X(z) = 1 for 
9EA < 2 < 1 • Then we can write 

R<yt, s) ^ (3X[C{t, s)] 6{t - s) ds C{t, s) (176) 

where C and R are the full correlation and response functions. For the FDT part the 
MCT-like equations (|157|) are recovered. To derive the equations for the non-FDT 
part the time derivatives are neglected since Cag and i?ag are slow varying functions: 

5tCag(i,s)~9,ag(t,s)~0 for large t,s. (177) 

As a consequence the solutions are time-reparametrization invariant, i.e., the solutions 
are invariant under the transformation 

Cag(t,s) ^ C^^{h{t),h{s)), R^git,s) ^ 

where h{t) is an arbitrary (well-behaved) function. The full dynamical solution 
obviously does not have such an invariance. This ambiguity stems from the fact 
that equations are only the first order equations of an asymptotic expansion. If higher 
order terms are included the ambiguity is removed, however we shall not discuss this 
problem. Motivated by the fact that the relevant time scale in slow relaxation motion 
is given by the waiting time, one looks for a solution of the form |52|: 

ag(t,s)^ag(/i(i)/Ms)). (179) 

The selection of the correct function h{t) is still an open problem that requires the 
matching of the short and long-time regimes. Numerical solution of the dynamical 
equations 52] suggest a power law h{t) ~ t^. 

An analysis of the correlation function near the plateau qea |144| , similar to that 
done in the high-temperature phase |132| . reveals the following scenario. As found in 
the high-temperature phase, the decay of the correlation function to the plateau ^ea 
is a power law with a temperature dependent exponent. The subsequent departure 



dhjs) 
ds 



i?ag(/l(t),/i(s)) 



(178) 
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from the plateau is still a power law with another temperature dependent exponent, 
but at difference with the high-temperature phase it is also t,„ dependent: 

C(T + t^,t^) - 9EA + Car-°, C>qEA (180) 

C{t + t^,,t^,) qEA- Cb (^Y-^ ' C* ^ 9EA (181) 

with = [dlnh{tw) /dtt^]~^ an effective waiting time. The exponent a and b are 
related: 

r^(i - «) _ ^ r^(i + b) _ iP~ 2)(i - gEA) , . 

r(l - 2a) r(l + 26) 2(jEA ^ ^ 

where ^ea is determined by the marginal condition 

1 

(1 - ^ea) 



^{p-l)qP,f = j--—^ (183) 



^^^^^ (p-2)(l-gEA) (184) 
9ea 

and X{C) ~ 1 otherwise. 

Here we did not consider the case of non-zero external field, this was studied in 
[ 1311 1132| . Another extension of the p-spin spherical model is the case of multiple 
phases treated in |145l 11461 1147L I148| . It is interesting to note that depending on the 
degree of non-linearity of the interaction three different scenarios for the transition can 
be observed |147j . Finally, it is worth to note that the p — 2 version of the model has 
the property of being solvable even for intermediate time-scales [H??! 11^(11 IT^ II II 521 ITKH) 
or finite sizes |154) . 



6.2.2. The Sherrington- Kirkpatrick model. The Sherrington-Kirkpatrick (SK) model 
[ 1291 1130j belongs to the class of continuous spin-glass models characterized by a low- 
temperature spin-glass (SG) phase of oo-RSB type with a continuous order parameter 
function q{x). The transition to the SG phase is continuous with a q{x) which grows 
continuously from zero as T is lowered below Tc- Other models in this class are, e.g., 
the case of a particle in a long-range correlated random potential |155| or spherical 
models with mixture of p = 2 and p > 3 interactions |145l I147| . For T > Tc and 
large times the correlation function decays to zero in the absence of external fields. 
However, at difference with the discontinuous SG models, it does not exhibit a plateau 
at gEA for T slight above T^. The SK model without external fields is defined by the 
Hamiltonian 

H = - ^ <J,aj (185) 

i<j 

where the interaction J,j are independent quenched Gaussian variables of zero mean 
and variance 

(186) 

As done previously J can be set equal to 1 by rescaling the temperature. The spin 
variables can be either Ising spins (cr = ±1) or soft spins in which case and extra term 
may be added to the Hamiltonian to control fluctuations. The first choice is used in 
static calculations, while the second is preferred in a dynamical approach. 
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The oo-RSB structure of the SG phase reflects a complete different topology 
of the phase space. Again below the phase space is broken into a large number 
(exponentially in N) states, but now the overlap between states can take any value 
from to qEA- The equilibrium states are organized in an ultrametric fashion with 
non-extensive barriers - 0{N") with a ~ 1/3 - between them. The TAP analysis 
shows that the TAP solutions tends to split as the temperature is lowered is a fashion 
similar to a second-order transition. Moreover the spectrum of the Hessian matrix of 
the solution extends down to zero, leaving the possibility of finite free-energy barriers. 
All these facts lead to a dynamical scenario quite different from that discussed above 
in section 15. 2. II 

The relaxational dynamics of the soft-spin version is given by Langevin equations 
similar to (|153|I - H154(I and the self-consistent dynamical equations for the two-times 
correlation and response functions were first derived by Sompolinsky and Zippelius 
[ 1561 1134j . These are more involved than those of the spherical case because spin 
variables cannot be integrated away and hence will not be reported here. However, 
near the critical point the dynamical equations can be written in the MCT form (|157|l , 
l|147|l with suitable vi and V2 |157| . The two parameters vi and V2 are not independent 
and their particular form leaves only the possibility for the type A transition with gsA 
growing continuously from zero at Tc. 

Above the critical temperature Tc no ergodicity breaking occurs, and the solution 
are TTI and satisfy FDT. Below Tc the ergodicity is broken and some scheme must be 
adopted for the long-time dynamics. In the Sompolinsky approach |141| . as discussed 
in the preceding section, one assumes that the initial time is sent to — oo keeping the 
system size large but finite so that the system can equilibrate, and only then N is sent 
to oo. Two-time quantities as correlations and response are then trivially TTI, but 
FDT may not be necessarily satisfied for the infinite size system due to the emergence 
of infinitely high (free)energy barriers for N ^ oo where freezing of some degrees of 
freedom confine the system to a portion of the available phase space. In the finite 
system, on the contrary, all barriers can be surmounted in a finite (but large) time, so 
the degrees of freedom can be frozen only for times smaller than the typical time-scale 
for barrier crossing. The large number and complex structure of states in the phase 
space led Sompolinsky to postulate a set of very long time-scales, eventually diverging 
for A'^ oo, for (free)energy barrier crossing. The times are organized hierarchically, 
i.e., denoting them with tx, where x is an index varying for convenience in [0, 1], then 

lim tx — oo, but — > oo if x > x'. (187) 

A'' — *oo tx 

to account for the ultrametric organization of states. When the system is observed 
at time tx ^ t <^ tx' all degrees of freedom with relaxation times tx shorter then 
t will have relaxed completely, while those with longer relaxation times will remain 
essentially frozen and cannot contribute to the response at time t to an external 
perturbation at time zero. The FDT must then be modified to account for the missing 
contribution of the frozen degrees of freedom. This leads to an anomalous response 
term in the response function which measures the degree of FDT violation. Since 
barriers with time-scale larger than t — tx cannot be crossed the correlation function 
cannot decay to zero but relaxes to 

q{tx) =^^dtx)^M) (188) 
In the thermodynamic limit these partially relaxed states will become stable states 
of typical size q{x) = q(tx) since barriers with time-scale larger than tx cannot 
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be surmounted anymore, while those with time-scale smaller than t^ have been 
already crossed several times. The overlap q{x) is a non-decreasing function of x 
and corresponds in the thermodynamic limit to the Parisi's order parameter function 
of the static calculation. The correlation time-persistent part q{tx) can be written as 
a sum of contributions from unrelaxed degrees of freedom |141| [cfr. H162() for one 
time-scale] : 

x' <x 

Since the system is equilibrated for time-scales smaller that t^ FDT must be satisfied 
on that time-scales. The presence of a time-persistent part in the correlation function 
and the requirement of FDT for time-scales smaller that tx leads to an extra term [cfr. 
I|162|l for one time-scale], called anomalous response term and denoted by A^, in the 
response function |141[I158] . 

Like Parisi's, Sompolinsky's derivation of the self-consistent equations for the 
overlap and the anomalous response term is heuristic but allows for a dynamical 
theory which presents many similarities with the static solution. At difference 
with the usual Parisi's solution, however, the Sompolinsky's solution is expressed 
in terms of two order parameter functions: the overlap q{x) and the integrated 
response function A{x), sum of the anomalous response function terms. This extra 
freedom, called "gauge invariance" , reflects the time-reparamctrization invariance 
of the Sompolinsky's solution: any reparametrization t^ — > h{tx) where h{t) is 
an arbitrary well behaved function preserving the relations (I187|l will lead to an 
acceptable solution. We have already encountered this properties when discussing 
non-equilibrium solutions of the spherical p-spin model in section 16.2.11 This fact 
should not be surprising since it just reflects our missing of knowledge on how the 
barrier-crossing time-scales diverge in the thermodynamic limit. Viceversa we may 
also say that this invariance is intrinsic in any mean-field solution since the details 
on how barriers diverge are irrelevant, the only important point is that they diverge. 
We note indeed that the gauge-invariance of the Sompolinsky's solution has its static 
counterpart in the invariance of the replica solution under replica permutations. The 
Parisi's solution corresponds to the "special" gauge |141| 

A'{x) ^ -Pxq'ix) (190) 

where the prime means derivation, which relates the anomalous response term to the 
derivative of the overlap (or time-persistent correlation) at time-scale tx [cfr. (|175|l ]. 
Equation (|19()|l known with the name of Parisi's gauge is actually an FDT relation. 

As already noted in section 16.2.11 the Sompolinsky's solution, and its variants, 
are equilibrium solutions and cannot account for aging phenomena found in spin-glass 
experiments. To deal with them Cugliandolo and Kurchan have proposed a non- 
equilibrium scheme in which the thermodynamic limit iV ^ oo is taken before any 
large-time limit, including the initial time limit ii — > — cx), to force the system to a non- 
equilibrium state j52| . This procedure drives the system to a non-equilibrium regime 
named weak ergodicity breaking in which TTI is lost and the system displays aging, see 
section l5.2.1l for more details. The TTI is recovered only for small time separations (the 
so called FDT regime) where the dynamics is described by Sompolinsky-like equations 
and the correlation function approaches the plateau C = <7ea with the power law form 
1)180(1 with a temperature dependent exponent a. 

The departure from the plateau, i.e., the aging or non-FDT regime, is 
more complex since the presence of many different time-scales related to different 
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(free)energy barriers must be taken into account, and representations like ()179|l or 
H181fl cannot be adequate. Adapting the Sompolinsky's picture of many time-scales 
to the weak ergodicity breaking regime the non-equilibrium relaxation from time s to 
time t is due to the crossing of (free)energy barriers with time-scales between s and 
t. Then, using the time-reparametrization invariance of mean-field solutions, H179() is 
replaced by the asymptotic form valid for large times [cfr. (|189|l ]: 

ag(t,s)^5]c«(/i,(t)//i.(5)) (191) 

i 

where each contributing term Cag will vary in each separate time sector defined by 
two successive barrier crossing and, as in the case of one time-scale, the functions hi{t) 
could be power law with a time sector dependent exponent hi(t) ^ i^'. Using the fact 
that correlations decrease as times become more and more separated, it is possible to 
show that for large times the following relation must hold 

C(t3,il) = /[C(i3,i2),C(t2,il)], ti<t2<t3 (192) 

where / is an associative function which defines the geometry of the triangles described 
by the trajectory in the phase space [2^1 ■ Next one defines the fixed point qi of / 
filiiQi) — li- The intermediate value of the correlation between two successive 
fixed points defines a time-sector. We note that triangles whose side belongs to 
different time-scales, e.g., C(i3,i2) < Qi but C{t2,ti) > qi, are isosceles with 
C{t3,ti) = min[C(f3, i2), C(t2, ii)]. This defines an ultrametric geometry analogous 
to what is found in equilibrium calculations: 

C(t3,ii) = min[C(i3,t2),C(i2,ti)], (193) 

if at least one C(tz,t2), C{t2,ti) is less then gEA- 

The set of fixed points Qi can be either discrete or continuous. In the latter case 
the correlation (|191|l is the limit case of a continuous sum of infinitely many scaling 
functions C^g. As for the Sompolinsky's approach, the term H191|) in the correlation 
function implies an analogous term in the response function [see (|171|l ] which can be 
related to the correlation function through (|175|l with an X{C) no more constant for 
C < 9ea but which coincides with the function x{q) of the static treatment | 159| . 

We note that for the spherical p-spin model x{q) evaluated from statics is different 
from that evaluated from dynamics [1311 \V6'Jl I52| . It can be shown that a sufficient 
condition for the equality between the static x{q) and X{q) is that the system is 
stochastically stable |16UI I16l| . i.e., the overlap probability distribution Pe{q) of the 
system in the presence of a small perturbation must smoothly converge towards the 
probability distribution of the unperturbed system when e — *■ 0. Moreover the limit 
e ^ must also commute with the limit of large times in the dynamics. If this is the 
case, x{q) and X{q) are then equal. This result holds for mean-field spin glasses with 
continuous RSB such as the SK model, but not for the spherical p-spin model where 
the dynamics is dominated by long lived metastable states. 

6.3. Random manifolds and diffusive models 

The basic ingredient of a glass behavior is the appearance of a multitude of long 
lived states, that prevents exploration of the whole phase space. This situation is not 
restricted to glasses but may be present in several, apparently unrelated, far from 
equilibrium problems. Typical examples have been discussed in the previous sections, 
here we shall comment a bit on random manifolds and diffusion models. 
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A typical situation where glassy behavior shows up is when studying the 
dynamics of an elastic manifold, with or without internal structure, in a random 
quenched environment. This problem appears, for example, in flux lattices in high- 
Tc superconductors |162j . interfaces in random fields ilB3,, charge density waves, 
and surface growth on disordered substrates |164l I165| . The competition between 
elastic stress and disorder produces a state with many characteristics of a glass: slow 
dynamics, nonlinear macroscopic response, aging and so on p^6 , 144 . In the mean- 
field limit a manifold in a random media is described by a field theory with a large 
number of components. In this case it is possible to derive a closed set of dynamical 
equations of the type discussed for spin-glass models. The model of a manifold of 
dimension d embedded into a random medium of dimension N is described by the 
Hamiltonian 

n[cl>{x)] - I d'^x (||V</.(a;)p + V^xlx] + ^(t>{x)^) (194) 

where the TV component field (p = (^ii 02, • • ■ , '/'at) gives the displacement of the 
manifold. The mass term fi constrains the fluctuations of the manifold to a restricted 
volume of the embedding space. The potential term V is a Gaussian variable of zero 
mean and correlations 

^' (195) 



y[<A, x\ V\(i)\ x'\ = -NS{x - x') V 
A common choice for V is 



N 



where 9 is the free-energy fluctuation exponent. The models are divided into "long- 
range" models if 7 < 1 and "short range" models for 7 > 1 since in the first case 
correlations grow with distance, while in the second case decay. 

The study of the static (equilibrium) properties of the d ~ limit |167j . i.e. the 
case of a particle moving in a random potential, reveals that the short-range case is 
solved by a IRSB Ansatz, while in the opposite case of long range the full RSB scheme 
is needed. 

The manifold dynamics is given by the usual Langevin equations: 

where 77 is a Gaussian random process of zero mean and variance 

{T],,{x,t)Tj,{x',t')) ^2T6^,5{x-x')5{t-t'). (198) 

To study the long-time dynamics one introduces two-time quantities, which for the 
simple case of d = are the usual correlation and response functions 



C{t, .) = 1 m) ■ 0(.)), Rit, s) ^ 1 (199) 

where h{s) is a small perturbation applied at time s < t. In addition one also considers 
the mean-squared displacement correlation function 



B{t,s) = l^{m)-4>{-s)f) 

^C{t,t)+C{s,s)-2C{t,s) (200) 
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Figure 8. FD plot in tlie model of directed polymer in random media revealing 
the IRSB or two-timescales character of the dynamics. Prom |168| . 



The analysis of the long-time dynamics for N ~> oo reveals that in both cases 
(long and short range models), two regimes can be defined in the relaxation from an 
initial random configuration: (i) FDT regime for large waiting time t,^ and not too 
large time difference; (ii) non-FDT regime for large t^, and time differences. Under the 
assumption of a weak ergodicity breaking scenario, the FDT in the non-FDT regime 
can be replaced by the generalized form 1)1 76|) with a function X which is different for 
short and long-range models jl44| . For short range models, X{C) is solved with the 
two-tinicscales Ansatz, i.e., it is 1 in the stationary sector while it is a constant smaller 
than 1 in the aging sector. For long-range models many time-scales are needed, and 
X{C) is a non-trivial non-constant function, as found for the SK model. This scenario 
has been extended to the d > case in reference fIBSj . 

Related studies have analyzed FDT violations in polymer models. Yoshino jl68j 
considered the directed polymer model in random media (i.e. the random manifold 
p94|) with c — and N = d = 1) and through numerical simulations confirmed the 
two-timescales character of aging dynamics, see figure |H1 Pittard and Shakhnovich 
| 169| have considered a heteropolymer model with uncorrelated monomer-monomer 
interactions. By analyzing the mode-coupling equations they found a two-timescales 
solution that violates FDT as reported for the random manifold model in the short- 
range case. 

The dynamics of the directed polymer reduces to pure diffusion in the absence 
of disorder. This is the well known random walk which in the continuum limit is 
represented by a stochastic variable x{t) satisfying the (stochastic) differential equation 

j^x{t)^ij[t) (201) 

where 77 is a Gaussian random process of zero mean and correlation 
(77(t)77(s)) =2T5(t-s) 



(202) 
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and T is the bath temperature. A simple calculation gives for the two-time correlation 
and response functions 

C{t, s) = {x{t) x(s)) = 2 T mm{t, s) (203) 

R{t,s)^^-^ = eit-s) (204) 

where h{t) is an external field added to the l.h.s of eq. H201|) . As a consequence we 
see that for any t and s 

X(t, s) = T R{t, s) I \dC{t, s) /ds] = 1/2 (205) 

a constant value different from 1 so that the FDT is violated for all times. Despite 
violation of FDT, the model is simpler than those discussed above, indeed the 
correlation and the IRF show a rather simple form as a function of the waiting time: 

C{t,s) = 2Ts, xit,s)= [ dt'R{t,t') ^ s. (206) 







They both depend on s but not on t (as required by causality). This is a rather extreme 
example, however other less trivial cases, such as spinodal decomposition, scalar fields 
at criticality and so on also do exhibit non trivial non-equilibrium behavior |170j . 
More complicated non- linear diffusion effective models have been shown jl71l I172| to 
display FDT violations compatible with a one-timescale aging scenario with a single 
valued FDR X(i^) < 1 that monotonically converges to 1 like in entropic models 
(Section EH). 



6.4- Trap models 

A successful family of models to describe the glass transition are phenomenological trap 
models. The dynamics in the aging regime can be understood in terms of jumps among 
different phase space components, each jump corresponding to a new rearrangement 
of a cooperative spatially localized region. The dynamics of the system can then 
be viewed as an intermittent motion where some regions remain inactive for a long 
time (and no net heat current is present between the system and its surroundings) 
until an activated jump occurs and thermal heat is released from the system to the 
surroundings, and from there, to the thermal bath. Phenomenological trap models, 
contrarily to mode-coupling theories, are based on the activated nature of glassy 
dynamics. Although trap models have appeared from time to time in the literature (see 

and references therein) the most recent and successful study is due to Bouchaud 
| 4U| who has considered its relevance to describe aging phenomena in glassy systems. 

The trap model corresponds to a set of unstructured energy (or free energy) 
traps that live in a "free energy space" without any explicit reference to real-space 
configurations. It corresponds, in many aspects, to the coarse-grained description 
developed in section 14.21 where activated processes are represented as transitions 
between different phase space components TZ that here could be visualized as traps. 
The number of traps, like the number of components TZ in the coarse grained 
description of the phase space, is exponentially large with the volume of the system. 
The model is defined by a set of traps of different depths E (with E > 0) with a 
density p{E) and a distribution of escape times given by the Arrhcnius expression 
t{E) = Toexp{E/T) where tq is a microscopic time and T is the temperature of the 
bath. Note that, in this last expression, the top-level for all barriers is fixed at zero- 
height. The dynamics of the trap model is then described by the ME lf75|l discussed in 
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section in terms of the probability function P{E,t) that specifies the probabihty 
that the system stays in a trap of energy E at time 

^^f/^ = E ^(^'' - E ^(^' (207) 

The rates Z(E\E') are assumed to be given by H82|) . 

z(i?|£;') = w(i?|£;') p(i?) = 4It! ^(^)^^ = 1 (208) 

nE') Jo 

where we have identified il{T',T) = p{E') and where the bare rate W{E\E') = 
1/t{E') has the dimensions of a frequency. Note that this bare rate only depends on 
the energy of the departure trap but not on the energy of the arrival trap. Other 
versions of the trap model include a dependence on the energy on the arrival trap, see 
for instance |173) . Inserting H208I) into (|207|l we obtain, 

where 

The rates il^Ugt satisfy detailed balance if P°'i{E) oc p{E)t{E), where 
\imt^^P{E,t) = The Bouchaud trap model (BTM) gO] is described 

by the distribution p{E) = (1/Tg) exp(— i?/Tg). The static formulation of such a 
model corresponds to the random-energy model (REM) of Derrida |174[ 1175] . Other 
trap models have considered a Gaussian distribution of energies |39l I176| . However, 
the main interest of the model proposed by Bouchaud is the existence of a critical 
temperature Tg where the distribution P°''(£') ceases to be normalizable. In general, 
for any distribution p{E) the temperature Tq that marks the onset of the non- 
normalizability of P°'^{E) is given by |177| . 

1 . _ Km (211) 

To E^oo E ^ ' 

the BTM corresponding to the case To ~ Tg. Let us note that, in the trap model, 
energies are not extensive but finite. Comparing the BTM with the REM, where 
energies are extensive, we observe that the finite T dynamics in the Bouchaud model 
occurs in a range of finite energies around Eo, the value at which the energy freezes 
in the REM below Tg. The same exponential distribution of states, over a finite 
free-energy interval, is found in the SK model 48 , . 

Dynamics in the trap model has been exhaustively investigated in many works. 
In particular, it offers a rather good explanation of magnetic relaxation phenomena 
observed in spin-glasses gUI 11781 1179| , viscosity anomalies in glasses [1771 11761 1180| 
and, more recently, it has been used as a test model to check whether FDT 
violations are well described by the Ansatz H91|l and whether FD plots are meaningful 
l^m]. Correlation and response functions can be defined in the BTM by assigning 
magnetizations to the different traps as is done to analyze the statics of the REM. The 
effect of the magnetic field on the traps is to modify the escape time by the relation, 
Th{E) = Toexp{{E + mh)/T) = t{E) ex.p{mh/T). The resulting FD plots have been 
analyzed by Sollich and coworkers |^ ESI ■ As remarked in reference |^ , equal 
time correlations can be unbounded so proper FD plots are constructed from the raw 
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Figure 9. Normalized FD plots 11041 for the BTM for a Gaussian distribution of 
trap magnetizations with zero mean and variance A^E) = exp{nE /T) for different 
values of n and T = 0.3 (tq = Tg = 1). For each n times shown are 10^, 10^ which 
are indistinguishable since the limiting FD plot has been reached. Note, however, 
that a temperature factor has been absorbed in x in such a way that the slope is 
-1 at equilibrium (dashed line). From 1541 1551 . 

plots (|l()2(l by normalizing correlations and IRF by the equal times correlation at the 
later time CA,B{t, t) as described in (| 104(1 . The ME (|209(l . modified to include the effect 
of the field, is described by a probability distribution P{E,m,t) — P{E,t) a{m\E) 
where (T{m\E) is the probability that a trap of depth E has magnetization m. This 
probability is assumed to be a Gaussian parametrized by its mean m{E) and variance 
A^(i?). As there is no specific meaning attached to the observable m one can 
think of the two functions rn{E),A'^{E) as describing different class of observables. 
Therefore, observable dependence in the BTM refers to dependence of the FD plots 
on the choice of these functions. The following cases have been considered |^ 
m{E) = e-Kp{nE/2T),A'^{E) = or rn{E) = 0, A'^{E) = exp{nE/T) with n > T - 1 
in both cases. Figure El shows some typical FD plots. 

Three are the main results of this study: 1) The FD plots strongly depend on 
the average and the variance of the Gaussian distribution, therefore the FDR and 
the effective temperature are observable dependent; 2) Most importantly, for a given 
choice of observables the effective temperature (|132|l smoothly changes along time 
within a given time sector. In fact, in the trap model the scaling t/t^ is fulfilled 
in the glassy phase T < Tg but two straight lines (typical of two-timescales glassy 
systems) are not observed in the FD plots; 3) For all observables Xoo = 0, see (I231f) . 
supporting the conjecture that this quantity is indeed universal and may have some 
physical meaning (see section 1^^ . 

Among these results 2) seems particularly interesting. Why FD plots do not 
display the characteristic two-step form of relaxational systems with two-timescales? 
Still for observables with zero mean and finite (but E independent) variance where 
C{t,t) = A^, i.e. for observables that can be considered neutral (see discussion in 
section . the one-step shape of the limiting FD plot is absent. The origin of this 
discrepancy is presently unknown and the finding of a trap model that shows a two 
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straight-lines FD plot remains an interesting open problem |. Let us finish these 
considerations by noting that FD plots, such as the one depicted in FigurelHI are more 
characteristic of systems with full RSB where X{C) is a non-trivial function. Quite 
interestingly, it has been shown 179 that the BTM at the critical temperature T = Tg 
has correlations that do not fulfill the simple t/t^ scaling but a more complicated 
dependence with many time sectors and ultrametricity. However, the resulting FD 
plot at criticality shows only very small deviations from the equilibrium behavior 
X = 1 being very similar to the FD plot observed in the Ising chain |182j (see section 

Just before ending this section, let us remind that these results have been endorsed 
by considering the corresponding driven version of the BTM introduced in [1831 1184| 
in the context of rheology. In this explained in section [7.4.11 the equivalent 

of the waiting time is the inverse of the shearing rate. In the stationary state TTI 
is satisfied but FDT violations persist. For the driven model |185j . as well as for the 
purely relaxational model, the same relationship between correlations and responses 
holds and the two models the non-driven and the driven one) show similar behavior. 

6.5. Models with entropy barriers 

Many of the results described in the previous subsections deal with disordered 
models with a complex thermodynamics. However, many aspects of the violations 
of the equilibrium FDT as well as the existence of an QFDT can be investigated 
in the framework of simpler models that are exactly solvable but still retain the key 
ingredients for the emergence of these new properties. In turn, this can help to identify 
the basic ingredient that any sensible general theory must incorporate. 

The scope if this section is mainly illustrative as it intends to present some of 
the basic ideas of previous section 14.21 applied to very simple examples. We will 
focus our discussion on the oscillator (OSC) and the Backgammon (BG) models. 
A comprehensive account of other results about these models can be found in a 
recent review |15j . Both models have a simple energy landscape and dynamics 
is determined by the presence of entropy barriers. The intuitive meaning of this 
term is the following. In general, relaxation in glassy models proceeds by activated 
jumps over energy barriers that allow the system to escape from a given trap after 
reaching a barrier of height B, the typical time needed in this process being given 
by the Arrhenius law r ~ exp(_B/T). Activated dynamics is strongly temperature 
dependent, and for T = the dynamics is completely arrested, the system remains 
trapped forever and correlations do not decay to zero anymore (ergodicity is broken). 
When the dynamics is dominated by entropy barriers the relaxational mechanism 
is different. The system escapes from a trap through a process which involves a 
timescale r which does not directly depend on the temperature but, for instance, 
on the typical energy E of the trap itself, t{E) which usually is a decreasing 
function of the energy E (see Figure [111)1 . At T = relaxation is not arrested but 
proceeds slower as the energy decreases. This type of dynamics is sometimes called 
marginal dynamics |186| as the system wanders around saddle point configurations, 

I After completion of this work, it has been shown that the influence of the dependence of the 
perturbed rates in a field (upon the observable value taken at the arrival trap) are crucial to get well 
defined limiting FDRs and effective temperatures )181l 

§ Now we assume the standard convention for energies being negative rather than positive as in the 
trap model. 
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{configuration space} 



Figure 10. A typical energy landscape in a glassy model determined by the 
presence of entropy barriers. The effective barrier log(T(E)) increases as E 
decreases. From I187| . 



hence is never arrested. Of course, the temperature dependence of the relaxation 
time in entropic models can appears as effectively activated: at finite (but low) 
temperatures the maximum value of t{E) corresponds to the lowest energy reached, 
i.e. T{Eeq{T)), Eeq{T) denoting the equilibrium energy. The temperature dependence 
of the relaxation time in equilibrium T{Ecq{T)) is activated in most entropy barrier 
models. A phenomenological description of these entropy models has been introduced 
by Barrat and Mezard |173j who have generalized the BTM to the case where the 
distribution of trapping times is itself a time-dependent function. The oscillator 
and Backgammon models described in this section are microscopic versions of this 
entropic trap model. According to the scenario presented in section 1?!^ the role of 
entropy appears to be important as relaxation in many glassy systems is accompanied 
by the emergence of a non-equilibrium microcanonical ensemble which determines 
fluctuations and responses in the aging state, leading to the existence QFDT and 
an effective temperature. This fact suggests that a deep comprehension of the glassy 
dynamics in exactly solvable entropy barrier models can be a first step toward grasping 
the leading aspects behind the behavior of more realistic systems, where both entropy 
and energy barriers simultaneously intervene. 

6.5.1. Oscillator models We begin our discussion by describing the OSC model in its 
simplest version. A review of some results can be found in _68j 188. .15) . Originally, 
oscillator models where introduced indirectly in the analysis of the Monte Carlo 
dynamics of the spherical Sherrington-Kirkpatrick model, which can be mapped to 
a set of disordered harmonic oscillators ^189. ,190) . The OSC model is obtained by 
simplifying the previous one to an ensemble of identical harmonic oscillators |191| . 
The OSC model is defined by the following energy function 




(212) 
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where the Xi are real-valued displacement variables for the N oscillators and K > 
is the Hooke constant. The equilibrium properties are trivial due to the absence of 
interactions, but the Monte Carlo dynamics couples the oscillators in a nontrivial way. 
Moves are proposed according to the following rule, 

(213) 

N 



where the are Gaussian random variables with zero mean and variance A^. Moves 
are accepted according to the usual Metropolis rule. Each move is a parallel update of 
the whole set of oscillators. Both the energy function (|212|l and the dynamics H213() are 
invariant under rotations in the iV-dimensional space of the Xi . This symmetry makes 
the dynamics exactly solvable, and many non-equilibrium quantities e.g., effective 
temperatures and FDT violations can be tackled analytically. 

The OSC model is a classical model where the equilibrium entropy is given 
by S{T) — ilog(T), thereby diverging when T — > as expected for a continuous 
model (the ground state is a zero-measure state corresponding to the absolute global 
minimum of (|212|l . i.e. the configuration Xi = 0, Vj). At T = only those moves 
that decrease the energy are accepted, therefore as the system approaches the global 
minimum the frequency of accepted moves H213|l dramatically decreases. However, 
that frequency never vanishes so dynamics is never arrested. Dynamics slows down 
because phase space directions where energy decreases are exceedingly difficult to find. 
The full solution of the OSC model has been presented in 

The main physical quantity containing detailed information about the dynamical 
evolution is the probability density of energy changes P{AE). This quantity expresses 
the probability density that a proposed move (|213|) changes the total energy of the 
ensemble by the amount AE. P{AE) was originally derived in |191| using standard 
integration tools. Here we present two other alternative derivations which help to 
understand the mechanisms behind slow relaxation. 

The first method relies on the Gaussian form of the distribution while the second 
one uses a microcanonical argument to count the number of accessible configurations 
from a reference configuration with a given energy E. The first derivation is rather 
simple as the distribution for AE can be easily obtained. Indeed using (|212ll and H213|l 

^E^^Y.'.'-' + ^JZrl (214) 

^ i i 

and the Gaussian property of r^, it follows that AE has a Gaussian distribution. 
The mean and variance of the distribution are Mab = AE = KA'^/2, (j^e = 
[AEY - (AEf = 2EKA^/N, yielding [TM] . 

(AF— KAl)"^ 

F(A£;) = (4^eXA2)"*exp ^ ^ ) 



AeKA^ 



(215) 



where e = E/N is the energy per oscillator. 

The second method to derive (|215|l is based on a microcanonical computation. 
In Figure 1111 we depict a schematic two-dimensional representation of the motion 
of a representative configuration in phase space. The reference configuration {x^} 
at a given time has an energy E and lies on the spherical hypersurface of radius 
R = \J2EIK (depicted as the point P in the figure) with center at the origin 
{xi = 0} (point O in the figure). The smaller dashed circle represents the region 
of points accessible from according to the dynamics (|213|l . All accessible points 
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Figure 11. Geometrical construction to compute P{AE). The thick hnes 
denote the departing and final energy hypersurfaces centered at O. The dashed 
circle indicates the hypersurface accessible from point P. The intersection region 
between the accessible hypersphere centered at P and the final hypersurface of 
energy E' defines a hypersurface I of radius C (the radius is represented by a 
thick line). See the text for explanation. 



{xi} satisfy J^ii^i ~ ^i)^ — i-^- lie at a distance A from {x^} which is the radius 
of the smaller dashed circle. The accessible configurations in a single move lie in a 
spherical hypersurface of dimension N ~ 2 corresponding to the intersection of the 
hypersurface of energy E' and the smaller spherical hypersurface of radius A. We call 
this region the intersecting region / as shown in Figure ITTI The final configurations 
contained in / lie at a distance R' = \J2E' jK to the origin O. The change in energy 
associated to this transition is Ai? = E' — E. The probability of this jump is therefore 
proportional to the surface of the intersection region, P{AE) cx C^^^ where C is the 
radius of the intersecting region. The computation of C is quite straightforward as 
can be deduced from the triangle including the points P,0 as vertices and whose three 
sides are R,R',A. In terms of R, R' and A, the distance C is given by the relation, 

2 



A^- 



K f2AE 



8E \ K 



(216) 



The surface ^{E, AE) corresponding to the region / of radius C, relative to the energy 
E of the reference configuration a;" is. 



n{E,AE) oc C 



N-2 



A'- 



K /2AE 
8E [ K 



Using the fact that E is extensive with N this expression can be rewritten as. 



n{E, AE) (X exp 



" A{E/N)KA'^ 



(217) 



(218) 



which is proportional to the probability distribution H215II . 

The distribution P{AE) is depicted in figure IT^ for different values of the energy. 
As E decreases the number of moves with AE < shrinks as the total area under the 
curve with AE < decreases. From (|215|l the dynamical equations immediately follow 
by defining a Monte Carlo step as a collection of M elementary moves (each elementary 
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Figure 12. Probability density of the energy of the proposed moves for different 
values of the energy E as defined in Figure ITT1 



move corresponds to a global change of all oscillator coordinates as described in H213|) ). 
Because the average change of energy AE is finite in an elementary move and the 
total energy (|212|) is extensive, the number of moves M in a Monte Carlo step must 
be proportional to N. For simplicity we will take M — N. In the limit N ^ oo, time 
becomes a continuous variable. This allows to write a closed equation for the energy 
E and acceptance rate A, i.e., the average number of accepted moves in a Monte Carlo 
step, 

arp poo 

— =/ xP{x)W{x)dx (219) 

^t J —oo 

BA 

— = y P{x)W{x)dx (220) 

where P{x) is the probability distribution H215|) and W{x) is the transition rate which 
ensures that detailed balance is obeyed. For instance, according to the Metropolis 
rule W{x) ~ min[l, exp(— /3a:)] or in heat-bath W{x) = + exp(/3a:)]. At zero 
temperature the transition rate rules W{x) converge to W{x) — 0{x). In this limit 
1^1^ . if^^ become 

dE 

xP{x)dx (221) 



dt 



A{t) 



P{x)dx 



(222) 



As shown in Figure El the acceptance is given by the shaded area enclosed in the 
negative tail of the distribution. At zero temperature, according to H215(l . both the 
mean and the width of the Gaussian decrease as well as the shaded area in Figure [T^ 
implying a systematic decrease of the acceptance rate. 

We do not want to dwell here on all results one can learn by solving the dynamical 
equations (UnD , I222l (see [roTHTUg) ). Interestingly, in the OSC model the FDR ifnTjl 

can be exactly computed at any temperature | 
1 dE{s) 



for the magnetization M = 

reff(s) = 2E{s) - 



f{s) ds 



(223) 
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and the QFDT is satisfied. At T — in the large time hmit one gets T^sis) — 2E{s) 
plus subleading corrections. This result shows that FD plots are straight-lines starting 
at C(s,s) = {2E{s)/K),xis,s) = and finishing at C = 0,x = l/K. The relation 
between the effective temperature and the dynamical energy is exactly the same as 
the equilibrium relation given by the equipartition theorem. The aging system is in a 
quasistationary state where relations between dynamical quantities are formally the 
same as in equilibrium. This allows to define a time-dependent configurational entropy 
Sc{E) through the relation (p^ . 

^ • S,{t)^S,{E{t))^lHE{t)) . (224) 



T^ff{E) dE ' ' y " 2 

Until now we have discussed some of the analytical results obtained by solving 
the dynamical equations of the OSC model. However, an interesting question is the 
following: can we determine the value of the effective temperature from the sole 
evolution of the energy E without having to analyze correlations and responses in 
the framework of the QFDT as described in (|223|l ? Ideally we would like to apply the 
ideas presented in section ll^^ to identify the value of T^^ solely from the off-equilibrium 
transition rates Z{F\F'). In that description dynamics proceeds by activated jumps 
over different states, whose dynamics is described by the free-energy ME (|79|l . What 
are the states in the reduced description of the OSC model? As the energy landscape 
is a single parabolic well it appears that a reduced description is not possible. The 
clue to this question is easy to understand if one realizes that at zero temperature the 
acceptance of the dynamics goes to zero with time, therefore each time a proposed 
elementary move is accepted we can effectively consider that the system has jumped 
from one state to another, the typical time for this jump steadily growing with the 
time elapsed since the system was quenched. In this view, each state corresponds to 
a configuration and the reduced dynamics simplifies. Moreover, the free energy of 
the state is simply the energy of the corresponding configuration. In the asymptotic 
regime {dE / dt) / E << 1 where finite size effects are not important, i.e. A <C y/2E/K, 
the probability distribution describing the energy change after the first jump is given 

by, 

P{AE) - exp(^)6i(-A£;) (225) 

where we used (I215|l and expanded it around AE = up to the linear order. Using 
relation (|224|l we can recast (I225|) in the following form, 

P{AE) = -— ^7— exp ( ) 0{-AE) (226) 



2T,s{E) '^\2T,s{E), 

showing that the statistics of energy jumps is an exponential with a width that directly 
depends on the effective temperature. This result has two consequences: 1) It shows 
that the OSC model is a microscopic version of the trap model proposed by Barrat and 
Mezard |178| : 2) The effective temperature could be computed from the statistics of 
the first free-energy jumps among components (here corresponding to configurations). 

Before finishing, let us note that a number of variants of the oscillator model 
have been considered, all sharing the feature that oscillators do not interact. For 
example, Nieuwenhuizen and Leuzzi |108[ll92Tll93[[TM| have considered a model with 
spin variables in addition to oscillators. The new variables are discrete and used to 
mimic fast relaxational processes not contained in the original formulation; aging, slow 
relaxation in these models, can be still described in terms of entropy barriers. 
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6.5.2. The Backgammon and urn models Another instructive model where relaxation 
is determined by entropic barriers is the Backgammon (BG) model |195j . The model 
belong to a large class of models under the name of urn models where N particles or 
balls are distributed among M urns or boxes. The BG model is defined by the energy 
function, 



where stand for the occupancies for each box and A^empty stands for the number of 
empty boxes. The model has different versions |196| according to whether particles are 
distinguishable (Maxwell-Boltzmann statistics) or not (Bose- Einstein statistics) . The 
easiest way to compute the thermodynamic properties of the BG model is by expressing 
the partition function in terms of the occupancies Ur rather than in terms of the set 
of boxes occupied by the particles. In the Maxwell-Boltzmann case, thermodynamics 
needs to be corrected by dividing the partition function by the usual Nl term to avoid 
Gibbs' paradox. In the model with Bose-Einstein statistics this is not necessary. The 
dynamical rules of the BG model directly depend on the type of statistics considered. 
In the Maxwell-Boltzmann case, a departure box d is selected with a probability 
proportional to the occupancy Ud of that box and a new arrival box Ua is selected 
with uniform probability. The proposed move is given hy Ud ^ l,na Ua + I 

and accepted according to the standard Metropolis rule. Instead, in the Bose-Einstein 
case, the proposed move and the transition rate are the same as for the Maxwell- 
Boltzmann case but both the departure and the arrival box d and a are selected with 
a uniform probability 1/M. 

The resulting dynamics of the model is strongly dependent on the type of statistics 
considered, the interesting one corresponding to the Maxwell-Boltzmann case which 
displays a strong separation of time-scales. In what follows, although otherwise stated, 
we will concentrate in the analysis of the model with M = N. 

The thermodynamics of the model is exactly solvable and, the entropy being given 
by S(e) — log(l -I- e) where e = E/N denotes the entropy per box (or per particle). 
To analyze the dynamics from the perspective of the scenario described in section [4.21 
we follow a similar reasoning as we did for the OSC model. Let us consider a system 
quenched down to T = and the aging regime reached in the asymptotic long-time 
regime where (dE/dt) / E << 1. In that limit, the system has a number of empty boxes 
iVonipty — —E and further decrease of that number by one unit Ai? = — 1 requires a 
time that exceedingly grows as E decreases toward its minimum value E = —N -\- 1 
(all particles condensed into a single box) [1971 1198j . Therefore, as relaxation slows 
down, for a long time the system wanders through the constant energy surface by 
exchanging particles among occupied boxes. When a new box is emptied the energy 
decreases by one unit. As we did in the OSC model, also in the T = dynamics of the 
BG model each component TZ corresponds to a single configuration with free energy 
equal to the energy of that configuration. 

The transition rate for going from a configuration of energy £^ to a configuration 
of energy £^ — 1 is sole function of the number of available configurations in the initial 
and final state. 



M 




(227) 



r=l 



Zt{E-l\E) _ n{E - 1) 
Zt{E\E) " n{E) 



(228) 
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where the rate frequency l/t{E) associated to (for its definition, see Section 
has canceled out from the numerator and denominator in the left hand side of (|228|) . 
As E is extensive with N, and using H83|l we obtain, 



Equation H229|l provides a way to estimate the effective temperature by looking at the 
number of moves required for a move to decrease the energy by one unit. Inverting 
yields 



Therefore the ratio Zt{E\E) /Zt{E — l\E) is simply given by the inverse of the number 
of accepted moves N* required to decrease the energy by one unit. We stress that N* is 
not the number of MC steps but rather the number of accepted elementary moves with 
AE = required until a move with AE = — 1 is found. This number iV* is independent 
of the size of the system. Note also that the rate of rejected moves (which gives the 
acceptance) does not enter into the expression of the effective temperature but rather 
into the value of the typical timescale t{E) associated to the jump. This is an essential 
ingredient required in any "quasi-equilibrium" or microcanonical description of the 
relaxation. The value of the effective temperature at a given time only depends on the 
number of accessible states with lower free energy (energy in the present model) and 
not on the time necessary to escape from that state. This expression is very amenable 
to numerical calculations. 

Note that the energy levels in the BG model are discrete, therefore the relation 
H226|l cannot be applied. However, one can consider a non-degenerate disordered BG 
model where boxes are assigned different energies |187) . In this case the distribution 
P{AE) has been derived and shown to have a exponential time-dependent behavior 
characteristic of trap models. However, the computation of the effective temperature 
in that model, and the verification of (|226|l has not yet been done. 

The Backgammon model as well as some of its variants has been solved by the 
technique of the generating function 11101 EOOI IMl EOS EHl EOl EOSI- A 

discussion of some of the main results and analytical techniques has been presented in 
|15|. The origin of the effective temperature as derived from the QFDT remains yet 
not completely understood. No neutral observable has been yet identified in the BG 
model. Although all studied observables show that the QFDT H94|l is verified in a one- 
timescale scenario, there seem to be different effective temperatures depending upon 
the observable considered [2021 121)41 l2U6j . On the other hand, none of the different 
effective temperatures associated to these observables appears to be linked to the 
effective temperature obtained within the adiabatic approximation |199) . An interest 
variant of urn models are zeta urn models introduced in the context of quantum gravity 
that show a finite-temperature Bose- Einstein condensation transition |207i i208| . see 
section 

6.6. Ferromagnetic models at criticality 

A new emerging line of research related to the mainstream of researchs on FDT 
violations in glassy systems, is the study of FDT violations in critical points. In critical 
points the relaxation time diverges and one can investigate, for instance through field 




(229) 




(230) 
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theoretical methods, the time-dependence of correlations, responses and the resulting 
FDR H91|) in the asymptotic regime s — > oo. These type of investigations are not 
different of those undertaken in the study of coarsening behavior in ordered phases. 
The main difference between the slow dynamics in a critical point and coarsening is 
that, in the former case, critical slowing down is at the origin of the slow dynamics. 
Interfaces have no stiffness tension and their motion is subdiffusive (or diffusive) 
and only consequence of the curvature of the interface. Growing domains are not 
islands of up or down spins but regions of spatially and temporally correlated spins 
of zero net magnetization. At Tc dynamics is described by the renormalization-group 
dynamical equations of the corresponding finite-temperature fixed point. However, 
in coarsening systems below Tc activated processes are important and interfaces have 
non-zero tension as they separate domains of up and down spins. Competition between 
the curvature of the surface and its tension leads to different growth laws, described 
by the zero-temperature fixed point. Because the origin of slow dynamics at Tc is 
different than in standard glassy systems, the main Ansatz (|91|l may not hold. An 
indeed, many studies of the FDR at criticality reveal that X{t,s) is not a single 
function of the correlation C{t, s). 

One among the first studies of ferromagnetic models at criticality is the Ising 
chain solved by Glauber in 1963 |2()9| . Strictly speaking the dynamics of this model 
is that of a coarsening model as the system orders at the critical point T = where 
the magnetization discontinuously jumps from (for T > 0) to 1 (at T = 0), i.e. 
the critical exponent /3 vanishes. We will see below that the ferromagnetic Ising 
chain has some peculiar properties. Starting from a random initial configuration 
the coarsening dynamics at T = has long been studied in j21()| and revisited in 
| 211) . In the asymptotic long-time limit, the aging part of the two-times correlation 
functions scales like Cag(t, s) ~ F{L{t)/L{s)) (|252|l with L{t) ^ t^^^ corresponding to 
a diffusive process of the interfaces (see section 1731 . More recently, and nearly at the 
same time, the FDR has been analytically computed |212II182]| for a random staggered 
perturbation finding X{C)^l/[2- sm'^{TTC/2)] in agreement with the Ansatz H91(l . 
Or, in terms of times, X{t,s) = (1/2)(1 -I- s/t) showing that X — > 1/2 if t ^ oo 
or C 0. This last result coincides with that found in the random walk or in the 
Gaussian model |170| . all models characterized by a diffusive dynamics. This has led 
to the proposal |2()51 1213] that, in systems at criticality, the limiting value X{t, s) for 
t ^ oo is a universal quantity, 

Xoo = lim lim X{t,s). (231) 

s — 'OC t — >oc 

This conclusion has been substantiated by the exact solution of the ferromagnetic 
spherical model in general d dimensions |2()5I I213j . The authors have noticed how 
is only a function of amplitude ratios describing the leading scaling behavior 
of correlations and responses. This conclusion is endorsed by the following result: 
Xoo — 1/2 for d > ducd — 4 and X^o = 1 — 2/d for 2 < c? < ducd — 4 where ducd stands 
for the upper critical dimension (below d = 2 there is no finite T transition in the 
model). Interestingly, for d = 2 the FDTR X{C) has a non-trivial form |214| . similar 
to what is found for the ferromagnetic Ising chain. Numerical results in the Ising 
model show that X^o — 0.26 in 2d || and preliminary results in 3d give X^c — 0.4. 
In figure El we show the the IRF corresponding to the TRM susceptibility H98() as 
function of C in the 2d Ising model at the critical point. The general scenario about 

II After completion of this review we learned from recent results finding a slightly different value 
Xoo ^ 0.34 in 2d ET51 . More work is necessary to accurately estimate that number. 
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Figure 13. IRF 19^1 corresponding to tlie thermoremanent magnetisation in the 
2d Ising model at criticality. The full and broken line correspond to the quasi 
stationary regime X = 1 and Xao = 0.26 respectively. From 12051 . 

ferromagnetic models at criticality is as follows |213j . Let us consider a system that is 
quenched from a random configuration to a temperature below Tc- If t = < — .s <C .s the 
system is in equilibrium, TTI holds and FDT is not violated {X{t, s) = 1). However, 
if t/s ~ the system ages, and both TTI and FDT are violated. A solution for 

the correlation Cag(i, s) and response i?ag(i, s) that matches the intermediate regime 
between the stationary X — I and the aging regimes is given by, 

ag(i, s) = ml^ C [^-^ , i?ag(t, s) = R (232) 

where a > is a coarsening exponent. From (|232ll X{t, s) ~ s^^/m^q, therefore X ^ 
for coarsening systems where rucq is finite. The same expression is valid for the critical 
point but replacing mcq by its time-dependence at Tc. Using mcq ~ (Tc — T)'',^ ~ 
(Tc — T)~'^, t ~ where /3, v, Zc are the, correlation length and dynamical exponents 
respectively. Substituting these relations into H232(l and using ac = 2(3/ ZcV gives at 
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leading to X{t,s) — X{t/s). This result has two consequences: 1) for t oo, 
X{t, s) and 2) only for Qc — 0, according to the left expression in H233I) . X{t, s) 

can be expressed as sole function of Cag. This is the case of the aforementioned Ising 
chain where the Ansatz is valid because (3 = 0. It has been suggested |216| that 
Rag (t, s) covariantly scales under conformal transformations of time leading to specific 
predictions for the scaling function R. An important implication of the result Ac = 
is that the aging part of the IRF (|97|) does not vanish in the asymptotic long-time 
limit. As noted in '217' the non- vanishing of the aging part of IRF in the large s 
limit is related to the failure of the scheme that links static and dynamical properties 
11601 ll()lj in the Ising chain. In fact, as discussed in section 16.2.21 the stochastic 
stability property links the equilibrium P{q) with the behavior of IRF, 



where we have used ()100ll . In the Ising chain P{q) = 6{q — 1). Equation H234|l 
is not fulfilled in the Ising chain as can be easily checked inserting the result 
X{C) = 1/(2 — sin^(7rC/2)) (discussed in the paragraph preceding H231(l ^. As soon 
as Tc is finite {d > 1), X(C) = 9{C - m^^) for T < and (^Ijl is again satisfied. 
Despite its simplicity, the Ising chain appears to be an interesting solvable example 
that allows to check many results. For instance, observable independence has been 
also recently addressed [SHI showing that X^o = for a large class of observables. 
Recent progress in the study of the FDR at criticality has been achieved by Calabrese 
and Gambassi |218[ 1219) who have considered the FDR in momentum space X^{t, s). 
The study of 0{N) models using field-theoretical techniques yields estimates for the 
value of Xoo H231() in an e = 4 — d expansion. Two-loops computations |2 1 9j give, for 
the Ising case = 1, the following values: X^{3d) = 0.429(6), Xoo(2d) = 0.30(5) 
compatible with the results obtained from numerical simulations 205, 215 . 

Other studies of models at criticality include the XY model with a Kosterlitz- 
Thouless transition [220 |. This model has a low-temperature phase where correlation 
functions decay algebraically, therefore correlations are critical below Tc. A failure of 
the stochastic stability property, similar to that reported in the Ising chain, appears 
in the 2d XY model at low temperatures where the density of vortices increases. 

Finally, let us comment on results for zeta-urn models that have confirmed the 
validity of the relation H231fl |207l 12081 ITS] . Zeta urn models show a quite rich phase 
diagram described in terms of the density p of balls (i.e. the number of balls divided 
by number of boxes) and the temperature. There is a critical line pdP) that separates 
a fluid regime (p < pdP)) from a condensed regime (p > Pc(P)) with glassy dynamics. 
Along the so called regular part of the critical line (/3 > 3) X^o = is temperature 
dependent. This number lies in the interval [4/5, 1] quite far from the values X^ < 1/2 
found in ferromagnets. In the condensed phase X^o T^l"^ at low temperatures. It 
vanishes at T = as in coarsening systems, however the condensation dynamics in 
urn models is totally inhomogeneous in contrast to the homogeneous character of 
coarsening in ferromagnets. 



P{q) = -T 



dC^ 
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7. QFDT: the numerical evidence 

Computer simulations have the great advantage, over the real experiment, of direct 
access to the microscopic level, even if only relatively small timescales and system 
sizes can be studied. Although this can be a serious limitations, the true fact is that 
the numerical study of aging phenomena and FDT violations has been successfully 
done during the last years for several systems. 

7.1. Structural glasses 

Using ideas developed in the field of spin glasses, many conjectures have been 
formulated concerning the structure of the phase space of glassy systems. However, 
obtaining direct information, either from experiments or from numerical simulations, is 
a difficult challenge. Relaxation times in a glass are so long as to preclude equilibration 
within experimental times. Numerical or experimental exploration of the phase space 
in these systems is therefore necessarily incomplete. The increase in computational 
power and the recent developments in the theory of disordered systems has pushed 
forward an approach which should not suffer from these limitations. The idea, which 
was actively developed in the study of spin glasses, is that relevant information on the 
phase space structure should be hard-encoded into the non-equilibrium dynamics of 
glassy systems. 

According to the conjecture of the similarity between structural glasses and some 
spin-glass models, X{C) for structural glasses is a two valued function with X{C) = 1 
at short times, and X{C) = rn < 1 in the long-time aging regime. This scenario has 
been largely studied using numerical simulations. 

In a numerical investigation of aging effects not only the waiting time must be 
changed over several order of magnitude, but for a given waiting time the subsequent 
dynamics must be studied over a long time. For these reasons aging phenomena have 
been studied for models that are simple enough to be simulated over long times, but at 
the same time still catch the essential features of real glasses. Moreover, to maintain 
the systems in a non-equilibrium state for very long times, crystallization must be 
strongly inhibited. This is obtained either with a particular choice of interaction 
potential parameters or by adding a (small) extra term in the potential. In the 
following we discuss results obtained for several models. 

7.1.1. Mixtures of soft particles of different sizes 1^ 1^ 1^ 1^ 1^ . 

This system consists of N particles half of which are of type A and half of type 
B interacting via the Hamiltonian 



where the radius depends on the type of particle. It is known that the choice of 
two different types of radius such that r^/r^ — 1.2 prevents crystalhzation and the 
system can be brought into a glassy state. 

Due to the simple scaling of the potential, the thermodynamic quantities depend 
only on r = p(3*, where p is the density which can be taken equal to one. This model 
presents a glass transition at about Tc ~ 1.45 |222| . In figurelHIwe report the response 
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where fi is a random Gaussian vector of squared length equal to the space dimension 
d, versus the self-diffusion function 



A{t^, + T, tu 



1 



N 

TV ^ 

4=1 



{\x^{t^ +r) - 



(237) 



The average is over different initial states at t^ and realization of /. Two linear 
regions with different slopes, one with X{C) = 1 and one with X{C) = m < 1 are 
clearly seen, in agreement with the two-timescales scenario. The dependence of m 
with T is well fitted by the spin-glass model prediction m{T) = T/Tc, see figure [131 
Similar results have been obtained in monoatomic Lennard- Jones glasses [22 7. i228] , 
see section (17.1.31). 



7.1.2. Lennard- Jones binary-mixtures [^l^l^l^ 

The system consist of a mixture of particles of type A and B of equal mass rn 
interacting via a 12 — 6 Lennard- Jones potential of the form 



(238) 



where Cafi and (7af3 depend on the particle pair type and are chosen to prevent 
crystallization. For a 80 : 20 mixture, and using eAA and ctaa as units of energy 
and length and (mCT^^/48eAA)^^^ as the unit of time, these are eAA — 1, ctaa — 1, 
eAB = 1.5, aAB = 0.8, cbb = 0.5 and aBB ~ 0.88. The atomic dynamics of this 
model is well described by the mode-coupling theory with a critical temperature of 
Tc = 0.435 in reduced units. 
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Figure 15. The quantity m = -g^^ as = 2048 as function of tlie temperature. 
The straight line is the prediction of the approximation m(T) = T/Tc. From 

Typical FD plots numerically obtained by Barrat and Kob |233l 12341 1235j are 
shown in figure [TCI The correlation is given by the incoherent scattering function for 
a wave vector k 

Ckit^+T,t^) = 1^ e^fe-[-.(*™+-)-'-.(*™)l (239) 
j 

while the response is measured by adding to the potential a term of the form 

SV ^VqY^ tj cos(fc • rj) (240) 

where Cj — ±1 with equal probability and Vb < T. Again a two-timescales scenario is 
clearly seen. Moreover the effective temperature in the slow regime where m < 1 is in 
reasonable agreement with the glass transition temperature Tg of the system. f 

Thus the FDT is broken as the system fails to equilibrate, as expected from 
spin-glass models. 

7.1.3. Monoatomic Lennard-Jones systems PSII^I^I^ . 

The system consists of equal particles interacting via the potential V — Vlj + SV, 
where Vlj is the usual 12 — 6 Lennard-Jones potential (expressed in reduced units), 
and SV is a many-body term that inhibits crystallization: 

= f E ^(^(9) - ^o) (S{q) - Sof (241) 

where S{q) is the static structure function, a = 0.8 and Sq = I. The sum is made 
over all q with q^aK — ^ < \q\ < Qmax + A, where qmax ~ 7.12 p^^^ and A — 0.34, p 
being the particle density. 

t The glass transition Tg is defined as the temperature below which the system fails to equilibrate 
on the experimental time-scale. For structural glasses Tg is defined as the temperature at which the 
viscosity is equal to lO'^^ Poise or, equivalently, a relaxation time of 100 s. 
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Figure 16. Parametric plot of the integrated response function M{tw +r, t^,) 
and the correlation function C'{tm +t, tm) for k = 7.25. Final quench temperature 
Tf = 0.3, t„ = 10 00. C ircles: = 10000. The straight lines have slopes -1.0 
and -0..45. From l234l . 



In figure 1171 we report the parametric plot of the mean square displacement A 
and IRF in a crunch experiment. It is important to note that the temperature below 
which m < 1 does coincide with the glass transition temperature of the system at 
the density reached after the crunch. This shows that the breaking of FDT does not 
depend on the initial state nor on the path followed in the (T, p) plane, but only on 
the final (non-equilibrium) state to which the system is brought. 

The FD plots for glass-forming liquids discussed here reveal the typical two- 
timescales (or IRSB) scenario found in some spin-glass models. This supports the 
original Goldstein's idea that the phase space of supercooled liquids is divided 
by high barriers into different valleys each with its own statistical properties. This 
picture has been recently confirmed by a a direct analysis of the motion of a glass 
forming liquids in terms of IS 62 . 

In the right panel of figure^] 2^S| we show the temporal behavior of the average 
energy of minima for a binary mixture Lennard-Jones system visited in the non- 
equilibrium motion following a quench to a low temperature. As a consequence 
when the system is quenched from an high-temperature state to temperature T, the 
fast intra component degrees of freedom will quickly equilibrate with the thermal 
bath temperature T. Applying the condition of minimum free energy to the 
system, constrained to stay in components of depth Eis, allows to define an effective 
temperature 

id/dE,s)F4T,E,s) 

which reflects the non-equilibrium net heat flow from the system to the thermal bath 
p.00 . This expression coincides with that proposed in in the context of p-spin 
models, once the components are identified with the Thouless-Andreson-Palmer states 



Violation of FDT in glassy systems 



75 



y2 

0.3 



* 0.2 
H 



0.1 
0.0 



1 1 1 1 1 


' \ ' \ y 


. O 2R/p . 




O - 

o 




- - 




-6 




■ 





CO. 




Figure 17. Mean square displacement A and integrated response function 2i?//3 
in reduced units at p|=1.24 for three temperatures T*=0. 48, 0.96, 1.76. The left 
side shows the log time dependence of the two quantities. In the right side the 
response function vs A. Dashed lines indicate equilibrium FDT, while full lines 
fit the off-equilibrium aging region. From 12271 . 



| 71l [7S| , see also the discussion after (|90|l . Inserting into H242|) the value of £"13 as 
function of time one finally gets the function Tcs{t). The definition of Tcff is shown 
graphically in figure El |238j . 

The two-timescales scenario is rather well confirmed by numerical results, as 
shown in the FD plots of figure El [2'^H\ , where the FD plot for the binary mixture 
Lennard- Jones system is reported. The full lines are the prediction from H242|l . the 
agreement is rather good. 

7.1.4- Finite-size mean-field glasses |65l Ititil 1^ 

We have seen in the previous sections that the essential features of MCT for 
glass-forming systems are also common to some fully connected spin-glass models, 
called mean-field p-spin glasses with p > 2. In the thermodynamic limit, the high- 
temperature paramagnetic phase is described by the schematic mode MCT for super- 
cooled liquids. At the critical temperature Tc an ergodic to non-ergodic transition 
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Figure 18. Left: Solutions of 12421 for several values of of the final quench 
temperature Tf for the BMLJ system. Right: eig as a function of time, following 
the temperature quench. The arrows show graphically the procedure which 
connects the eig(t) value to Tgff value, once Tf is known. [Data courtesy of 
F. Sciortino and P. Tartaglia, see also Ref. 1^^ ]. Prom 




Figure 19. Response (p^(t)} versus the dynamical structure factor S^"{t) = 
(Pk (*)Pk* where is the Fourier transform component of the density of a = 
A, B particles at wave-vector k, for the binary mixture Lennard-Jones particles 
system for a quench from initial temperature Ti = 0.8 to final temperature 
Tf = 0.25 and two waiting times = 1024 (square) and = 16384 (circle). 
Dashed lines have slope T^^ while thick lines have slope . [Data courtesy of 
P. Sciortino and P. Tartaglia, see also Ref. ESHl]- From 
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takes place. In mean-field models the relaxation time diverges at T^. as barriers 
separating different ergodic components become proportional to the system size, 
thereby diverging in the thermodynamic limit. In real systems, or glass models just 
described, the barriers are of finite height and the transition to a glassy state appears 
at the glass transition temperature Tg < Tc, where the typical activation time over 
barriers is of the same order as the observation time. 

Despite these differences mean-field models, having the clear advantage of being 
analytically tractable, are a very useful tool to study the phase space structure of glassy 
systems, especially between the dynamical temperature Tc and the static temperature 
Trsb (Kauzmann temperature Tk in glass language). The main drawback is that, since 
activated processes are not captured by mean-field models, the picture that emerges is 
not complete. To go beyond mean-field it is necessary to include activated processes, 
a very difhcult task since it implies the knowledge of the excitations involved in the 
dynamics. A simple approach is to include finite-size effects in the dynamics of an 
infinite mean-field system just extending the analysis to finite-size mean-field. This 
approach has been suggested by Nieuwenhuizen |49| and is somehow reminiscent of the 
dynamical approach of Sompolinsky jl41j , see also section 16.2.21 We stress that the 
assumption that finite-size mean-field models capture the physics of glasses beyond 
MCT is not trivial. In fact, activated process in finite-size mean-field models could 
be different from those of supercooled liquids, making the behavior different. This, 
for example, seems to be the case of the Potts Glass model, where recent studies on a 
finite-size version indicate some differences with the fragile-glass scenario p39, 240 . 

A spin-glass model in the p-spin universality class, that displays a fragile-glass 
behavior, is the Random Orthogonal Model (ROM) piTUT^ . The model is defined 
by the Hamiltonianf 



where Ui are N Ising spin variables {a — ±1) and Jij is a random N x N symmetric 
orthogonal matrix with zero diagonal elements. We note that at difference with 
previously discussed spin-models, the condition of orthogonality leads to a strong 
correlation among the matrix-elements. In the limit N oo this model has the 
same thermodynamic properties as the p-spin model. The dynamical transition is 
at Tc = 0.536 with threshold energy per spin eth — —1.87. A static transition 
occurs at Trsb = 0.256 and the critical energy per spin is eirsb — —1.936 where the 
complexity vanishes |241l 17^ . The analysis in the mean- field limit gives a rather clear 
"geometrical" interpretation of the two transitions. The phase space is composed by 
an exponentially large (in N) number of components, separated by infinitely large (for 
A'^ — > oo) barriers. Each component is labeled by the energy density e of its minimum 
and the largest allowable value of e is eth [ZHl I242| . Components with e equal to eth 
have the largest (exponentially with N) statistical weight and become dominant, in 
thermodynamic sense, at the dynamical transition T — T^. Since components with 
e smaller than eirsb have negligible statistical weight [ZSlI^j the static transition is 
ruled by components with e = eirs&, i.e., the lowest accessible ones 02 ESI- 

For finite the scenario is different since not only basins with e < eth acquire 
statistical weight, but basins with e > eth with few negative directions |76| may 

t The factor 2 in l'H'A\ is set only for convenience to match the values of all relevant temperatures 
with those reported in the original paper '241 . The Hamiltonian studied in 1741 differs in a factor 4 
from the present definition. To compare the results discussed here with those in |74| temperatures 
and energies must be properly scaled by a factor 4. 




(243) 
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become stable, simply because for finite N there are not enough degrees of freedom 
to hit them. The ROM for finite N has been largely studied during the last years 
| 65l I66L 1571 136L I243j and its behavior has been compared with that of supercooled 
liquids finding a remarkable agreement. 

In the left panel of figure [501 it is shown the average energy- minima of basins (IS) 
as function of temperature for the ROM with = 300 obtained from a Monte Carlo 
simulation. It can be shown [HDIIS] that if the density of states il{E) is Gaussian and 
if the basins have approximately the same shape or are, to a good degree, harmonic, 
then the IS energy density eis ocl/T. The data in the figure can be well fitted by 

eis=e^+eiT-i + 0(r-2) (244) 

indicating that for a relatively large energy range the basins are roughly of the same 
shape. This means that the contribution fy{T, eis) to the free energy density of the 
system is of the form fv{T, eis) = eig + Sfi,{T), with the second term independent of 
eig, i.e., of the component This in turn implies that the effective temperature for 
the ROM is completely determined once the complexity (density) Sc(e) is known 
Indeed from j^ . lf^I^ we have 

1 ^ ds,,{eis) ^245) 
T'cff(eis) deis 

Furthermore, in an aging experiment Tcff depends on time only through eis(t). For 
each time t the effective temperature Tes can be obtained graphically as shown in 
figure 1201 The left panel in that figure shows the average eis energy as a function 
of time in a typical aging experiment. We note the two-regimes decay also observed 
in supercooled liquids |244| . The two regimes are associated with different relaxation 
processes. In the first part the system has enough energy and relaxation is mainly 
due to path search out of basins through saddles of energy lower than T, where T is 
the temperature after the quench. This part depends only on the temperature of the 
equilibrium state from which the system has been quenched. This process stops when 
all barrier heights become of 0{T) and relaxation slows down since it must proceed 
via activated inter component processes. In figure ETI it is shown the response versus 
correlation plot for the ROM. Correlations and responses were computed by projecting 
over the IS, the corresponding FDT also holds in equilibrium as discussed in section 
13.81 The figure clearly show the two-timescales scenario with X = 1 at short times 
and X = T/Tctf < 1 at later times, with Tes in very good agreement with the value 
predicted by H242|l . Also notorious is the fact that the effective temperature shifts 
with time as expected. 



7.2. Spin glasses and other random systems 

As we have explained in section [HI spin glasses represent the most important motif of 
many results regarding FDT violations. Particularly, numerical simulations have been 
the most widespread tools to investigate many aspects of the equilibrium behavior of 
spin glasses that cannot be tackled by analytic means (for a review see |245) '). It is 
usually said that the advantage that numerical simulations offer in the study of non- 
equilibrium properties, as compared to equilibrium ones, relies on the fact that systems 
do not need to be equilibrated. However this observation is naive and deceitful as many 
dynamical aspects cannot be observed in the range of accessible time-scales. Indeed, it 
is widely believed that many dynamical results in spin glasses are suspect because the 
asymptotic dynamical regime, defined as that regime where the dynamic correlation 
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Figure 20. Panel (a): Equilibrium average eis a function of temperature. The 
arrows indicate the construction of the effective temperature Tcfl (eis). Panel (b): 
Average inherent structure energy for the ROM as function of time for initial 
equilibrium temperatures Ti = 3.0 and final quench temperatures Tf = 0.1, 0.2, 
0.3 and 0.4. The average is over 300 initial configurations. The system size is 
N = 300. The lines denotes the two regimes, [sec also Ref. 66;] From EHl- 
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Figure 21. Integrated response function as a function of IS correlation function, 
i.e, the correlation between different IS configurations, for the ROM. The dash 
line has slope = 5.0, where Tf is the final quench temperature, while the 

full lines are the prediction fl^ : T^g{2^''-) ~ 0.694, reff(2l^) ~ 0.634 and 
Tj.g(2^®) ~ 0.608. The dot-dashed line is T^g for t„ = 2^^ drawn for comparison, 
[see also Ref. El]. From 
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length ^ is many lattice spacings, is usually not reached. Establishing whether the 
range of simulated time-scales reaches the asymptotic long-time regime is at the heart 
of a present controversy in the field. Indeed, not by chance, this controversy is 
quite reminiscent of another parallel ongoing discussion concerning the magnitude of 
finite-size effects in the equilibrium properties. Stochastic stability arguments linking 
non-equilibrium with static properties \W{]\ confirm that any strong finite-size 
corrections to the equilibrium properties should manifest as strong finite time-scales 
corrections in dynamical experiments. Unfortunately, a precise theory that quantifies 
(even in an approximate way) these corrections is presently unknown. We will not 
deal here with the difficult issue of ascertaining in which cases simulations do reach the 
asymptotic time regime, but present the evidence on FDT violations for the accessible 
simulated time-scales. 

7.2.1. Spin glasses We begin our tour by reporting the first numerical evidence of 
FDT violations in 3d EA spin glasses j24(ij and their representation in the form of FD 
plots |247j . In |246| it was shown that deviations from the equilibrium FDT appear at 
timescales comparable or larger than the age of the system. In figure|22]we show TRM 
measurements by Franz and Rieger 1247' on the 3d EA model. In those measurements 
the systems starts from a random initial configuration and a magnetic field h is applied 
for time t^,. The field is cut off at t^ and the subsequent decay of the MtrmI^j^iu) 
recorded. In the linear response regime this experiment is equivalent j248j to a ZFC 
setup where the field is initially zero and switched on at t^, the resulting MzFc{t^ tw) 
being given by AfTRM(i, + MzFc{t, tw) — Mpc for t large enough (see where 
Mfc is the equilibrium magnetization. Most of the numerical simulations use the 
ZFC procedure. The most extensive simulations and the most clarifying FD plots, as 
described in sectionOl have been done in d = 3 and d = 4 j^HESUlESlESSESSl- In 
those papers, the authors consider the IRE associated to the global and the spin-spin 
autocorrelation function as described in section 14.31 The system is quenched at low 
temperatures for a time t^ and a small magnetic field is subsequently applied and the 
M{t) measured. Typical FD plots in the 3d Edwards- Anderson (EA) model using this 
construction are shown in figure 1231 There we show S{C) at two different magnetic 
field intensities as well as two different waiting times. FD plots reveal that a constant 
slope X{C) for C < qea is a good approximation to the data (although more accurate 
data in 4d hint at the existence of a curvature in S{C) PSl , 250 ). This constancy of 
the slope X{C) (i.e. the linearity of S{C) in the region where FDT is violated) can 
be also interpreted as complementary evidence of the accuracy of the t/tw scaling in 
the correlation function j254j . Stochastic stability arguments 160, 161 state that the 
dynamical X{C) is related to the static fimction x{q) by the relation P{q) ~ x'{q) or. 



where P{q) is the probability distribution of overlaps between replicas of the same 
system. This identity offers a way to obtain S{C) from equilibrium data. The 
derivative of the relation l|101(l respect to C yields P{C) = —<PS{C)/dC'^. Inverting 
this identity and inserting an estimate of P{C) as obtained from equilibrium 
simulations allows an alternative way to compute S{C). The applicability of this 
method is shown in figure l^ A quantitative evidence of the mean-field character of the 
FDT violations has been also reported by checking the accuracy of the Parisi- Toulouse 
approximation in spin glasses (this approximation states that the order-parameter 
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Violation of FDT in glassy systems 
T = 0.7, U = 10^ 



81 




o 

I 



X 

I 



.6 



0.7, U 



10^ 



1 1 1 IIIM| 1 1 1 IIIM| 1 


1 1 lllll| 1 1 1 IIIM| 1 1 I^IIM 
• 


- 


• 

• 

• 

• 














1 1 1 





I 



.4 



X 

I 

o P 1-9 
X 



T = 0.7, 



10=^ 



I — I 



I I 



ocaaai 



10' IC^ 10'^ 10^ lO^' 10° 



I 

1—1 



.4 



X 

I 

X 



T = 0.7, t„ = 10^ 

TTTTT] 1 I I lllll| 1 I I ■lll| 1 I I lllll| 1 I I III 



ooQcssEsocaiiffl 



I I I I 



10' lO'' lO'^ 10^ 10= 10° 



Figure 22. Zero-field cooled IRF XZFc(t, <™) = MzFc(t,tw)/h = Xdc - 
XFc(*i*tii) (smpty circles) in the 3d EA model plotted as function of the time 
t = T for different values of tw at T = 0.7. For times t = r > deviations from 
the FD relation XZFc(*i tw) = — C{t, t^)) (filled circles) are noticeable. From 



function q{x) is a function of the argument x/T, see |255l I256| for an exposition). 
Within this approximation (which works pretty weU in MF spin glasses) it can be 
shown [2491 1253j that x{t,tw) = x(C) = PS[C) is independent of temperature in the 
region C < qea where FDT is violated. In the SK model it can be proved that this 
function x(C) ~ — C). In general, for any short-range system one can assume 
the following behavior: a) xi^^) = — C) for C > qea when FDT holds and b) 
X(C) = A{1 - C)^ for C < qsA- Multiplying x(C) by T^-'t' with 4> = 1/(1-5) 
one finds that the resulting quantity is a sole function of the argument (1 — C)T~'^: 
T^~'^x{C) = x((l ~ C)T^'^) thereby showing that data for different waiting times and 
temperatures should collapse on a single master curve. Bot in d = 3, 4 a best collapse 
is obtained taking B = 0.41 |249[ I253| which is quite close to the mean-field result 



Violation of FDT in glassy systems 



82 



0.6 



0.5 i 

0.4 ^ 

0.3 - 

0.2 - 

0.1 - 



t„=10**ho=0.05 - 



0.2 



0.4 



0.6 



0.8 



Figure 23. S(C) versus C in tlie 3d EA model at T = 0.7 < Tc ~ 1 for L = 64. 
The continuous line is the prediction obtained from equilibrium data for L = 16 
(averaged over 900 samples) as explained in the text. The straight line is the FDT 
prediction. From ,249^ . 




Figure 24. Scaling plot as described in the text for a cubic lattice, d = 3 with 
L = 64 showing the existence of two dynamical regimes described by two different 
scalings for the scaling function x(^) where x = (1 ~ C)T~'^ . In the FDT regime 
X{x) ~ X while for x > 1, x{^) ~ with B ~ 0.41. From 12531 . 



B = 1/2, see figureEl 

Apart form the EA model many other results have been obtained studying short- 
range versions of the disordered p-spin model. Two different models have been 
considered. On the one hand there is the so called disordered plaquette model |257) 
where spins occupy the vertexes of a finite-dimensional lattice and the interaction 
occurs between p spins belonging to a given plaquette. The Hamiltonian of this model 
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Figure 25. Schematic figure of tlie disordered plaquette model in d = 2. 



reads, 

n = -Y,Jnl[<J. (247) 

□ iea 

where '□' denotes a plaquette (not necessarily a square plaquette) that connects 
different spins. As usual, Jn are quenched variables of zero mean and finite variance 
and fTi = ±1 denote Ising variables. The simplest case, the one considered in |257| . 
is a disordered version of the p = 4 model and consists of a regular lattice of side L 
and dimension D where spins occupy the vertexes and plaquettes correspond to the 
different faces of the lattice. Each face contains four spins and each spin belongs to 
4(^) plaquettes. A schematic picture of the lattice in d = 2 is shown in figure!^ The 
study of the static and dynamics properties of this model revealed that, although there 
was no compelling evidence in favour of a finite T spin-glass transition, the relaxation 
time shows superactivation effects and stretching of correlation functions characteristic 
of fragile glasses. The relaxation time can be fitted both to a VTF law with Tq — 0.65 
or to an exponential inverse temperature squared law with Tq = 0. The study of the 
equilibrium properties confirmed both possible scenarios (Tc ^ Tq — 0.65 or Tc ~ 0) 
but show that, whatever scenario holds, the relaxation time r and the equilibrium 
correlation length ^ are linked by the relation, t = A exp(B^/T) supporting a scenario 
of cooperative dynamics. Also, the trap-like character of the dynamics was confirmed 
by studying the overlap among to identical replicas at t^ but evolving with different 
noises, Q{tw,tw + t). This quantity should coincide with C{tu],tw + 2t) if jumps 
among configurations are uncorrelated or entropically driven. Numerical results are 
compatible with this prediction. Accordingly, the FD plot (figure I26|l measured at 
T = 0.7 (just above the suspected finite Tc) showed strong deviations from the 
equilibrium line x = /9(1 ~ C*) with a FDR x ^ 0.4, and in agreement with the 
one-step pattern characteristic of structural glasses (see section [^J). 

On the other hand there is another short-range version of the p-spin model 
OSgi una EHOI HHD where M spins (sj, ..,s^) occupy the different sites i of 
a cubic lattice. For each two adjacent sites i and j one considers all possible groupings 



Violation of FDT in glassy systems 



84 




Figure 26. FD plot for the disordered plaquette model (with p = 4) in a cubic 
lattice of lattice size L = 20 at T = 0.7 and three different waiting time values. 
Prom 2571 . 



of different p spins that can be formed by taking k spins from site i and p — k spins 
from site j. In an obvious abuse of notation we can write, 

where the sum runs over all possible nearest neighbours {i,j) and all different groups 
of p spins as described above. Again, the J's are quenched variables with zero mean 
and finite variance. Two cases have been considered: M = 2,p = 3 [2581 1261] and 
M = 3,4,p = 4 I^SD) (for M = 2,p = 4 the model reduces to the standard 
EA model |259p . In the first case the model is not time-reversal invariant while 
it is in the second case. FD violations have been measured in this last case |26U| . 
Finite-size scaling studies of the model shows some evidence of a second order phase 
transition at Tc — 2.6 characterized by a divergent spin-glass susceptibility. This result 
is confirmed by a study of the FDT violations in this model that show the existence 
of a non-trivial X{C) characteristic of a full RSB scenario. As for the EA model, 
the X{C) appears to be linked to the static P{q) via the relation P{q) = X'{q), see 
figure 123 The main message conveyed by most of these results is that FDT violations 
are qualitatively and quantitatively well described in the framework of MF theories of 
spin glasses. However, the implications of these similarities must not be taken too far. 
In particular, the already old but recurrent issue about the validity of the many state 
picture in finite-dimensional spin glasses cannot be answered from such a point of view. 
As discussed previously, the precise link between statics and dynamics proposed by 
stochastic stability arguments confirms that the dynamic X{C) is related to the static 
x{q) after introducing a coherence length l{t^) (related to the spin-glass correlation 
length ^{tw) obtained from the two-point replica correlation function) which depends 
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Figure 27. FD plot for the short-range version of the p-spin model (with 
p = 4,M = 3) in a cubic lattice of lattice size L = 16 at T = 2.0 = 0.77Tc and 
tw = 2^7. The continuous line with error bars is the shape of X(C) derived by 
numerical integration of the relation P{q) = x'{q) and measuring the equilibrium 
P{q) for L = 5, the dotted line is the equilibrium line S{C) = 1 — C and the 
isolated cross at C = correspond to the FC magnetization. From 12601 . 

on the waiting time tw according to the relation x{q, l{tw)) = X{C{t^ tw)). A numerical 
study of FDT violations in d = 2 |2()2j shows that such an assumption is indeed true 
and l{tw) oc £,{tw)- More important, the resulting FD plots are extremely similar 
to those found in d = 3,4 and the scaling Ansatz for S{C), as derived from the 
assumption that the low T phase has many states, also works pretty well. However, 
as in d = 2 there is no finite T transition these results show that statements in 
favour of the validity (or not) of the mean-field picture in the accessible range of 
timescales (in off-equilibrium experiments) or sizes (in equilibrium measurements) are 
inconclusive. The stochastic stability property is well satisfied also in the XY model 
at not too low temperatures ^220. where critical fluctuations dominate, the overall 
resulting behavior being quite similar to that of the 3d EA model. From a completely 
different perspective, the overall presence of mean-field aspects in the analysis of off- 
equilibrium data suggests that the many state picture is effectively valid and that only 
for experimentally inaccessible sizes or timescales (therefore irrelevant from a practical 
point of view) the true scenario (whatever mean-field or droplet) is recovered. 

We finish this subsection by commenting some recent results |263| aiming to 
identify and quantify the low-energy fluctuations that locally result in deviations 
from the average QFDT as measured in the bulk. Local deviations from the bulk 
QFDT curve x(C) are the equivalent of fluctuations from the average magnetization 
in a Heisenberg magnet, where transverse fluctuations correspond to low-energy spin- 
wave excitations and the longitudinal fluctuations that modify the length of the 
magnetization vector being the massive ones. Numerical results in the 3d EA model 
show that local correlations Ci{t,tw) and IRFs xii^^'tw) measured over local boxes 
spread over the whole lattice generate a two-dimensional surface p{Ci,xi) with a 
prominent maximum centered around the bulk curve. The contour lines of this density 
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map gently deform along the QFDT bulk curve x(C') and deviations far away from 
that curve appear to be penalized. This study offers the possibility to understand the 
connections between the mean-field character of the FDT violations and the existence 
of deviations due to short scale cooperative processes in an eventual (but yet unclear) 
heterogeneous scenario. 

7.2.2. Other random systems Apart from spin glasses other lattice models with 
quenched randomness have been considered in the literature aiming to elucidate 
whether off-equilibrium studies can tell something about the character of the low- 
temperature phase. Many of the conclusions of these numerical studies need to be 
taken cautiously as no conclusive evidence in support of a given scenario or in refusal 
of other ones is ever reached. 

Let us start the discussion with the ferromagnetic diluted and random-field 
Ising models (RFIMs). The 3d version of both models has been investigated in 
|250J. Simulations in the low T phase and in the Griffiths phase (i.e. the region of 
temperatures between the critical temperature of the pure system and the transition of 
the random system) show that FD plots in ferromagnetic diluted and RFIMs are very 
similar one each other but quite different to those measured in spin-glass systems (see 
the preceding section [7.2.1|l . The former ones are characteristic of a ferromagnetic 
phase with X — while the latter are described by a non-trivial function X{C). 
These studies exclude the possibility of a spin-glass and Griffiths phase in both models 
described by a mean- field like RSB solutions. 

Other finite-dimensional model with interesting behavior is the frustrated Ising 
lattice gas (FILG) defined by the Hamilt Oman, 



where the = ±1 are Ising spins and = 0,1 are occupancies which may take 
the value 1 or depending whether site i is occupied by a spin or not. The sum 
is over near-neighbours on a c?-dimensional lattice. The are quenched random 
variables that may take the values ±1 and /i stands for a chemical potential. The 
average particle density p — y Ui is a monotonically increasing function of fi. 
In the limit J — > oo the model converges to the site frustrated percolation problem 
|264| (a variant of the standard percolation problem where clusters are made out 
of sites connected by non frustrated links). This model has been simulated in 3d 
EMI where different regimes have been singled out. The percolation transition 
occurs at a given value of the chemical potential /ip and manifests in the onset of 
two different relaxational regimes (a fast exponential relaxation followed by a slow 
stretched decay). A second transition is observed at a higher value of /i, p,d > p,p 
where the relaxation time grows extremely fast and dynamics arrests. However, it 
is unclear whether the relaxation time diverges at /i^. Less clear is, in the case 
of the existence of a dynamical singularity at ^d, whether this is associated to a 
thermodynamic singularity. For the FILG 1)249(1 different correlation functions can be 
constructed depending on whether spin variables a or occupancies n are considered. 
A study of FDT violations in this model leads to the following conclusions |266| : 1) 
In the glassy regime fi > dynamics is one-step like, i.e. a two-timcscale scenario 
with two temperatures describes pretty well the relaxational behavior; 2) The effective 
temperature Tcs, as derived from the slope of the FD plots, is pretty independent on 



V 




(249) 
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Figure 28. FD plot for the 3d FILG with lattice size L = 30 at 

T = 1, J = 10, n ~ 5.5. Spin-density and spin autocorrelations (see the text) 
yield FD plots compatible with the same effective temperature. From l^bbl . 



the observable, whether this corresponds to spin variables {C{t,t^) — {a{t)a{ty^))) 
or mixed spin-occupancy variables {C{t,tw) — {<y{t)n{t)fj{tw)n{tu]))) and 3) Teff is 
apparently independent on the waiting time. However, this last result has to be 
taken with caution as the range of waiting times considered in |266j may not be large 
enough % to display such a small effect (compare for instance with the results described 
in section rm figure I^T|l . The observable independence of Toff is shown in figure ESI It 
should be noted that the above scenario is reminiscent of f RSB behavior and occurs in 
finite D rather than in the mean- field limit D ~* oo. Actually the mean- field version 
of the model does not have a IRSB low temperature phase j267| . 

7.3. Coarsening systems 

Although coarsening has been briefly sketched in section Ffi. 61 here we present a more 
detailed account of results. Coarsening systems are the paradigm of systems which 
do not reach equilibrium. In such systems TTI does not hold, and all time-dependent 
correlation functions for large times are of the form C{t,s) — C{L{t)/ L{s)) where 
L{t) (X t^^^ is the typical size of the coarsening regions |268j . The dynamic exponent z 
is characteristic of the universality class of the system and its value depends on whether 
dynamics conserves or not the value of the order parameter. This functional form is 
similar to that found for the long-time correlation functions of glasses in the aging 
regime, and indeed a certain type of coarsening has been advocated as responsible for 
slow relaxation in glasses |269l 12701 1271L I5T] . The difference between the two systems 
only becomes manifest when one also considers the response functions associated 
to the correlation functions. Glasses, such as spin glasses or molecular glasses, are 
characterized by long term memory which results in a non-zero FDR X. On the other 
hand, for systems stochastically stable |16(J[ 1161] the FDR X{C) coincides with the 
static Parisi's function 12461) . For a ferromagnetic system P{q) is trivial 

P{q) ^ 6{q - m^) (250) 

X In fact, the range of waiting times explored does not even cover one order of magnitude. 
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where m = m{T) is the magnetization a temperature T, and hence X is 1 if 
1 > C > TO^ and if C < ni^. The argument can be easily extended to the case 
of few separate phases. Therefore we expect that in systems in which two (or few) 
phases separate, X should vanish for long times, signaling the presence of weak long 
term memory. 

The simplest model displaying domain-growth is a ferromagnetic Ising model on 
a square or cubic lattice of linear size L with a single-spin-flip Glauber dynamics. 
When the system is quenched at time t = from a random configuration {T = oo) to 
a finite temperature T below the critical temperature Tc domains of "up" and "down" 
spins start to form and grow. This is well reflected by the behavior of the two-times 
spin-spin correlation fmiction 

1 ^ 

= E^'^^W^^^^)) (251) 

2 = 1 

which for times t — s s (assuming s < t) is TTI and rapidly decays from 1 to , 
m being the average magnetization at temperature T. Later, for more separate times 
t — s » s the TTI is lost, the aging part of the correlation scaling like 

C,s{t,s)^F 

where L{t) is the typical size of the domains at time t. The calculation of the linear 
response proceeds as usual, i.e., at a certain waiting time tw a small magnetic field hi 
is applied and the induced magnetization is computed. For disordered systems, such 
as spin glasses, the applied field can be either uniform or random. The advantage 
of an uniform field is that averaging over different realizations of the fleld is avoided. 
However, for systems without disorder, such as ferromagnetic systems, an uniform 
fleld would favor one of the phases making it growing faster. In this case a random 
fleld must be used and the correct quantity to measure is the staggered magnetization 

Mit,s)^-J2(''^(*^^^) (253) 

1=1 

where hi is the local (quenched) random fleld, and the over-bar denote overage over 
the fleld realizations. 

In figure we show the curves x(t, i^,) = M{t^tw)/h versus C{t,t^^) obtained 
with a bimodal field distribution hi = zLh \2T2\ . The FDT region and the fiattering 
of the curve are well evident. Two aspect of these curves are worth to be noted. 
First of all the value of the plateau reached by the magnetization decreases as t^ 
increases. Moreover, for fixed t^ the magnetization first grows in the non-aging part 
as h{l — C)/T, then saturates and eventually goes down again. Indeed, the comparison 
|^73j with the equilibrium response function shows that the equilibrium value of the 
response lies rather below the plateau. The study of a soft spin version with Langevin 
dynamics |272l 1273] leads to similar results. 

There are two contribution to the staggered magnetization: one from the domain 
walls, the other from the domain bulks. The difference between the plateau and the 
equilibrium value of M{t, tw) can be attributed to the domain wall response. After a 
time tw the domains have reached a certain typical size, and the domain walls have 
a certain total length. The effect of the random field is to try to fiip some spins. 



L{t) 



Us) 



(252) 
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Figure 29. TM{t,t^)/h versus C{t,t^) domain growth. Left Panel: 2d 
case with Tc = 2.27, at temperatures (from top to bottom) T = 1.7 and 
tw = 200, 400, 800, 2000, T = 1.3 and = 800, T = 1 and t» = 800. Right 
paneh 3d case with T = 2.5 (Tc ^ 3.5), = 100, 300, 600, 1000, 1500. The 
straight hne is Af = I — C: we see that FDT holds at short times t, and the 
violation of FDT with X = at longer time separation. From 12721 



Clearly the flipping is easier at the domain walls where the spins are less constrained 
by their neighbors. As the time proceeds the domains grow increasing the bulks at 
the expense of the total domain wall length. Therefore the contribution from the 
interfaces decreases with time. On the other hand, the contribution of the bulk is 
almost independent of t^ since the effect of a random field on "up" and "down" 
domains is the same on average. Therefore, after the initial quasi-equilibrium growth, 
the total staggered magnetization decreases as t^ increases and, for a fixed value of 
tw it decreases as t increases. Analytical results '2731 show that the aging part of the 
IRF is of the form 



M,g{t,s)r^ A{s)F 



L{s) 



Lis) 



(254) 



where 



A{t) 



1 



m 
m 



d>2 



(255) 



Because X{s) 



I ^g^^y ^ I the results (|255ll explain the slower decrease of the IRF for 
the two-dimensional case observed in figure 1^ The numerical test of the scaling law 
H255|l in 2d is shown in figure QUI 

This scenario leads to the conclusion that coarsening systems do not display a 
non-trivial X{C). Recent results |217| . however, indicate the possibility of non-trivial 
X{C) also in these systems. The motion of the domain wall in the presence of an 
external random field follows from two competing processes: the tendency to reduce 
the interface curvature due to surface tension and the pinning of the domain wall in 
favorable positions introduced by the external field. This introduces a dependence 
on the space dimension since the curvature process, which dominates at large enough 
dimensions, weakens as the dimension decreases. When dimension reaches the lower 
critical dimension the curvature process disappears and the response of the system 
becomes non-trivial |217[ 1214) as it is indeed seen for the ferromagnetic Ising chain 
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Figure 30. Test of the scaling 

of L{s)AIag{t, s) w.r.t. ln(s) with s = 

very well. From HT^ . 



12551 which predicts a linear dependence 
tw The dashed line fits this dependence 



The domain walls may give a large contribution to the response also at the early 
FDT part, but almost exclusively from their deformation on relatively short lengths. 
These fluctuations can be considered thermalized and hence do not spoil the 1/T 
behavior but, on the contrary, contribute to make longer the 1/T slope of the initial 
part in the FD plot. A generalization of the Ising model to include some frustration 
via long-range anti- ferromagnetic interactions has been studied in |274j . Models of 
this type have been proposed to study, among others, the avoided phase transition in 
supercooled liquids |275| and charge density waves in doped antiferromagnets |276j ||. 
They are described by the Hamiltonian 

n = -6j2a.a, + J2^ (256) 

where cr^ are Ising spins, the first sum runs over all pairs of nearest-neighbor sites of the 
lattice, the second over all distinct pairs, and r^j is the distance between sites i and j. 
The parameter 5 represents the local/non-local exchange ratio. It is known |278| that 
for a two-dimensional square lattice the ground state of the model is anti-ferromagnetic 
for 5 < 0.85. For S > 0.85 the anti-ferromagnetic state becomes unstable with respect 
to formations of striped domains. The study of the non-equilibrium dynamics reveals 
a crossover from a logarithmic decay for (5 < (5c ~ 2.7 to an algebraic decay for 6 > Sc 
|274| . However despite this richer scenario, the FD plot leads in both cases to a 
vanishing X, and hence to a coarsening scenario. Similar results are found in other 
coarsening-like systems such as the Migdal-Kadanoff spin glass |279| although in that 
model the definition of the dynamics appears rather tricky as most of the spins occupy 
the deepest hierarchical layer in the model. 

§ See also the discussion in section [6.6l . Similar effects are also observed in kinetically constrained 
models described in Sec. 17.51 

11 These models are of interest for the information storage in ultra- thin ferromagnetic films |277| 
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We conclude this section on coarsening discussing a hard-sphere lattice gas model 
in the spherical approximation originally introduced by Lebowitz and Percus (LP) 
| 28U| . The interest in this model is twofold. On the one hand, its dynamical 
behavior can be solved analytically, hence allowing a detailed investigation of the 
non-equilibrium behavior. On the other hand, it is known that hard-sphere models 
have a fragile-glass behavior |9()j . 

Lattice-gas models are defined on a lattice of finite dimensions. On each site there 
can be a density p{x) of particles, which in the limit of hard-spheres can take only 
the values (empty) or 1 (occupied). In the LP model this restriction is relaxed and 
p takes any continuous value allowed by a spherical constraint: 

ai = ^p(a=)2-^p(a;)- (257) 

X X 

There is the additional restriction that the density-density correlation function 
between nearest neighbors vanishes. This is added to mimic some kind of extended 
hard-core: 

a2 = 5]^p(a;)p(a; + q) = (258) 

X q 

where q are the vectors that join a lattice site x to its nearest neighbors. In [^3 the 
following Langevin dynamics for an open system has been studied 
d d 

- p(x\ t) = p- ^-^n[p(x'),t] + j^(x', t) (259) 

where p is the chemical potential, 77 the thermal noise, and TL the Hamiltonian 

H = Ao(t)^[p(a;)2-p(a;)] +Ai(t)^^p(a;)p(a; + g) (260) 

X X q 

Strictly speaking in this model there is no energy, but only entropy. The role of 
Xo(t) and Ai(t) in the Hamiltonian is to make the dynamics to fulfill the constraints 
H257|l and H258|) at all times. The non-equilibrium dynamics of this model shares a 
large number of features with that of the spherical SK model |281lll49Tll5U| where the 
dynamics is driven by the macroscopic condensation of the system onto the disordered 
ground state. A positive chemical potential would increase the local density, thus 
starting from an empty state the local density relaxes toward the equilibrium value. 
The relaxation can be divided into two regimes, the first one where the system is filled 
in a spatially uncorrelated way. The typical time for this process is order t* ^ 0(1). 
It is only later that slow relaxation starts when the system is spatially correlated 
and needs to reorganize large regions in order to increase its density. The model has 
no built-in disorder, and the slowing down is a purely entropic, direct consequence 
of the decrease of the number of available configurations imposed by the the short- 
range constraints. When the temperature is below the critical temperature Tc the two 
motions have well separated timescales and the two-times correlation function shows 
the usual two-steps form with a first part TTI and the second part scaling as t/t^, 
see figure 121 (that plot corresponds to the 3d model at T = 0.1. Different waiting 
times from top to bottom are t^ = 10000,1000,300,100,30,10,3,1) The study of the 
FDR reveals that there is no anomaly in the response function |31| and X vanishes 
for values of C below the plateau value of the correlation. The glassy scenario of 
this model corresponds to that of phase-ordering kinetics with non-conserved order 
parameter. Similar results have been reported in spherical models with long-range 
ferromagnetic interactions ,282| . 
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7.4- Non-relaxational driven systems 

We have been underlining throughout this review that aging systems are characterized 
by a non-equihbrium behavior with lack of TTI and by the presence of FDT violations. 
The relevant parameter which controls the aging non-equilibrium state is the waiting 
time or time elapsed since the system was quenched. However, there is another way to 
generate a non-equilibrium state that can be characterized by a given timescale in the 
same fashion as the waiting time characterizes the aging state. For instance, adding a 
time-dependent perturbation of frequency w to a time-independent Hamiltonian. In 
the regime where Lotw ^ 1 the perturbation oscillates slow enough to probe only slow 
processes occurring at time-scales t ^ tu,. While in the opposite regime wi^, S> 1 
the oscillatory perturbation probes fast relaxation processes occurring at timescales 
T <^t^,. The line Lot^^ ~ 1 marks the onset of glassy behavior, the shape of this line 
depending also on the intensity of the perturbation. 

Driven systems have advantages when compared to aging systems. One of the 
most important differences is that driven systems, in the stationary state, are TTI but 
FDT is still violated. Again the concept of an effective temperature can be introduced 
as a measure of these violations. However, as the non-equilibrium stationary state of 
driven systems can be described by the intensity of the driving force, they are more 
experimentally accessible than aging systems, where the waiting time appears as an 
external parameter difficult to control. For this reason, it has been advocated that 
experimental measurements of FDT violations and the effective temperature should 
be done in driven systems rather than aging systems. 

There are many ways to put a system into a driven stationary state and these 
have been investigated in the literature for different types of models. Driven systems 
can be classified into two main groups: 

• Sheared systems. In this case one considers systems where, in addition to 
conservative forces, other non-conservative forces (i.e. that cannot be derived from 
a potential function) act upon the system. In these systems the non-conservative 
forces maybe time-dependent or not. In both cases, the non-conservative forces 
do net work along a given closed dynamical path. This implies that energy power 
is continually supplied to the system by the driving force. The parameter which 
describes the stationary state is the intensity of the shearing or driving force 
which we will identify by the symbol e or 7. These systems include models that 
violate the action-reaction principle (such as models with asymmetric couplings) 
and sheared fluids. These systems are described in section 17.4.11 

• Tapped systems. In this class, systems are driven to a non-equilibrium state by 
a time-dependent force which, however, derives from a time-dependent potential. 
This means that the driving force, if constant in time, does not exert work 
upon the system whatever its intensity e. This class of systems includes spin- 
glass models in a oscillating magnetic field and tapped granular systems where 
the relevant parameter is the frequency of the driving force. These systems are 
described in section 17.4. 21 

7.4-.1. Sheared systems Studies of models described by non-conservative forces go 
back to the study of neural network models described by synaptic interactions that 
are non-symmetric |283| . This has inspired future investigations of disordered models 
where couplings among spins include an important degree of asymmetry. In |284l I285| 
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it was considered the relaxational dynamics of the p — 2 spherical spin-glass with 
pairwise interactions given by = jfj + eJ^ where jfj = J^^ denotes the 
symmetric (therefore conservative) part and = —J^ stands for the antisymmetric 
(therefore non-conservative) part. It was shown that any finite amount of asymmetry 
is enough to destroy the spin-glass phase. The relaxational time of the system was 
found to diverge as for e ^ 0. Importantly, this result suggests the following 
general scenario: whatever the intensity of the non-conservative force, the stationary 
state has a finite relaxation time, therefore the stationary state, although of non- 
equilibrium nature, must be TTI (and therefore, correlations and response functions 
do not display aging). Subsequent investigations have confirmed this result showing 
that this is a generic feature of driven systems. 

Among this family of asymmetric spin-glass models, one which has been 
intensively investigated in the past years is the p-spin spherical spin-glass Hamiltonian 
H150() . H151() with asymmetry in the interactions. The model is defined according to 
the following Langevin dynamics P^IW) 

^ = -r{t)<j,{t) + F,{{a}) + rj,. (261) 
This dynamics is the same as (|153|l but now the force Fi{{a}) is replaced by, 

0(Ti 



5H 

5(Ti 



E Kr'-"'-''^n-^j.-. (262) 



e 

ji<--<jk-i 

where TC is the Hamiltonian given in H150|l and the driving force -F7'"™({o'}) describes 
a fc-spin interaction term where the couplings ^^I'^^'- '^'^-i ^re uncorrelated among all 
permutations of the k different indexes {i,ji,j2, ••, jfe-i), 



j^n,32,--,3k-i j^ji j2,--,jr-i,*,jr+i,-- jfc-i _ Q (263) 

These models show the following behavior. In the regime above the mode-coupling 
temperature Tc, the relaxation time is finite for the unsheared model e = 0. Therefore 
both TTI and FDT hold in the stationary state. A small driving force e > puts 
the system in a new stationary state where TTI holds but FDT is violated. The 
numerical analysis of the mean-field equations 287 reveals that both FD plots and 
the value of the FDR are very similar to those found in aging systems. Figuremishows 
these quantities above Tc for the model with parameters k — p = 3. Below Tc a new 
phenomenon, called "shear thinning" occurs. At e = the relaxation time diverges 
(as is common in mean-field (MF) models where activated processes are neglected, 
see discussion in section fG.lll . However, as e > the relaxation time becomes finite 
and decreases with e (shear thinning). Again, numerical analysis of the mean-field 
equations reveals that for finite e the resulting FD plots are the same as for an aging 
system with a waiting time t^ given by - e-"(^) where a{T) is a temperature 
dependent exponent that takes the value 2 at Tc and slowly increases as T decreases. 
This scenario, as derived from the study of MF spin-glass models, has been confirmed 
in numerical studies of binary mixtures of Lennard- Jones (LJ) sheared fluids (see 
section !?. 1.2(1 . In a series of papers, Barrat and Berthier ,2.88 , 289 , 2Q2I have shown that 
driven short-range systems display the same features as their equivalent MF disordered 
models. These similarities have been confirmed by other studies that measure the 
temperature dependence of the shearing rate at which the average potential energy 
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Figure 31. FD plot for the model j^MJ , , OMJ with k = p = 3 and 
T = 0.613 > Tc = 0.612. The full Une is the e = equihbrium curve. Dashed 
lines correspond (from bottom to top) e = 0.333,0.143,0.05,0. The inset is the 
FDR as function of the intensity of the perturbation. From |287| . 



deviates from its equilibrium value |291| . The similarity between LJ models and MF 
spin glasses is striking concerning rheological properties. In a sheared fluid the velocity 
field V = 'jyex induces a stress a^y that is well described by the phenomenological 
law: axy ~ (Tq + 07". This law is commonly found in rheological systems 292 . The 
viscosity is defined by 77 = cr/7: for Newtonian fluids n — 1 and (Tq = 0, the viscosity 
is therefore independent of the stress. However, this is known to be inaccurate as 
many complex fluids show transport coefficients (such as the viscosity) that depend 
on the shear rate (for a discussion of these effects in the framework of non-equilibrium 
thermodynamic theories see |293l 129^ '). In figure 1^ we show the flow curves for the 
viscosity as function of the shear rate for both the LJ fluid and the MF spin-glass, the 
viscosity rj and the shear 7 in the LJ model corresponding to the terminal relaxation 
time Ta and e/Ta in the MF spin glass respectively. The similarity is noteworthy. There 
are two regimes depending whether T > Tc or T < Tc- In the first regime the fluid is 
Newtonian at low shearing rates so the viscosity is shear independent. In this regime, 
standard non-equilibrium thermodynamics 295^ is applicable. However, for T < Tc 
the fluid is non Newtonian and the viscosity diverges at zero shear, rj ~ 'y~"(^) with 
a(T) between 2/3 (at Tc) and 1 (for T — > 0) in agreement with the aforementioned 
results found in the MF spin glass. The central question addressed in j288l 12891 1290) 
was the dependence of the resulting FD plots on the type of observable used as a 
perturbation. They considered 5 different classes of observables: 1) The incoherent 
part of scattering functions (corresponding to single particle density fluctuations) , 2) 
the coherent part or the correlations, 3) The "chemical" observables associated to 
correlations of a single species of particles in the binary mixture, 4) The mean-square 
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Figure 32. Viscosity as a function of the shear rate for the LJ (top 
panel) and MF spin-glass models (bottom panel). The viscosity rj and 
shear rate 7 in the LJ fluid correspond to the terminal time Tc and e/ra 
respectively in the MF model. Tc ~ 0.435 in the LJ model and ~ 
0.612 in the MF model. Temperatures are (from bottom to top): (bottom 
panel) 0.9,0.8.0.7,0.64,0.62,0 .613, Tc,0.6115,0.58,0. 45,0.3 ,0.01 and (upper panel) as 
indicated in the box. From (left panel) and 1^^ . 

displacement of particles associated to a constant small force transverse to the flow and 
5) the stress in the transverse direction after compression of the box. 1), 2) and 3) were 
measured at different wavevectors. In all cases, the effective temperature was found 
to be the same within numerical accuracy. A compendium of their results is shown in 
figure 1231 There have been other studies on driven systems that have investigated the 
shear thinning effect, i.e. whether a driving force of rate 7 stops aging up to a timescale 
proportional to I/7 by restoring TTI. Corberi et al. |29tij have investigated coarsening 
models in the presence of a driving force. In particular, they have considered a N- 
component ferromagnetic model with non-conserved order parameter in the large N 
limit. The authors find that, above the ferro-paramagnetic transition temperature Tc, 
the inverse of the driving rate I/7 sets the time scale after which aging stops and the 
system becomes TTI. However, below Tc, contrarily to what is observed in MF spin 
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Figure 33. 14 FD parametric plots in the LJ model showing the independence of 
the effective temperature on the observable. The numbers in the key of the figure 
indicate different values of the wavevectors. 1.87 refers to class 1 in the text. 
coh stands for coherent or global scattering functions (class 2). A and B refer 
to the two different chemical species, (A+B) standing for incoherent scattering 
observables independently measured for A and B species (class 3). se// indicates 
the use of the mean-squared displacement (class 4). stress refers t class 5. The 
value of the effective temperature is compatible with Teg = 0.65 > T = 0.3 for 
all cases. From nM . 



glasses, aging never stops even in the presence of shearing. The origin of this difference 
is unclear. It could be related to the fact that in coarsening systems the mechanism 
responsible of aging is different than in MF spin glasses or LJ glasses. This difference 
can be traced back to the absence of complexity in the coarsening model as compared 
to the other cases. Moreover, this behavior has to be contrasted with what is found 
in tapped systems such as MF spin glasses in an AC field |297| where aging survives 
below Tc in a certain range of magnetic field values (see further discussion on this 
model in section 17. 4. 2|l . Other, studies have considered models such as vortex glasses 
with random pinning centers in two dimensions |298| where the external driving force 
is uniform over all vortices. Also in this case FDT violations appear to be described 
by an effective temperature related to the slow motion of vortices. Trapping dynamics 
due to pinning defects has also been considered in the study of a driven classical 
particle submitted to a force [2991 l3U0j 

A common feature of all studies reported here, that has not been emphasized 
enough in the existing literature, is that measured correlations must be transverse 
to the direction of the shear flow to yield a meaningful effective temperature. 
FDT violations for longitudinal observables are apparently not described by an 
effective temperature. This suggests that neutral observables describing an effective 
temperature are restricted to the transverse direction, a fact that lacks yet a clear 
explanation. 

7.4-2. Tapped systems As already said, another class of driven systems correspond 
to those where the external force derives from a time-dependent potential that pumps 



Violation of FDT in glassy systems 



97 



energy into the system. This includes also driven granular media that have recently 
received considerable attention. 

In I297| . the authors have studied the MF p-spin glass model in an oscillating 
AC field. This problem is interesting as it shows a behavior different to that observed 
in sheared MF spin glasses or LJ glasses (see section r7.4.1|) . The model considered 
was again H15()|l , p51|l but in the presence of an uniform AC field of frequency uj and 
intensity h. The phase diagram of the model is described by three parameters T, h, oj. 
Below Tc, at fixed w, the glassy phase survives below a critical field h^{T) meaning 
that, contrarily to sheared systems, a small driving force does not destroy the glassy 
phase and aging never stops, hence TTI is not restored. This justifies the use of AC 
fields in experiments to explore the aging regime within the linear response region. A 
striking result in these studies is the presence of reentrant behavior at constant field 
h, as OJ is varied, indicating that liui^^o /i* (T, w) < /i* (T, oj = 0), a result that needs 
yet to be clarified. All over the glassy phase of the driven model FDT is violated with 
the characteristic FD plots of the corresponding relaxational model. 

Granular systems may present very slow processes, analogous to what is seen in 
glasses. Similar to what has been done for other glassy systems one can try to describe 
the dynamics of the slow degrees of freedom through an effective temperature defined 
from the FDR The first example we consider is the kinetically constrained 

Kob- Andersen model |3l)l| that, even if very schematic, reproduces rather well several 
aspects of glasses I302j and of granular compaction |303j . Although next section 1731 
is devoted to kinetically constrained models we prefer to describe the Kob- Andersen 
model here since it is a good model for granular media. This model consist of N 
particles in a cubic lattice, with periodic boundary conditions. There can be at most 
one particle per site. Apart from the hard-core repulsion there are no other static 
interactions among particles. At each time step a particle can move to a neighboring 
empty site of a three-dimensional lattice only if it has less than four neighbors in the 
initial and final position. In its simpler version there is no gravity, but the system is 
subject to a constant pressure on its surface, obtained by adding or destroying particles 
on the topmost layer with a chemical potential /i. The dynamic rule guarantees 
that at equilibrium all configurations of a given density are equally probable. Indeed 
when the density are of the order of the the jamming density pg ~ 0.88 the particle 
diffusion becomes extremely slow, due to the kinetic constraints, slowing down the 
whole compaction process. As done for glasses, one can try to describe the slow 
non-equilibrium motion through an effective dynamical temperature Tee which can 
be defined via a generalization of the Einstein-Stokes relation between the diffusion 
coefficient and the viscosity to a non-equilibrium (aging) situation: 

^(r{t) - rit^)) ^ -^±-{{r{t) - r{t^))') (264) 

where r is the particle position, / a (small) perturbing field, t > t^ two widely 
separated times, and B{t^tw) the mean square displacement: 

B{t,t^) = {{r{t) - r{t^)f) 

, 3 N 

= ^EE([^'"w-^"(^-)]') (265) 

a— 1 i—1 

The linear response function can be computed numerically by applying a small random 
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Figure 34. Einstein relation in the Kob-Andersen model: plot of the mobility 
x{t,iw) vs. the mean-square displacement B(t,tw) (data shown as circles). The 
slope of the full straight line corresponds to the equilibrium temperature (T = 1), 
and the slope of the dashed one to Edwards' prescription obtained from figure l35l 
at = 0.848. From WEI . 



perturbation at time t^ of the form |304| : 

3 N 

^^^ = ^EE/>" (266) 

a— 1 i—1 

where /f are independent quenched random variables which take the value ±1 with 
equal probability. With this choice the IRF is defined as: 

3 N 

x(i,i.«) = ^EE((/"^'-^"W) (267) 

a— 1 i—1 

where Ar{t + tw) is the difference between the displacement of the same particle in two 
identical copies of the system one evolving in presence of the external perturbation 
and one without. From (|2fi4|l x is related to B through; 

Xit,t^) = ^ X{B)dB 

^-^B{t,t^), if X{B) = m ^ const (268) 
2Tcff 

with Toft — T/m. The numerical simulations show a two regime scenario |3U4| similar 
to what is seen in glasses, i.e., for t — t^ smaller than t^ m = 1 and FDT holds, while 
for t — t.uj tyj FDT is violated with m < 1, figure|^ As discussed in section [^31 the 
value of the observables attained dynamically in a granular system could be computed 
from the usual equilibrium microcanonical distribution at the corresponding density 
p restricted to the subset of blocked configurations (Edwards ensemble), i.e., only to 
those configurations in which every grain is unable to move. Much alike to the IS 
analysis in glassy systems f section 14.511 . the Edwards ensemble leads to the definition 
a temperature, called the Edwards' temperature, as 

1 1 dsEdwip) 



7Edw(p) P dp 



(269) 



Violation of FDT in glassy systems 



99 



where SEdw(p) is the Edwards' entropy density obtained from the logarithm of the 
number of blocked configurations of given p, see (|135|l . The chemical potential fixes 
the dimension |3U6[ I3U5| . 

The Edwards' entropy for this model has been computed in |3U6[I3(J5| through the 
use of an auxiliary model in which each particle has energy equal to one if the dynamic 
rules allow it to move, and zero otherwise. The auxiliary energy i?aux is hence equal to 
the number of mobile particles. The configurations of the auxiliary model are sampled 
with a Monte Carlo procedure with non-local moves at the auxiliary temperature 
l//3aux- These non-local moves have nothing to do with the true dynamics of the 
original model, and the auxiliary model is not glassy. This allows to obtain equilibrium 
properties such as the auxiliary energy density Caux from which the entropy density 
can be obtained via thermodynamic integration: 

(^aux,p) 

eaux(/3Lx,P)rf/3aux (270) 

since Saux(0,p) = Scquii(p)- The blocked configurations at a given density can be 
computed performing a simulated annealing down to /3aux — > oo of the auxiliary model 
at fixed particle density. The Edwards' entropy density is then obtained as 

SEdw(p) hm Saux(/?aux,p) 

/3aux— ^OO 

/"OO 

= = Scquil(p) - / eaux(/3aux,P)rf/3aux- (271) 

Jo 

since e(/3aux7p) vanishes for /3aux — > oo. The Edwards' entropy is shown in figure 
1351 as function of density. We are now in a position to compare the long-time non- 
equilibrium results with those obtained from the Edwards' measure. In figure 1341 it 
is shown the mobility versus the mean square displacement B{t,tw)- The 

agreement between Tcff and TEdw is clearly satisfactory. 

Similar results have been found ^306 305 307, also for models with geometrical, 
rather than kinctical, constrains, the so called "Tetris" models |308[ I309| . These 
models are defined on a two-dimensional lattice with particles of randomly chosen 
shapes an sizes. The only constraint is that particles cannot overlap: for two nearest 
neighbor particles the sum of the arms oriented along the bond connecting them has 
to be smaller that the bond length. 

Tapped granular matter (310k i311) has been considered in several models such as 
one-dimensional model ferromagnet |312| showing that the Edwards measure provides 
a very good description of the stationary state. A similar agreement is found in 
simulations of sheared granular matter |313j . Real compaction occurs in presence of 
gravity j314j . The gravity introduces an extra term in the Hamiltonian 

N 

Hg ^mg^h, (272) 

1=1 

where g is the gravity constant, m the particle mass and hi the eight of i-th particle. 
The simplest way of including gravity in the above models is assuming that particles 
can move up and down (if they respect the geometrical or kinetical constraints) 
with different probability, p = min[l, exp(— mgA/i/T)] where Ah = —1,0,1 is the 
elementary vertical displacement. A closed boundary is situated at the bottom of 
the system. The control parameter x = exp(— mg/T) represents the "vibration". 
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Figure 35. Edwards entropy per particle of the Kob- Andersen model vs. density 
(full curve) . For comparison we also show the equilibrium entropy (dashed curve) . 
At high enough density the curves are indistinguishable, and join exactly only at 
p = 1. The slope of the tangent to SEdw(p) for a generic p allows to extract 
TEdw(p) from the relation Tsdw (p) ''"fp " = P- From ,305^ . 

The presence of gravity introduces a preferred direction in the diffusive motion 
(downward) which in turn may produce inhomogeneities in the vertical density profile 
making the situation more complicated, moreover horizontal (transversal) and vertical 
(longitudinal) quantities must be treated separately. 

Inhomogeneities are only along the vertical direction, so transversal observables 
when measured well inside the bulk are not too sensitive to the detailed form of 
the density profile. Indeed if homogeneity of the bulk is imposed, the dynamical 
temperature obtained from the FDT ratio for the horizontal displacement-mobility 
coincides with the Edwards' temperature |3(J7| . 

The analysis of the FDT ratio for the longitudinal motion [3151 I316| also 
shows a two-slope scenario, however the comparison with the Edwards' measure 
is more complex due to the presence of inhomogeneities. In |316| . based on the 
observation that the density profile far from the top and bottom layers is rather flat, 
the effective dynamical temperature for longitudinal observables and the Edwards' 
temperature have been compared assuming an homogeneous density. The comparison, 
however, reveal strong deviations when inhomogeneities are stronger indicating that 
inhomogeneities of the density profile must be included. This could be done by using 
the recent introduced restricted Edwards measure |116ll317l| . a route not yet explored. 

We stress that the Edwards' measure is constructed from a white sampling of 
blocked configurations, hence it reproduces the physical quantities at large times only 
if this condition is satisfied by the long-time dynamics. In other words the Edwards' 
measure can be inappropriate even though the system presents a slow dynamics. This 
is for example the case of the three-dimensional Ising model in a weak random magnetic 
field |306:. 305 . Here the long-time configurations at low temperatures are made of 
domains of "up" and "down" spins of similar volumes, so that the global magnetization 
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is zero. This is quite different from either the equiUbrium or blocked configurations 
since both of them are magnetized. Other studies have extended this method to the 
study of kineticaUy constrained models (KCMs) (see section I7.5|l such as the Kob- 
Andersen model |318| . 

7.5. KineticaUy constrained models 

A category of statistical models that has received considerably attention during the 
last years are kineticaUy constrained models (KCMs). In a nutshell, KCMs are models 
with trivial thermodynamics but complicated dynamics arising from a, put by hand, 
set of forbidden transitions in configurational space. Allowed transitions are selected 
according to a given transition rule that constrains the dynamics of the system. There 
are different ways to justify KCMs as valuable models for glassy dynamics. One 
can think of KCMs as effective models in which the slowest degrees of freedom are 
idealized as quenched variables that manifest as dynamical constraints. Therefore, 
the statistical variables of KCMs can be seen as the fast dynamical variables that are 
slaved, through the dynamical constraints, to the motion of the slowest ones. These 
constraints are of local nature as also is the interaction between the original degrees 
of freedom. 

KCMs display most (if not all) of the features characteristic of glassy systems 
including slow relaxation, activated behavior, cooperativity and non-equilibrium 
phenomena such as aging and FDT violations. Recently, a review on KCMs has 
been written that covers all these aspects. For this reason, here in this review we will 
not dwell much on discussing FDT violations on these models but content ourselves 
on underlining some of the most important results. We refer the reader to for a 
comprehensive and exhaustive survey. 

The most representative families of KCMs are spin- facilitated models |319l I320j , 
lattice gases [3011 EHl EOS E12 EOll, topological cellular models 3231 EH and 
plaquette models |325l 13261 132 7j . We already discussed in section 17.4.21 the Kob- 
Andersen model as an example of kineticaUy constrained lattice gas for granular 
matter. There are two aspects of KCM's that makes them specially attractive from 
a theoretical point of view. On the one hand, the thermodynamics of KCMs is 
straightforward, there is no underlying thermodynamic singularity, and even more 
astonishing, often the model is purely non-interacting from the point of view of its 
equilibrium properties. On the other hand, the non-interacting character of the energy 
function, entails that the slow dynamics is described by representing the original 
variables in terms of a new set of effective variables. Relaxation can be visualized 
as a dynamical process where defects diffuse (in either a free or cooperative way) and 
annihilate each other, leading to a more tractable problem from the analytical point 
of view. 

We begin our tour presenting the simplest among these type models. Our 
intention is to illustrate with an example what type of models KCMs are. Maybe 
the simplest KCM is the spin facilitated model (SFM) introduced by Fredrickson and 
Andersen (FA) model |319l I320| consisting of free "spins" in a field. The model is 
defined by, 



where rii = 0, 1 corresponding to two possible orientations of the spin {n — 1 up, n — 



N 




(273) 



i=l 
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down) f that occupy a D-dimensional lattice. This model would be trivial if it were 
not by the dynamical rules that describe how spins can flip. In the standard SFM 
model transitions rij — > 1 — are allowed if at least / among the possible nearest 
neighbours are up or n = 1. This model defines the /, d-SFM that shows different 
behaviors according to whether f — I (diffusive) of / > 2 (cooperative). Variations 
of the SFM include KCMs with directed constraints such as the Id East model \[^'2H\ 
or the 2d North-East model |82Hj . defined again by (j2ZS|), but where a spin Ui can 
flip only if the spin on its left is up (East model) and if both spins, that are nearest 
neighbours in two fixed orthogonal directions (North-East model), point up. In Id, it 
is also possible to define a model that interpolates between the 1, 1-SFM and the East 
model pn| . 

In what follows we concentrate our attention in FDT violations in diffusive KCMs 
where some understanding has been recently gained. A generic aspect of KCM's is 
that dynamics is determined by the motion and annihilation of isolated defects. Slow 
dynamics strictly occurs at temperatures close to T = (and timescales larger than 
an initial fast transient). The relevant relaxing variable is the number of isolated 
defects c{t) (e.g. in the 1, 1-SFM this corresponds to the number of up spins). This 
has interesting consequences [^'27\ 12171 \',VM\ in the behavior of the IRF. The IRF in 
the presence of an external random staggered field x{t, tw) is then proportional to the 
product of the number of isolated defects c{t) and the individual local response typical 
of one defect XcS- This local response is TTI as the defect is isolated leading to. 



where r(T) is the relaxation time of the defect. This relation has two important 
consequences: 1) The function x{t,tw)/c(t) is TTI, a result which is characteristic of 
coarsening systems (see below) but not of other (cooperative) glassy models where the 
IRF is a genuine aging function depending, for instance, on the ratio h{t)/h{tw) and 
2) x(^j^to) has a maximum as function of t — t^, when the defect concentration (that 
decreases with time) compensates the growth of the monotonically increasing effective 
response XcS- 

Non-equilibrium measurements for the 1, 1-SFM |7(JII332| and related models such 
as the 2d triangle model 333 , 334 or topological cellular models _323_. 324, show an 
IRF (|?7^ displaying a maximum as a function of t — t^ and leading to awkward FD 
plots. One example is shown in figure Buhot and Garrahan j327| have explained 
how to recover well defined FD plots by using H274|l and plotting x{t, tw) as function 
of the difference Cc{t,t) — Cc{t,t^) where Cc is the standard connected correlation 
function C^t, s) = {{n{t) - c{t)){n{s) - c(s))) with c{t) = {n{t)) for the 1, 1-SFM. In 
this case, the resulting FD curve corresponds to a straight line (corresponding to the 
equilibrium result) casting doubts on the usefulness of FD plots for non-cooperative 



The origin of the maximum in the IRF has also been considered in the Id RFIM 
with infinite- ferromagnetic coupling J — oo |331j . Actually, in this limit the model 
turns out to be a KCM as transitions are only allowed on the spins sitting on the 
interfaces of the ferromagnetic domains, i.e. a spin can flip only if its left and right 

t The variables n cannot be considered as occupancies and 12731 does not define a lattice gas model. 
In lattice gas models the total number of particles (i.e. the energy E in 12731 ') is conserved while in 
the SFM it is not. See 1151 for a throughout discussion. However, we will continue denoting the spins 
by n. 




(274) 



models. 
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Figure 36. FD plots for the 1, 1-SFM 1^7^ with = 10^, quench temperature 
T = 0.3 and different waiting times (see the box). The straight line is the FDT 
relation. From 1701 . 

neighbours point in different directions (meaning that the ferromagnetic contribution 
to the local field acting on that spin vanishes). If quenched at low enough temperature 
this model displays a coarsening behavior with two well separated regimes. If the 
average distance between interfaces (or average domain length) L{t) is smaller than 
a length scale Lg ~ (T/ho)"^ (where ho is the mean square deviation of the intensity 
of the random field) then dynamics is diffusive. However, if L{t) > Lg dynamics 
becomes activated of the Sinai type. In the Sinai regime the IRF is well described by 
the relation H274(l . xit^^w) = Xa.git,tw) = c{t)xes{t,tw) leading to a maximum of the 
IRF at intermediate times where L{t) ~ Lg. This relation sets a crossover timescale 
where pure diffusion takes over to Sinai diffusion. Comparing this result with (|7()|l we 
note that there is no stationary contribution to the IRF because for J = oo thermal 
fiuctuations within domains are suppressed. Again this leads to FD plots similar to 
figure|2ni For the Id RFIM in the asymptotic long-time limit L{t) ^ Lg the xi^jtw) 
decays to zero, a property required to establish a link between static and dynamic 
properties. As remarked in section [6.51 however, this property does not hold for the 
Id Ising model. 

A description of the IRF for the 2d plaquette model along the same lines of ()274f) . 
but modified to account for the diffusion and annihilation of oscillating pairs of defects, 
has been presented in |327| . The resulting FD plots have been shown to display the 
characteristic two slope curves of two-timescale systems. However, it is unclear what 
is the physical meaning of the piece of the curve with slope X < 1 and whether indeed 
Teff = T/X can be considered as a thermodynamic temperature. Future research will 
show what is the true meaning of these violation factors in KCMs of the coarsening 
type. 

FDT violations in the glassy regime of KCMs are not simply described within a 
thermodynamic IS formalism, the Kob- Andersen model perhaps being an exception. 
In fact, the blocked states in the 1, 1-SFM and the East model generate identical 
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configurational entropies [70i albeit they show very different dynamics (diffusive 
and cooperative respectively). The applicability of thermodynamic non-equilibrium 
concepts to models with trivial equilibrium thermodynamics remains an open question. 



8. QFDT: the experimental evidence 

Any valuable physical theory must be successfully challenged by experiments. 
Traditionally, the most direct way to experimentally access FDT violations is through 
noise measurements. In the frequency domain the FDT 1)601) corresponds to the 
Nyquist formula that relates the power spectrum of an observable to the imaginary 
part of its susceptibility. With the same notation we used in sections l2l3l we can write 
define the power spectrum Sa,b{^) and the complex susceptibility xa,b(w), 

1 

Sa.b{oj)^- CA,B{t)exp{iLut)dt; (275) 

...M^ilf (2™, 

where A{uj) (and analogously B) is given by, 

A{uj) ^ / A{t) exp{iujt)dt (277) 



If Xa,b{^) = x'a b(^) + b(^) the Nyquist formula (|133|l reads 

SA.B[i^) = (278) 

The power spectrum can be experimentally measured considering the case where the 
external perturbation couples to the measured observable A = B. In the experimental 
protocol, the time evolution of the observable A{t) and the out-of-phase susceptibility 
x'a Bi^) are recorded. The power spectrum is given by SAi(^) oc (|A(aj)p). This 
allows to verify the Nyquist relation 1)278(1 . Typical noise experiments are the 
observation of electric voltage fluctuations or the motion of a Brownian particle. 
Other measurements included sample to sample fluctuations of the resistivity in small 
samples 13351 133|ij. Observations of magnetic noise were successfully undertaken in 
spin-glasses |337ii338i i339 , more than ten years ago, and despite of the extremely low 
noise signal characteristic of magnetic systems. Globally, these experiments show that, 
within the accessible window of frequencies and times, no systematic FDT violations 
are observed. The out-of-phase susceptibility x"('^) associated to x(u;) is practically 
frequency independent and the power spectrum shows the characteristic 1// behavior. 
The negative result of these experiments points out one of the most important 
difficulties encountered in these type experiments. As they cover the frequency range 
10~^, 10^ Hertz, these frequencies are much larger than the inverse of the aging time, 
therefore only the locally equilibrated regime cot ^ 1 is explored. Only seldom the 
"slow" regime ~ 1, where FDT violations are expected, is measured. Very recently, 
direct measurements of FDT violations in the time domain have been reported for 
insulating spin glasses |340) . FDT violations have been measured beyond the quasi- 
stationary regime and experimental FD plots have been found to be consistent with 
results obtained by numerical simulations in spin glass models (see section I7.2.1|) . 
However the value of the effective temperature obtained in these experiments is much 
larger than the annealing temperature at which the system is equilibrated before the 
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Figure 37. FD plot obtained by measuring voltage autocorrelations and 
relaxation susceptibility in the insulating spin glass CdCr\,T Ino.^Si quenched 
from T = 1.2Tg down to T = O.STg (Tg = 1G.2K). Effective temperatures are 
around ZOK and much larger than the annealing temperature 19AK . From 13401 . 



quenching takes place, casting doubts on the meaning of the effective temperature as 
a thermodynamic temperature in these experiments. Moreover these measurements 
reveal that, within the experimentally accessible time window, FD curves x(C) H103|l 
are time-dependent and quite far from the expected asymptotic curve. This may 
explain previously reported discrepancies among FD plots obtained by analyzing the 
magnetization data of several spin glass systems |341| . These indirect measurements 
of FDT violations do not show any clear evidence of a universal curve x(C) and 
suggest that experimental measurements are quite far from the asymptotic regime. 
Experimental FD plots are shown in figure 03 

In structural glasses recent measurements show also the existence of FDT 
violations, however the physical interpretation of these experiments is still unclear. 
Voltage noise measurements |342| in an electric resonant circuit formed by a capacitor 
containing glycerol and an inductance show that the Nyquist formula H278|) is 
violated depending on the waiting time and the frequency. By defining the effective 
temperature Tf,s{ijJ,ty^) as in (|133|l (i.e. the temperature that satisfies the Nyquist 
formula (|278|l ) experiments reveal that the capacitance ages and FDT violations 
appear also in the range uotw ^ 1. This seems to be in contradiction with previous 
old experimental results in spin glasses. The origin of this discrepancy is presently 
unclear. Other recent experiments by Ciliberto and coworkers 18431 1344] on Laponite, 
a synthetic clay of charged particles that have the shape of a disc, show strong FD 
violations. If let evolve inside water, Laponite generates a colloid glass consisting 
of a packed irregular structure of disks due to the complex pattern of quadrupolar 
interactions. It is found that the value of the effective temperature Tcs{L0^tw)^ for 
low enough frequencies, is up to 3 or 4 orders of magnitude larger than the bath 
temperature. Although this has been interpreted as an evidence for coarsening 
dynamics (where T^s oo) a conclusive explanation of the origin of these high values 
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is still unknown. One interesting aspect of these experiments is the scaling behavior 
ujt^^"^ observed in the power spectrum indicative that FDT violations persist even in 
the regime Lotyj :» 1, in agreement with the previous results on glycerol. Contrarily to 
what would be expected, noise measurements in a rheological experiment for Laponite 
do not detect significant FDT violations. 

Experiments that successfully clearly demonstrate the existence of effective 
temperatures related to FDT violations are certainly needed. These preliminary 
account of results shows that still much work has to be done in aging, driven or 
granular systems to provide a safe ground to many of these ideas. 

9. Conclusions 

In this review we have presented an overview on what is now an active area of research, 
i.e. the study of FDT violations in glassy systems. Glassy systems are widespread in 
nature and found in many different areas covering physics, chemistry or biology. These 
are non-equilibrium systems which are either in a non-stationary slowly relaxing aging 
state or in a weakly driven stationary state. Equilibrium systems arc often described 
by a set of intensive parameters such as temperature, pressure or density. In a similar 
way, in glassy systems an important role is played by the waiting time (aging systems) 
or the intensity of the driving force and/or its frequency (driven systems). These 
parameters describe how far from equilibrium the system is. In fact, it is becoming 
steadily clear that a thermodynamic description of glassy systems can be partially 
rescued. In this description the glassy regime can be rationalized by using some of the 
concepts of equilibrium statistical physics such as the existence of a modified version 
of the FDT, the so called quasi-FDT (QFDT). 

Associated to the existence of the QFDT there is the concept of effective 
temperature. Rather than being a useful parameter to describe the behavior of non- 
equilibrium systems, the effective temperature has a deeper physical meaning. It 
could be tli(^ vc^stige of the existence of some dynamical measure underlying the non- 
equilibrium regime. However we do not know how to prove the existence of the QFDT 
from first principles in the same way we do not know how to prove Boltzmann equal 
probability hypothesis in equilibrium theory. Therefore, establishing the existence of 
this dynamical measure is somehow equivalent to assuming the existence of a QFDT. 
Then the crucial point is which are the specific predictions, that can be experimentally 
tested, one can derive from the existence of a QFDT. The measurement of the effective 
temperature itself appears as the most direct way of challenging the QFDT. However, 
direct measurements of this non-equilibrium temperature appear quite difficult, results 
arc still very preliminary and much progress is yet to be done to reach a convincing 
and clear evidence. 

The study of several families of models, as described in this review, appears as a 
very fruitful source of inspiration for new concepts and ideas that could be eventually 
exported to different classes of problems. However this path is not free of challenges 
and ambiguities. Indeed what are natural concepts for some families of models appear 
to be quite artificial in other families. As an example, key concepts in a thermodynamic 
formulation of the glassy state are the existence of an effective temperature associated 
to the configurational entropy or complexity, in the same way the bath temperature 
is associated to the Boltzmann entropy in equilibrium theory. However, what appears 
to be an interesting quantity describing glassy systems with two-timescales cannot be 
easily translated into systems with many timescales. Even more, what appears to be a 
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meaningful quantity for models having a complex thermodynamics (such as spin-glass 
models) appears to be meaningless in models described by a trivial Hamiltonian (such 
as kinetically constrained models). 

From the point of view of theoretical studies, our understanding of the existence 
of a QFDT appears also quite problematic. Most of the solutions we have described 
in this review are only valid in the case of mean- field interactions, however their 
validity beyond that limit remains speculative. The use of numerical simulations has 
aided to bridge the gap. It is quite interesting that most of the results predicted 
in mean-field models are qualitatively also observed in short-range systems. This 
tendency to rationalize the behavior of real systems within a mean-field scenario, i.e. 
far beyond their natural domain of applicability, has become quite standard in the 
study of glassy systems. The description of the equilibrium properties of spin-glass 
systems has for many years followed a similar route. Although the knowledge we can 
gain from numerical simulations in glassy systems is always quite qualitative (cither 
due to the limited range of timescales or to the inherent simplicity of the simulated 
model) the accumulated evidence, as reported along this review, points toward the 
emergence of a QFDT in the non-equilibrium regime of glassy systems, reminiscent 
of how the Boltzmann measure emerges in equilibrium systems. The same conclusion 
holds for the Edwards measure in granular media. 

What will be the future in the research of FDT violations? Although modeling 
promises to offer new ideas and will clarify our understanding we feel that more 
progress is certainly needed in basic theory and experiment. In theory we need 
to imdcrstand the origin of the existence of a QFDT right from the microscopies. 
Cooperative processes in glasses involve a few tens of atoms and occur along 
nanometric length-scales. How to link the microscopic activated processes to the 
emergence of macroscopic properties (such as the effective temperature) is a real 
challenge. In this context, it appears quite interesting to pursue the investigation of the 
so called fluctuation theorems recently proposed to quantify transient violations of the 
second law of thermodynamics. From the experimental side the current accumulated 
knowledge is still too poor and more experiments are certainly needed for this field 
of research to grow. A future line of progress is the use of nanotechnology devices 
to do noise measiircmcnts over spatially localized regions of nanometric sizes. These 
devices could be used as a microscope to measure activated processes occurring in 
small length scales. 

Certainly we will sec upcoming developments in this exciting area of research. 
A continuous exchange of ideas among theory, simulation and experiments is highly 
desirable and certainly needed to improve our current understanding in this field. 
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10. List of abbreviations 

IRSB One-step Replica Symmetry Breaking 

BG Backgammon 

BTM Bouchaud Trap Model 

EA Edwards- Anderson model 

FD Fluctuation-Dissipation 

FDR Fluctuation-Dissipation Ratio 

FDT Fluctuation-Dissipation Theorem 

FILG Frustrated Ising Lattice Gas 

IRF Integrated Response Function 

IS Inherent Structure 

LJ Lennard-Jones 

ME Master Equation 

MF Mean Field 

MCT Mode-Coupling Theory 

OSC Oscillator 

QFDT Quasi-Fluctuation-Dissipation Theorem 

REM Random Energy Model 

RFIM Random Field Ising Model 

ROM Random Orthogonal Model 

RSB Replica Symmetry Breaking 

SK Sherrington-Kirkpatrick 

TAP Thouless-Anderson-Palmer 

Tc Mode-Coupling critical temperature, Dynamical critical temperature 

Tcff Effective temperature 

Tg Glass transition temperature 

Tk Kauzmann temperature 

Trsb Static transition temperature in replica calculations 
TRM Thermoremanent Magnetization 
TTI Time Translation Invariance 
ZFC Zero Field Cooled 
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